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PREFACE 


Tuis book is a translation of a review by Academician B. B. Kadomtsev of 
the I. V. Kurchatov Institute, Moscow, which was written at my suggestion 
in order to have available in concise form the essence of the large body of 
papers on this difficult yet important subject. It was originally published 
in Russian in Volume 4 of “ Problems in Plasma Theory”, edited by 
M. A. Leontovich (1964). 

Academician Kadomtsev is a leading authority on plasma turbulence, 
and he and his colleagues have contributed many of the original papers on 
the subject. It is now well known that some form of turbulence is very 
frequently and often disastrously present in experiments on the confinement 
of hot plasma. It is not impossible that in the quest for a thermonuclear 
plasma we shall ultimately have to deal with a more or less turbulent state, 
and the understanding of this state is therefore important to nuclear fusion 
research. I believe this book will contribute greatly towards such an under- 
standing and help to stimulate further work on the subject. 

The translation and editing of the book have been carried out by Mr. L. C. 
Ronson and Dr. M. G. Rusbridge of the Culham Laboratory. I am grateful 
to them and to Mrs. Mary Hardaker who typed the manuscript. 


J. B. ADAMS 
April 1965 


United Kingdom Atomic Energy Authority 
Culham Laboratory 
Abingdon, Berkshire 
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NOTATION 


THE notation conforms as closely as possible to that of the original, which 
is in general conventional or explained in the text. The only major exception 
is the substitution of curl for rot, which is not so familiar to English readers. 
Vectors are shown by bold type. The vector product is denoted by square 
brackets, thus [hk]. No special notation is used in general for the scalar 
product, but round brackets, thus (kk’), may be used if required for clarity. 
Otherwise two adjacent vectors imply the scalar product. Note the appearance 
of the scalar triple product in equations such as (IV. 128) and (IV. 129). 
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INTRODUCTION 


IT is well known that a real plasma is rarely quiescent; as a rule many forms 
of noise and oscillation arise spontaneously in the plasma. Langmuir 
pointed out that these fluctuations represent more than just harmless oscil- 
lations about an equilibrium position and often wholly determine the charac- 
ter of the phenomena occurring in the plasma. In particular, interaction 
between electrons and oscillations accounts for the “ strong ” scattering of 
electrons in a gas discharge, first described by Langmuir (1, 2). The collec- 
tive interaction of particles also accounts for the well-known Langmuir 
paradox; namely, that even at very low gas pressure the velocity distribution 
of the electrons in a glow discharge is, to a high degree of accuracy, a Max- 
wellian distribution. The decisive part played by oscillations in the inter- 
action of an electron beam with a plasma has been demonstrated in the 
papers of Merrill and Webb (3) and Looney and Brown (4). 

Experiments with plasmas in a magnetic field and in particular experi- 
ments on magnetic containment of a high temperature plasma in connection 
with controlled thermonuclear reactions have revealed further unexpected 
phenomena essentially connected with oscillations in the plasma. Prominent 
amongst these is the “‘ anomalous ” diffusion of a plasma across a magnetic 
field. This effect was first observed by Bohm, Burhop and others while 
investigating the operation of ion sources (5) and later the enhanced diffusion 
of a plasma, related to its instability, was observed in a series of experimental 
devices. 

Following the work of Bohm (5), who suggested that the enhanced dif- 
fusion of a plasma is due to random oscillations of the electric field set up 
by an instability, the term “ turbulence ”? has been increasingly applied to 
this process. At present we understand by turbulence the motion of a plasma 
in which a large number of collective degrees of freedom are excited. Thus, 
when applying the term “ turbulence ” to a plasma, it is used in a broader 
sense than in conventional hydrodynamics. If hydrodynamic turbulence 
represents a system made up of a large number of mutually interacting 
eddies, then in a plasma we have together with the eddies (or instead of them), 
also the possible excitation of a great variety of oscillations. Depending on 
the degree of freedom which is excited, the character of the interaction 
between the excitations may vary considerably. 

During the eddy motion of an ordinary fluid the separate eddies, in the 
absence of their mutual interaction, do not propagate in space. When their 
interaction is included the eddies ‘‘ spread out ” in space with time, though 


the corresponding velocity is not large and therefore each separate eddy has a 
I 1* 
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considerable time available to interact with its neighbours. In this case we 
are faced with a strong interaction of excitations and correspondingly with 
a strong turbulence. On the other hand, during a wave motion the separate 
wave packets can separate from one another over large distances. In this 
case the interaction of separate wave packets with one another is weak, and 
we can therefore refer to a weak turbulence. The motion of the plasma in the 
weakly turbulent state, constituting a system of weakly correlated waves, 
shows greater similarity to the motion of the wavy surface of the sea or the 
oscillations of a crystal lattice than to the turbulent motion of an ordinary 
fluid. 

The theoretical consideration of a weakly turbulent state is considerably 
facilitated by the possibility of applying perturbation theory, i.e. an expan- 
sion in terms of a small parameter such as the ratio between the energy of 
interaction between the waves and their total energy. The problem of the non- 
linear interaction between waves in a plasma has been considered by Sturrock 
(6) for the special case of interaction between Langmuir waves in a cold 
homogeneous plasma. 

For the case of very small amplitude, when the interaction between the 
oscillations can be neglected, one can use the so-called quasi-linear approxi- 
mation in which only the reaction of the oscillations on the average velocity 
distribution function of the particles is considered. The quasi-linear approxi- 
mation was referred to by Romanov and Filippov (7) and has been further 
developed in papers by Vedenov et al. (8), as well as by Drummond and Pines 
(9). Section I of the present review is devoted to the quasi-linear approxima- 
tion. 

Unfortunately the quasi-linear method has only a fairly narrow field of 
application, since non-linear interaction of the oscillations already begins to 
play a considerable part at not very large amplitudes. In the paper by 
Camac et al. (10), the non-linear interaction of Alfvén- and magneto-sonic 
waves is described by the kinetic wave equation, which is well known in 
solid state theory (11). Camac et al. (10) applied this method to describe the 
structure of a collisionless shock wave. The problem has been considered in 
somewhat greater detail by Galeev and Karpman (13), while the interaction 
of Langmuir waves has been studied in reference (19). 

In the simplest variant of the kinetic wave equation only three-wave 
processes are considered, namely the decay of the wave k, w into two waves 
k’, œw and k”, œ”, and the merging of two waves into one. Such processes 
are important only for dispersion relations @, = œ (k) for which it is possible 
to satisfy simultaneously the laws of conservation of energy and of momen- 
tum: k” = k—k’, a, = @,—a@,. When these conditions are not satisfied, 
scattering of the waves by the particles is a more important process and can 
be taken into account only on the basis of a full kinetic theory. Such a 
theory has been published by Kadomtsev and Petviashvili (14), who ob- 
tained a kinetic wave equation including the thermal motion of the particles. 
This equation was obtained later by a slightly different method by Karpman 
(15). The kinetic wave equation is deduced and considered in Section II. 
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Unfortunately, in numerous practical cases one is faced not by weak 
but by strong turbulence. In particular, strong turbulence is related to an 
anomalous diffusion of the plasma across the magnetic field. To determine 
the fluctuation spectrum in a strongly turbulent plasma and the effect of these 
fluctuations on the averaged quantities, it is sometimes possible to use the 
analogy with ordinary hydrodynamics and, in particular, to apply a pheno- 
menological description of the turbulent motion. 

Such an approach has been used by this author in two specific problems, 
referring to turbulent diffusion of a plasma in a trap with magnetic mirrors 
(16) and in the positive column of a glow discharge (17). In these cases the 
concept of the mixing length was used. The results of this discussion are in 
good agreement with experimental data, which indicates the success of this 
approach. 

However, in a plasma other strongly turbulent motions which are dif- 
ferent from the eddy motion of an ordinary fluid may develop. It is therefore 
desirable to have available more systematic methods for describing strong 
turbulence. In our view, such a method may be the weak coupling approxima- 
tion discussed in Chapter III. In this approximation, which ought preferably 
to be called the intermediate coupling approximation, the turbulent motion 
is described by a system of non-linear integral equations for the spectral 
density Jke and the Green’s function Gko describing the response of the 
system to an external force. As the coupling between the oscillations 
decreases, this system of equations goes over into the kinetic wave equation. 

In conventional hydrodynamics, the weak coupling equations have been 
obtained by Kraichnan (18) who showed that in their simplest form the weak 
coupling equations lead to a spectrum which is different from Kolmogorov’s 
spectrum in the region of large k. As will be shown in Section III.2, the 
reason is that in Kraichnan’s equations the adiabatic character of the 
interaction of the short wave with the long wave pulsations is not 
taken into account. The consideration of this adiabatic interaction makes it 
possible to obtain improved weak coupling equations. 

In Section IV, specific examples of turbulent processes in a plasma are 
considered. In particular the interaction of Langmuir waves and the excita- 
tion of ion oscillations by an electron current are considered, but the main 
attention is devoted to the turbulent diffusion of a plasma in a magnetic 
field. 

As we have mentioned earlier, the turbulent diffusion problem goes back 
to Bohm (5), who put forward the hypothesis that an inhomogeneous 
plasma in a magnetic field must always be unstable because of the presence 
of a drift current of the electrons relative to the ions. If this be in fact so, 
the corresponding instability must lead to a turbulent ejection of the plasma 
with a velocity of the order of the drift velocity. According to Bohm, this 
process can be considered phenomenologically as a diffusion with coefficient 
of diffusion of the order 10¢T 


e= H 
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where T is the electron temperature in electron volts and H the magnetic 
field in kilogauss. 

Bohm’s argument gave rise to the illusion of a universal validity for this 
coefficient and as a result attempts to obtain Bohm’s coefficient from more 
general considerations have continued to this day. It has now become 
evident, however, that the coefficient of turbulent diffusion cannot be ob- 
tained without a detailed investigation of the instability of an inhomo- 
geneous plasma and in particular of its drift instability. 

Investigations of the drift instability of a plasma in a magnetic field were 
started by Tserkovnikov (19) who limited his investigations to perturbations 
constant along the direction of the magnetic field. Rudakov and Sagdeev (20) 
went a step further by considering instabilities at oblique angles with a trans- 
verse wavelength considerably larger than the mean Larmor radius of the 
ions. The complete investigation of the drift instability has been carried out 
only quite recently, following the work of Rosenbluth, Krall and Rostoker 
(21), with the investigation of perturbations with transverse wavelength of 
the same order as the mean ion Larmor radius. The principal results in this 
field were obtained by Mikhailovskii (22). In a dense plasma where an 
important part is played by collisions between the particles, the drift instability 
changes to a drift-dissipative instability first observed by Timofeev (23). Ina 
note by Sagdeev and Moiseev (24), it was stated that this drift-dissipative 
instability may lead to an escape of plasma of the order of the Bohm 
diffusion. 

The principal results of the investigations of the drift instability of a 
plasma are discussed in Section IV.3. The following section considers the 
problem of the type of diffusion resulting from the drift instability in specific 
conditions. 

In Section IV.5, various turbulent processes in specific experimental 
devices are described and briefly discussed. Turbulent processes are con- 
sidered in toroidal discharges and in magnetic traps and brief reference is 
also made to experimental data on turbulent heating and diffusion of plasmas. 

In writing this review, we have aimed at maximum clarity of expression; 
where rigour appears to conflict with simplicity, simplicity is given preference. 
Attention has been given chiefly to processes occurring under laboratory 
conditions; problems of astrophysical application are completely omitted. 
In particular no reference is made in the present review to a broad group of 
investigations in which the turbulence of an ideal conducting fluid (i.e. 
magnetohydrodynamic turbulence) and turbulent shock waves are con- 
sidered. 


|. A QUASI-LINEAR APPROXIMATION 


1. INSTABILITY AND TURBULENCE 


As a rule, turbulence develops as a result of an instability of an initial laminar 
state. In order to visualise the transition from the laminar to the turbulent 
State, it is convenient to examine the behaviour of the system while changing 
some parameter, R, the increase of which results in a loss of stability. In 
conventional hydrodynamics such a parameter is the Reynolds number, 
whilst in a plasma there are a number of cases where the magnetic field 
strength plays a similar role. 

Let us recall first what occurs during the loss of stability of a system 
with one degree of freedom, for instance, a valve oscillator, where the part 
of the parameter R is played by the feed-back. If the feed-back R is smaller 
than some critical value R, then all small oscillations are damped. When 
R > Ra, on the other hand, the state with zero amplitude is unstable and 
oscillations of finite amplitude are excited. In other words, the value R = R, 
is a bifurcation point above which the oscillator is in an excited state. 

In these conditions two types of excitations are possible—a soft type and a 
hard type. In the soft regime the amplitude of the oscillation varies continu- 
ously with R, vanishing for R = R, (see Fig. 1). In the hard regime the ampl-, 


Ae 


R 


c 


Fic. 1. Soft excitation 


tude increases abruptly to some finite value as soon as the value R cxceeds 

the value R, and when R decreases it drops abruptly to zero for R = Ro < R. 

(see Fig. 2). In the region Rọ < R < R, the circuit is stable with respect to 

infinitely small perturbations, but unstable to perturbations of sufficiently 

large amplitude. The unstable equilibrium state is represented in Fig. 2 by the 
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Fic. 2. Hard excitation 


dashed line. In the case of a small non-linearity of the system this picture of 
the excitation is fully described by means of the widely known Van de Pol 
method, i.e. a small amplitude expansion. 

Landau (25) has shown that the transition to the turbulent state also 
constitutes an excitation of the system as a result of an instability. To 
examine Landau’s argument let us consider a continuous medium, 1.e. a 
system with an infinite number of degrees of freedom, and assume that the 
excitation is soft. Then, as R is increased the following picture will be ob- 
served. For small R we have the laminar state in which all quantities are 
completely defined by the initial and boundary conditions, i.e. the system 
has no superfluous degrees of freedom. For R greater than some critical value 
R, a normal mode of the system is excited whose amplitude increases mono- 
tonically with R; in other words, an additional degree of freedom appears 
in the system. As R is further increased further degrees of freedom may be 
excited and ultimately we arrive at a turbulent state. 

When R only slightly exceeds R,, the amplitude A can be determined by 
expanding with respect to R—R,, similarly to Van de Pol’s method. This 
method was first applied to the problem of thermal convection of a fluid 


Ae 


R; R2 Rx 


Fic. 3. Mild excitation of a system with many degrees of freedom 
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(Sorokin (26)) and was later used by Stuart (27) to describe the eddy cells in a 
fluid between rotating cylinders. This method is now known as the quasi- 
linear method and leads to satisfactory results for R— R4 < R4. 

The behaviour of the system as R is increased also depends on its specific 
properties. We may find that numerous strongly interacting modes are 
excited and in this case the transition to turbulence takes place fairly rapidly. 
Alternatively we may find that all the higher modes are harmonics of the 
fundamental mode first excited and are synchronised with it. In other 
words, we have a non-linear oscillation of finite amplitude which really 
represents only one degree of freedom. (The losses in the positive column of 
a glow discharge are an example of this situation.) Finally, there is the pos- 
sibility that the excited modes will interact only weakly with one another 
and a weakly turbulent motion will develop. If the interaction between the 
modes can be neglected, weak turbulence can be described by the quasi- 
linear approximation. 


2. LAMINAR CONVECTION OF A PLASMA 


Let us consider two simple examples where the instability of a plasma 
leads to the appearance of a convective flow. 


(a) Convection of a Weakly Ionized Plasma in an Inhomogeneous Magnetic 
Field 


As a first example, let us consider the convection of a weakly ionized 
plasma in an azimuthal (toroidal) magnetic field (28). As is well known, a 
fully ionized plasma in such a field is convectively unstable; because of its 
diamagnetism it is pushed out radially. If the magnetic field is sufficiently 
high, then a weakly ionized plasma is subject to a similar instability. Let us 
consider here what happens to a plasma when the magnetic field exceeds the 
critical value H, at which the instability first appears. 

Suppose that the plasma is located between two ideally conducting 
cylinders of radius R and R+d respectively, (d < R), and ionization is 
achieved in such a way (for instance using heated grids in caesium vapour) 
that a constant density n is maintained at the inner cylinder and a density 
n—O0n at the outer cylinder, where ôn <n. 

We shall assume both the electrons and the ions are “ magnetized ”, 1.e. 


eH 
Qt; =——1; > 1, where q; is the mean free time and m; the mass of particles 
mC 


of type j. Then from the equation of motion for the electrons 
en mn 
T.Vn = enVo — yeH] aN (I. 1) 


e 


we obtain the following 


T, 
nv, = = |h, nVo + eyn] — D,Vn+nb,V@ (I. 2) 
e 
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D T, ee l 
: - gz is the transverse diffusion coefficient, g the 
e 


where Di = ->> = —— 7 
Q772 Mete 


eD 
potential of the electric field, T, the electron temperature and b, = a the 


transverse mobility of the electrons. Since we assume Qt, > I, the last 
term in I.2 which contains b, can be neglected. 

Substituting I.2 into the continuity equation and neglecting the curvature 
of the magnetic field, we obtain 


0 
A + vVn = D, An (I. 3) 


where v = = (hV¢) is the drift velocity in crossed electric fields and h = H/F. 


On the other hand, from the ion equation of motion, assuming 7; = 0, we 
have 
c eT; 
v; = —[hVo] — ———“ V I. 4) 
l H [ o] m;(Q; w) p ( 
We now substitute the expressions for nv, and nv; in the quasi-neutrality 
condition div n(v;—v.) = 0. Considering that 


g Onl a (r) 2 an 
HH ézrér\H) = RH @z’ 
we obtain 
c 2D, ôn 
— di = — ——— I.5 
7 div (nVo) R os (1. 5) 
where D, = 5 is the ambipolar diffusion coefficient in the absence of a 


i 

magnetic field. 
In the equilibrium state n = n (r), ọ = ọ (r), and consequently there are 
no convective fluxes. Let us consider the value of the critical magnetic field 
H, at which convection starts. Since ôn < n, the density n can be assumed 
constant on the left hand side of relation I.5 and the perturbation of the 
TX 


density and potential can be chosen in the form sin — exp (ikz— iœt) where 


d 
x = r—R. The linearised eqns. I.3 and I.5 then take the form 
(1) 
; co’ dn 
— iw + D (k? +k) Jn ® = — ik — — I. 6 
[-i@+D,(k o)|n l H dx (I. 6) 
Ch ey2 4 LA (1) o Teti a) 
—(k*+k = — ik —— ; 
H' EOP i m, Rn” (Sa 
where kọ = n/d. From this we obtain 
2D, k? On 


—iw =—D,(k*+k~) +—* s 
AERE pa ka id 
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2D, k? ôn 
R (k (k?+ k2)? nd 
maximum value of D, is attained for k = kọ and is equal to 


_ Dad on (1. 8 
€ WR n 8) 


Clearly D, < D,. The critical field H, is defined by the relation D = D,; 

for H > H,, i.e. D < D,, laminar convection develops in the plasma. 
When H is only slightly larger than H, the convective flow can be deter- 

mined using the quasi-linear method. The total diffusion flux is defined by the 


For a = 0, ie. D, = the plasma becomes unstable. The 


; dn 
relation q=qọ+ôq, where go = D, ie ôq = <n and the angular 


brackets indicate averaging with respect to time. In the first approximation 
v can be expressed in terms of n“ with the aid of (1.7) and for n“ we 
insert the dependence on x and z given above. Integrating the resulting 
expression for the flux with respect to x with q = const gives 
1D, nd , 


D ôn + =—— d i 
1ôn +5 4A =a (I. 9) 


where 4 is the amplitude of the density oscillations defined by 
(1) 
n TOG E N 
See E A sin — e'*?7 iat 
n d 
A second relation connecting q and A? is obtained by substituting go? 
obtained from (1.7) into (1.6), setting œ = 0, and multiplying the result by 
n“), Representing the derivative of the density in the form 


dn q A Cie) 
dx Dy, D, 
averaging the result with respect to x, and using (1.8), we obtain 
D? 3 D nd 
— ôn —“— A’? = I. 10 
D.” TI R ae re 
From the relations (1.9) and (1.10) we obtain 
ôn\? 2 D 
Aes eed ee I. 11 
ye g 
26nD , D, 
= f= I. 12 
a= ao + 7+ (1- 5) (1.12) 


Thus, the flux q depends on the magnetic field as shown in Fig. 4. For 
H > H, the flux increases instead of decreasing with increasing field. 

This picture is completely analogous to the picture of the development 
of convective cells in thermal convection between two parallel plates; here 
the onset of convection is also sharp. In more realistic conditions of a chamber 
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Fic. 4. Plasma flux against magnetic field during development of convection 


bounded along the magnetic field, there may be no well-defined critical field, 
the flux g varying smoothly with the field. 

The increase of the effective diffusive loss due to convection has been 
clearly observed by Bostick and Levine (29) and has been investigated more 
recently by Golant et al. (30), who investigated the decay of a plasma in a 
toroidal magnetic field. The plasma was contained in a toroidal glass vessel 
of major radius R with a circular cross section of radius a. A fairly long inter- 
val of time was allowed for the density distribution to become established 
and the decay of the plasma density was then observed and found to follow 
an exponential law n = nye “^. 

The loss can easily be determined by assuming approximately that the 


. @ 
convective velocity parallel to r can be neglected, 1.e. T =~ 0. We then 


e, and eqn. (1.3) takes the form 


a 


obtain from (1.5) v = R 


on j 2D, 
ot R 
where e is the unit vector parallel to the x axis. 
Neglecting the field curvature in the expression for the Laplacian, eqn. 
(1.13) can be considered as the equation for diffusion in a straight cylindrical 
tube in the presence of a transverse flux. It is easily seen that the solution of 


eqn. (1.13), with the boundary condition that the density at the wall vanishes, 
is given by the expression 


a t vo XoF 
n = Ng exp EuT Jo P (1. 14) 
i 


where x is now measured from the tube axis. In (1.14) my = const, a) = 2°4 
is the first root of the zero order Bessel function Jọ and the decay constant is 


given by the relation 
cone (+r L. 15 
t 2to\Do D, Ge) 


eVn = D, An (1. 13) 
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where 
a av 2D 


Se ee —= 


9 0 = V 
Oe v 205 R 


According to (1.15), the decay constant as a function of the magnetic 
field must pass through a minimum at D, = Dy. In Fig. 5, which was 


To = 


0-2 04 05 08 1 2 4 


D 
Do 


Fic. 5. Lifetime of a decaying plasma in a toroidal tube 
(helium, R = 28 cm, a = 0:3 cm) 


@ p= 0-025 mm Hg 
O p= 0-04 mm Hg 
x p= 0055mm Hg 
O p=012 mm Hg 


taken from (29), a comparison is given of the experimentally determined 
relationship between t,/t and D,/Do (where 1/tọ and Do are the decay 
constant and the diffusion coefficient respectively at the minimum) with the 
theoretical formula (1.15). As may be seen, the experimental results are in 
good agreement with the theory. 


(b) Convection of the Plasma of the Positive Column in a Magnetic Field 


Let us now consider the convection of a plasma in a homogeneous 
magnetic field in the presence of a longitudinal current. The development of 
convection of a weakly ionized current-carrying plasma was described by 
Lehnert (31) while investigating the diffusion of charged particles from the 
positive column of a glow discharge located in a strong magnetic field. It has 
also been studied in greater detail in a paper by Lehnert and Hoh (32). 

In these investigations the relationship between a longitudinal electric 
field E and the magnetic fìeld H was studied. In the positive column a 
decrease of the diffusion coefficient leads to a decrease of the electric field, 
so that if the diffusion is classical, the electric field E should monotonically 
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decrease as the magnetic field is increased. This relationship between E and 
H is in fact observed experimentally, but only for not too strong magnetic 
fields. As soon as H exceeds a critical value H,, the dependence of E on H 
changes markedly: the electric field begins to increase, ultimately attaining 
a saturation value E, (Fig. 6). 


H 


Fic. 6. Dependence of longitudinal electric field on magnetic field in the positive column 
of a helium discharge: 1. Region of classical diffusion. 2. Region of laminar convection 
(helical discharge). 3. Region of turbulence 


The general nature of this effect was explained in a paper by the author 
and Nedospasov (32) where it was shown that the positive column of a glow 
discharge loses its stability for H > H,. The mechanism of this instability, 
which is known as the current-convective instability, can be illustrated by 
considering a plasma filament slightly curved along a helical line. To fix 
ideas suppose the discharge is vertical, then in the presence of a longitudinal 
electric field E, which we suppose directed upwards, the electron “ corkscrew ” 
of the curved filament is displaced downwards relative to the ion corkscrew 
and charges appear at the surface of the filament, namely a positive charge 
at the upper surface and a negative one at the lower surface. These charges 
set up an azimuthal electric field E, which leads to a drift of the plasma 
along the radius, i.e. to an increase of the initial perturbation. In a weak 
magnetic field this effect is suppressed by the diffusion. In a strong magnetic 
field, when the transverse diffusion decreases to such an extent that it cannot 
eliminate the density perturbations, this mechanism leads to instability of 
the plasma filament with respect to a helical distortion. 

To estimate the value of the critical field, let us consider the stability 
problem in the W.K.B. approximation; this is accurate for short wave 
perturbations only, but describes qualitatively correctly perturbations with 
a longer wavelength, of the order of the tube radius a. In this approximation 
jt is assumed that the density of the plasma is a slowly varying function of a 
co-ordinate x, and the perturbation is selected in the form of a plane wave 
exp (—iwt+ikr). At a neutral gas pressure of the order of 1 mm Hg, the 
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instability begins when 0,7; < 1, so that the ion velocity can be assumed to 
be given by v; = —b;Dọ. Substituting this expression into the continuity 
equation we obtain for the perturbation of the ion density the following 
expression: 

n; bi k? 


n @ 


p (I. 16) 


where @ is the perturbation of the electric field potential. 

The perturbation of the electron density n, can be determined from the 
continuity equation in which the electron velocity v, is given by the relation 
(1.2). In the W.K.B. approximation the linearised continuity equation for the 
electrons takes the form 


ik cd 
LN bel? + bik?) @ = 0 


H dx 
(I. 17) 


where b., D, are the longitudinal, and b,, D, the transverse mobility and 
diffusion coefficient respectively of the electrons, u = b,E the mean velocity 
of the electrons in the longitudinal electric field E, and E, the transverse 
electric field in the equilibrium state, which in the case of ambipolar diffusion 
is given by 


E 
(-io+ 2D, + 2D, + ik ik, Z2) ne — ( 


dl 
E, =—D,b; G 
dx 


b. 
( y= ziaz). The z axis is parallel to the magnetic field and the x axis to 


the direction of decreasing density. Using the quasi-neutrality condition 
it, =n; we can substitute (I.16) for n,/n in (1.17), obtaining the following 
dispersion equation 
b; k?Do+kł D, +ik,u— ikk, D (1 +y) (Qet) + (L. 18) 
= nin Aaaa : 
O= bdb, +Q t) 2+ K2/k2 — ik, K/Q, tek? 


dl a ie 
where K = -S According to (I.18) the instability condition Im œ > 0 


can be represented in the following form: 


Kk? b, uk,k 
X*4+(2+y)X*+1+y+—[ l 


k* lty b,” Dk? 
where X = k,Q,1,/k,. 

According to (1.19) the instability occurs only in the presence of a longl- 
tudinal current, u 4 0. Since the right hand side of (1.19) increases when k 
decreases, a perturbation must develop in the first place with the minimum 
possible k, i.e. k ~ x. The corresponding perturbation has the form of a 
curvature of the filament as a whole. Since usually b,/b; ~ 107—10° > 1 
and ux/D, ~ 1, the instability may arise even for y < 1. In this case X ~ 1, 
ie. k, ~ K(Q,2,)71 < «. In other words the longitudinal wavelength of the 
perturbation is considerably larger than the tube radius. 


x (I. 19) 
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All these qualitative conclusions correspond adequately to the experi- 
mental data and more exact calculations lead to quantitative agreement 
with the experiment also. Figure 7, for instance, shows a comparison be- 
tween the theoretical relationship between H,/p and ap (a being the tube 


Š 


ap 


Fic. 7. Comparison of theoretical and experimental dependence of critical 
magnetic field on neutral gas pressure, discharge in helium 


radius, p being neutral gas pressure) and experimental data of Lehnert and 
Hoh (32). The appearance of a helical curvature of the filament for H > H, 
was demonstrated directly by Allen et al. (34) using streak photography. 
Paulikas and Pyle (35) have published a detailed investigation of the transition 
through the critical field. Figure 8 shows their experimental relationship 
between the pitch of the twisted filament and the frequency of its rotation, 


Fic. 8. Comparison of theoretical and experimental dependence of pitch of helical 
discharge A and frequency of rotation œw on neutral gas pressure (discharge in helium) 
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and the neutral gas pressure for a discharge in helium. The theoretical 
relationship is shown by the full lines. As may be seen, the agreement 
between theory and experiment is very satisfactory. 

The helical distortion of the filament for a magnetic field slightly exceeding 
H, is simply the result of the convection of the plasma. For small H— H, it 
is again possible to use the quasi-linear approximation to determine the 
convective loss of the plasma and the amplitude of the distortion. Such a 
calculation (33) shows that the outward flux of plasma increases with H— H, 


and the amplitude increases as JH — H.. The curved section 2 on Fig. 6 
represents the dependence of E on H calculated from this theory. As may 
be seen, for a magnetic field only slightly greater than the critical field, the 
theoretical and experimental dependence agree very well. For H > H, the 


positive column goes over into a turbulent state which will be considered in 
Section (1V.4). 


3. QUASI-LINEAR APPROXIMATION IN KINETICS 


We now turn to the opposite extreme of a collisionless plasma and con- 
sider first the application of the quasi-linear method to the excitation of 
Langmuir oscillations by an electron beam (7) (8) (9). 


(a) Electron Beam in a Plasma 


Suppose that the electron distribution function has the form shown in 
Fig. 9 by the full line. The second maximum represents a diffuse beam of 


Fic. 9. Initial (full curve) and final (dashed curve) distribution functions 
during excitation of Langmuir oscillations 


electrons of mean velocity V, superimposed on the main group of thermal 
electrons. | . 
As we know, such a velocity distribution is unstable: Langmuir waves 


d, A PEO 
with phase velocities in the region where A > 0 will increase with time, 


since the number of electrons overtaking a wave and transferring energy to it 
will be larger than the number of electrons withdrawing energy from it. 
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Thus, as a result of Landau’s inverse damping mechanism a group of waves 
will be built up, the phase velocities of which cover the whole interval where 
df 

We 0. 

If the number of resonance particles is small, i.e. the particle density in 
the beam is considerably smaller than the density of the thermal electrons, 
the growth rate of the waves y will be considerably smaller than the frequency 
w. For y/w < 1 the interaction between the waves can be neglected and we 
can use the quasi-linear approximation, in which the only non-linear effect 
considered is the reaction of the oscillations on the ‘‘ background ’’—the 
averaged distribution function. 

In the unidimensional case the motion of the electrons will be described 
by the following system of equations: 


OF F eEOF 
ee (I. 20) 


aa = 
Ot Ox mv 


r= = 4ne} | F dv—n} (I. 21) 
where n is the ion density and F the electron distribution function. 

Let us split up the distribution function into two parts: F = f+f 
where fis the mean distribution, regarded as a slowly varying function of f, 
and f is the oscillating part which averages out to zero. The function f® 
represents a system of oscillations with randomly distributed phases 


FMS (je dk (I. 22) 
where œ, is the characteristic frequency of the k*® mode. For the Langmuir 
‘lat e : oe nw) Of 
oscillations œ, œ= w , where wo = , | —, Y; = —— . 
a A 2 m T 2n av Sere 


From the linearised form of eqn. (1.20), we determine the relation be- 
tween fy and the amplitude of the electric field E, in the k mode 
e i of 
= — ——— >E, +A,6(@,—k ' 
fr ee p+ A, 0(@, — kv) (1. 23) 
Here the second term allows for the possibility of the injection of weakly 
modulated beams into the plasma. For simplicity we assume that such 
beams are absent, i.e. A, = 0. 
The equation for the averaged function fis obtained by writing F = f+f" 
in equation (1.20) and expressing f? in terms of E, from (1.22) and (1.23). 
Averaging with respect to the statistical ensemble, i.e. with respect to the 


d 
random phase, and noting that a = 0, Ey = E_,, we obtain 


S2 (yO 
at dv Q a (en) 
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where 
YE) 
~ m? ) (@,—kv)?+y2 
For a constant mean distribution function the dependence of E? ont 


would be simply exponential, EZ ~ e?”*', but for a weak dependence of f on 
time Ex is given by 


(1.25) 


D, 


2 
ae = 2y,E? (I. 26) 
dt 

Equations (1.24) and (1.26) are the basic equations of quasi-linear kinetic 
theory. 

With the aid of these equations we can now examine the development of 
unstable oscillations of a plasma. Let us consider first the group of resonance 
electrons. Since y < w, for thts group of particles the coefficient of diffusion 
in velocity space D, can be written sufficiently accurately in the following 
form 


re? 
D; = =| 6(w, — kv) E? dk (I. 27) 


This coefficient of diffusion is different from zero only in the region of 


velocities where y > Q, i.e. z > 0. Because of the diffusion the distribution 
function will be flattened in this region until the growth rate vanishes, i.e. 
until a “ plateau ” appears in the distribution function. At this time a 
stationary spectrum of supra-thermal oscillations will be established in the 


plasma. To S the spectrum we use eqns. (1.24), (1.26), (1.27) and the 
w? 


o to obtain the following equation 


ô ô e*n 
AFIT E2\\=0 
ôt f+ dv (Fs 2) 


a 2 
s+ see ; E2) = const 


dv \m*Wov 


where instead of k we use the variable v = w ./k. From this relation we 
determine the oscillation amplitude in the steady state 


relation y, = 


whence we obtain 


EX(t = 0) = ste o° | {f(t = 00) f(t = 0)} dv (I. 28) 


where f(t = œ) = const. In (1.28) the value v = v, represents the lower 
boundary of the region of velocities where the plateau is established. It is 
determined by the condition of the conservation of the total number of 


resonance particles, namely fjdv = const. 
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Let us now consider what happens to the non-resonance thermal electrons. 
For such electrons (ow—kv)* > y? and, consequently, the corresponding 
coefficient of diffusion, broadly speaking, is @/y times smaller than the 
coefficient of diffusion for the resonance particles. But since the number of 
thermal electrons is just w/y times greater than that of resonance electrons, 
their diffusion cannot be neglected. From (1.25) the corresponding diffusion 
coefficient can be represented en in the following form: 


Yk Ek 
E 1.2) 

Since the coefficient of Pra of the non-resonance particles is small, 
we can replace f by f(t = 0) in this region on the right hand side of eqn. 
(1.24). Considering (1.26), the integration of (1.24) with respect to time is 
straightforward and we obtain the following expression for the change in the 
distribution function of the — (non-resonance) electrons: 


3 4 1 ôf =| 
f(t = 0) — f(t = 0) = < | E= 0) =| 55, | dk 

Thus, the full distribution ne for t = œ will have the form repre- 
sented on Fig. 9 by the dotted graph. All the momentum and half of the 
energy lost by the beam in building up the oscillations are transferred to the 
thermal electrons, which leads to a distortion of their distribution function. 
The remainder of the energy lost by the beam is stored in thc electric field. 

We can show for instance that the total momentum of the beam and of 
the plasma are conserved. By multiplying eqn. (1.24) by mv and integrating 
with respect to v we obtain 


2 
4 
ral mofdv = — E f atwo- ko)ERS dk dv — : — {42 d 12 ak dv 
(@o— kv) ô 
where in the second integral the range of integration covers the region of the 
2T, 


thermal electrons |v| ~ v, = m 
, l ! , sd 1 2kv 
Since in this region (@ọ— kv pd me = a we can integrate with 
0 
respect to v and, substituting for y, its expression given above (following 
(1.22)), it can be seen that the integrals on the right hand side of this expression 


exactly compensate one another. Hence the total momentum is conserved. 


(b) Waves in a Plasma 


In the general case the problem of describing a weakly turbulent state 
in the quasi-linear approximation can be split into two parts: the determina- 
tion of the wave field and the consideration of its reaction back on the par- 
ticles. Let us consider the first of these problems in greater detail. Since in 
practice the plasma is almost always inhomogeneous, we must understand 
the propagation of waves in an inhomogeneous plasma. Since we have in 
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mind a turbulent plasma in which a large number of waves has been excited, 
the wavelength for the main part of the oscillations must be considerably 
smaller than the characteristic dimensions and consequently the inhomo- 
geneity can be considered weak. 

We recall that in a homogeneous unbounded plasma all characteristic 
oscillations are plane waves of the form exp (—/wft+ikr). Maxwell’s equa- 
tions for such oscillations have the following form: 


[kH] + < Eye =( (I. 30) 


[kKE,..] -~ Hro = 0 (I. 31) 


where é is the dielectric permeability tensor of the plasma. 
Expressing H in terms of E from (1.31), we obtain a single vector equation 
for E, which in Cartesian co-ordinates has the following form: 


2 
p 
so that the frequency w of the characteristic oscillations must satisfy the 
dispersion equation 
D(k, œw) = Det |k ap — 0c eng — ka ke|| = 0 (I. 33) 


In general the characteristic frequencies are complex: from the point of 
view of the quasi-linear approximation we are interested mainly in the special 
case in which the imaginary part of the frequency is small compared to the 
real part and is due to the interaction of the wave with resonance particles. 
The growth rate of the wave is determined by the anti-Hermitian part of the 


* 
e ° eqe yf Ez — E . e . 
dielectric permeability lesg ae Se which is small compared with the 
be , bap t baa 
Hermitian part £g = a ae 


If eg, is neglected the characteristic functions represent a plane wave of 
constant amplitude. When these waves are uniformly distributed in space 
and are statistically independent, the correlation function for the electric 
field has the following form 


CE (k, o)ER(K',@')y = Val al 1(kK)5(k—k')5(@—')(— Oj) (I. 34) 


where the summation is over all characteristic frequencies œx belonging 
to the wave vector k, a; is the unit polarisation vector (which js in general 
complex) defined by the relation E = aE, and I(k) is the spectral distribution 
of the electric field. 

In the presence of a magnetic field, several different characteristic fre- 
quencies correspond to one value of k in general. In this case all waves 
can be considered uncorrelated, in contrast to conventional optics where it is 
necessary to consider the correlation of the polarisations, and for a complete 
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description of two transverse oscillations it is necessary to introduce the four 
Stokes parameters. 

When the small anti-Hermitian part €e% is taken into consideration, 
the amplitude of the oscillations varies with time and the correlation function 
will not be simply proportional to 6(#—w’). But when considering a large 
number of waves in a wide frequency interval, it is still possible to use the 
relation (1.34) as an approximation, by allowing J (k) to be a slowly varying 
function of time. Differentiating this relation with respect to time gives 


—i(a—w')E(k, E(k’, 0’) = 
=) (a! aI (k, ))5(k—k)0(o—o')5(m—wi) (1.35) 


from which it is immediately apparent that the variation of J with time leads 
to a “ broadening ” of the ô function of w—w’ in the expression for the 
correlation function. 

A spatial inhomogeneity must lead to an analogous broadening of the 
function 6(k—k’), and this can be considered approximately by allowing 
a weak dependence of J on r. In this case the operator V corresponds to 
i(k— k’). 

The approximate representation of the correlation function in the 
form (1.34) corresponds physically to the concept that the electric field of the 
oscillations forms a set of statistically independent wave packets. If the 
characteristic wavelength of the oscillations is considerably smaller than the 
dimensions of these packets, then over limited intervals of time we can neglect 
their dispersive spreading. To obtain an equation describing the behaviour 
of such packets we use the following averaging operation: we multiply 
eqn. (1.30) by Eğ, eqn. (1.31) by Hý, subtract the one from the other 
and then anti-symmetrise the resulting relation, i.e. subtract from it an 
analogous complex conjugate relation with k, w and k’, œw transposed. 
In the relation so obtained the value of the quantity ej, (k, œw)— eag (k’, w’) 
can be expanded in series for small k—k’ and m—q’. Finally we insert 


0 l 
—iV andi ar respectively for k—k’ and w—w’. An analogous exchange must 


be carried out in the other components where corresponding differences 
appear. For a homogeneous plasma we thus obtain the following energy 
balance equation 


oW 
ot 


where W is the energy of the ktt wave given by 


+ divS = — ra Im (E*, cE, ,) (I. 36) 
7 


a at 0 KA 1 x 
a (wE? EY + J <H*H> 
1 0 e7k? c? 
— oe E*?’E eee, $ a $ 
87 dw Ko 2 — 7 SE E> + — (KE VKE) | (I. 37) 
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and S is the energy flux 


TEE = (we E) = 


c? 


4 — <E*E) + — = <(Ke") KE} (I. 38) 


ê’ is the Hermitian part of the tensor ê. In the expressions for W and S 
the differentiation with respect to w and k can, in view of (1.32), be carried 
out at constant E*, E. 


Multiplying eqn. (1.32) by E, and summing the result with respect to «, 
we obtain 


2 
GE -z7 EXE k?E*E + (kE*)(kE) = 0 (I. 39) 


and it is clear that W and S are given by the derivatives of the mean value 
of the function G, defined by 


<G> = a (a*ta)—k?+ [kat = D,I (I. 40) 


The quantities D, and D (see (1.33)) vanish at the same frequency o, 
so that at some other œ = œ, they can be supposed proportional to one 
another. We then obtain 


ImDi, y, f% __ oD, (2 j 
ôD, ’ k dk ôk \ do 
0@ 


y= (I. 41) 


where y, is the growth rate of the k'* wave, U, its group velocity. Using 
these relations we transform the energy balance eqn. (1.36) to the following 
simpler form 
aI ol 
A + U5 = 2y, 1, (I. 42) 
Let us now consider an inhomogeneous plasma. The inhomogeneity 
leads to several additional features. Most important, in the tensor £g some 
small additions appear which are proportional to the gradients of the mean 
distribution functions for the electrons and ions. In the presence of a magnetic 
field these additional terms represent the effect of the drift currents. More- 
over, the permeability €p becomes a slowly varying function of the co- 
ordinates and therefore, in the expression for the determinant D(k, œ), Eag 
must be differentiated with respect to r wherever it appears multiplied by k, 
which is equivalent to —iV. Let us assume that all these corrections have been 
made; we retain the previous designation D(k, w) for the resulting determinant. 
The other new feature in the propagation of the wave packets in an 
inhomogeneous plasma is related to a correlation of waves with wave vectors 
which are close to one another induced by the inhomogeneity. In a homo- 
geneous plasma, the components of the tensor £g in co-ordinate representa- 
tion depend only on the distance between the point of observation r, where 
the current is calculated, and the source point r’, the field at which excites the 
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current. Therefore in the k representation the components é,, are simply 
numbers. In a weakly inhomogeneous plasma, however, &é,, must be a func- 


: _/rt+r’ 
tion of both r and r’, but can be represented in the form &,, (r-r tr) 


where the dependence on the second argument is weak. Transforming to a 
Fourier representation, the dependence on the first argument changes into a 
dependence on k, while in the case of the second argument we may expand 
about the point r, where the wave packet is located. To first order we have in 
the co-ordinate representation 


r -tr 


r+r’ 
Exp (r-r 5) = &,,(r—r',r) + 5 ay Cael — r's E) 


where the differentiation of ¢ in the second component is performed only 
with respect to the second argument. Now transforming to the Fourier 


representation r’—r goes over into — and the second component will 


_ 6 
tak? 
have the following form 
B 1 ĝep(k, @, T) ; KA 
2 or ok 


Thus in the relationship between ¢,, and r, ¢ in eqn. (1.32) we must add 
some small terms proportional to the derivatives of the electric field with 
respect to k and œw. If we now repeat the derivation of the energy balance 
eqn. (1.36), we find that these additional terms lead to a term containing the 
derivative with respect to k of the spectral function. As a result eqn. (1.42) 
for a weakly inhomogeneous plasma assumes the following form (see refs. 
(12), (36a)) 

Ol, 0a, 01, da, oh, 

ot E Ok or or ok T Cera 
where ©, = @,(r, t) is the characteristic frequency obtained from the dis- 
persion relation D(k, œ, r, t) = 0. Its derivatives with respect to k and r 
can be determined from the rule for the differentiation of implicit functions 


oy, 2D [2D Gey _ _ aD [aD cs 

ôk  ôkj/ôðw ê ô (ae 
Note that for longitudinal oscillations the determinant D becomes simply the 
dielectric permeability e. 


Equation (1.43) has a simple physical meaning. The second term on the 
left describes the motion of the wave packet in space with a group velocity 


OW, 

Ux E ôk , 
distorted in such a way that the oscillation frequency remains constant. 

Note that we have assumed throughout that the wave vector and the 


frequency @, vary continuously. In fact the characteristic oscillation fre- 
quencies of an inhomogeneous plasma must be “‘ quantised ”, i.e. must pass 


and the third shows that simultaneously the wave packet is 
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through a discrete series of values (see ref. (36)), but since the difference 
between adjacent characteristic frequencies of short wave perturbations is 
very small, this effect can be neglected. 

It must be remembered that eqn. (1.43) corresponds to the zero order 
W.K.B. approximation. To this degree of accuracy it is completely irrelevant 
whether the term representing the divergence of the energy flux is written 
in the form shown in (1.43) or in the form div (Ux A). When we consider a 
large number of modes covering a broad band of wave numbers this approxi- 
mation is obviously quite sufficient. 

From the point of view of the quasi-linear approximation, we are in- 
terested in the case in which yx is small. In this case the time derivative of 
I, can also be considered small, so that in zero approximation we shall have 

Owy, Cl, Ca, Ol, 
ok or Or ok j ea) 

It follows, therefore, that i, = J (a,) i.e. the dependence of J, on r and k 
is mainly determined by the relationship of J(w) and @,(r). In other words, 
I, represents a set of waves with different frequencies, the amplitude of each 
of these waves is independent of r, and the corresponding wave vector 
k = k (r) is defined by the relation œx = œw, = const. In the next approxi- 
mation we must consider the dependence of J, on the time and a possible 
weak dependence on r. For waves with frequency w, we shall have 
h = I, (r, t) ô (@,—@,). Substituting this expression in (1.43) gives 


— 4 U.— = 2y,], (I. 46) 


Consider for simplicity a one-dimensional case when the plasma is 
inhomogeneous along one co-ordinate only, say x. Suppose moreover that 
@, is a symmetrical function of k,, i.e. to each frequency correspond two 
waves propagating in opposite directions. Suppose further that the waves 
are reflected without absorption either from the walls or from the turning 
points (points at which k, = 0). 

For small y and at the function J, can, according to (1.46), be considered 
independent of x. Multiplying eqn. (1.46) by U, ' and integrating with respect 
to x, we eliminate the second term and obtain 


I 
C = Ky, (L47) 


d 
where <y,> = MOA ~— is the mean growth rate of the wave a,. 
J Ue) | J U 


Clearly the stationary state will be reached when < y, > = 0, but for 
this to occur it is not necessary for y, (x) to vanish everywhere. It is now 
trivial to include in this scheme absorption at the walls, and if this absorption 
is large enough we may have a steady state or even damping, even when the 
local growth rate is positive everywhere in the plasma. 
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(c) Absolute Instability 


In the preceding section, we did not take into account the dispersive 
spreading of the wave packets. In several problems such an approximation 
is fully justified, but this is not always true. We shall consider here a specific 
example where the spreading of the wave packet must be taken into 
account. 

In Section 3(a) we considered a beam instability in an unbounded plasma. 
However, in practice all experiments with beams except experiments with 
toroidal systems are conducted in bounded tubes. The problem therefore 
arises as to whether an instability can develop in such a tube. If we approach 
this problem from the point of view of geometrical optics, within the frame- 
work given in the preceding section on non-spreading wave packets, we ought 
to describe the build-up of perturbations with the aid of eqn. (1.43). Accord- 
ing to this equation, in the absence of a feed-back between the entrance 
and exit of the beam, which can be realised, for example, by waves propa- 
gating against the beam (“‘ backward ” waves), any perturbation will move 
along the tube with a certain group velocity. The plasma and beam then 
Operate as an amplifier rather than as an oscillator. Such an instability is 
called a convective instability, in contrast to an absolute instability where the 
perturbation grows at every point in space (refs. 25, 37-41). 

We shall now show that the necessary feed-back can be produced by the 
dispersive spreading of the wave packet. We shall find that there is a bounding 
group velocity U, such that for U < U, the instability is absolute, while for 
U > U, it is convective. We shall follow the treatment of ref. (41). 

Suppose the growth rate y, as a function of k has a maximum at k = ko. 
Obviously, after the lapse of a sufficiently long time, any given initial per- 
turbation will be deformed so that it will have a sharp maximum at k = kọ. 
It is therefore sufficient to consider a wave packet consisting of plane waves 
with wave numbers close to kọ. It is then possible to make a series expansion 
of y, and œw, with respect to the small difference k-— ko: 


Oo 
Ve = Yko gz 5 (k— ko)’, CO, = Oka t U(k—ko) + É (k— ko)? (I. 48) 


Ory, 
a aae M . We now transform to 
Ok” lisi OK |k=ko Ok? |p ako i 


moving co-ordinate system in which œ,ọ = 0. The evolution with time of 
the wave packet is then given by 


where « = 


E(x, t) = Aro EXP (Yeo! — iko xX) 


| exp lic- Ut)(k—k,) — (T 


where Æo is the initial amplitude of the electric field. Provided that the 
difference x— Ut does not increase too rapidly with ¢, this integral can be 
evaluated and we obtain for the dependence of the electric field E on time the 
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following expression 
, =U a 
E(x, t) = Ax, exp- iko X+ Yrkt = Xa + B21 (a+ if) 
where the last component under the exponent sign describes the spreading 
out of the wave packet. As may be seen, for 


a? +p? 


U? < 2y,, =U? (I. 50) 
the amplitude of the wave packet will increase with time at each point 
x = const. In the opposite case U > U, the amplitude increases with time 
only in the moving system of co-ordinates. 

Since (1.49) is accurate only near the maximum, i.e. not too far from 
the point x = Ut, the expression (1.50) for the critical value of the group 
velocity is only approximate. 


(d) Resonance and Adiabatic Interaction between Waves and Particles 


Let us now consider in greater detail the question of the interaction 
between waves and particles within the framework of the quasi-linear 
approximation. As we have already stated in Section 3(a), it is convenient to 
distinguish a resonance and an adiabatic interaction of particles with waves. 
Each separate particle has a resonant interaction with those waves whose 
phase velocity coincides with the particle velocity. During resonant interac- 
tion the electric field of the wave is constant in a system of co-ordinates 
related to the particle, and therefore during such interaction a considerable 
energy exchange takes place between the particle and the wave. 

In the electric field of non-resonance waves, a particle performs oscil- 
lations with an amplitude governed by the electric field of the wave. A slow 
variation of the field amplitude leads to an adiabatic variation of the ampli- 
tude of the particle oscillations. 

In the example 3(a) which we considered, the adiabatic interaction was 
described by a diffusion coefficient (1.29). In order better to visualise the 
effect of the adiabatic interaction and at the same time to clarify its representa- 
tion in the form (I.29), we consider a simple example where the external 
field Ee~ ‘®' with k = 0 is imposed on the plasma. For a slow variation of the 
amplitude E, according to (1.29) we obtain 
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i 
Suppose that the square of amplitude of the field varies linearly E* = E3 F 


where T > œw™t. Equation (1.51) becomes a diffusion equation according 
to which the distribution function of the electrons, assumed cold at ¢ = 0, 
will be at time ¢ = T given by 


(I. 52) 
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An exact consideration would obviously give 


@ 


1 
flv) = <ô — ot) = m8 = TE (I. 53) 
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The difference between (1.52) and (1.53) demonstrates the inaccuracy of the 
quasi-linear approximation when it is used for a single wave. When the 
number of modes of oscillation is increased the accuracy of (1.52) will im- 
prove, since the distribution function f(v) as the probability density of the 
sum of a large number of random quantities must tend toward the Gaussian 
form. 

4ne*n 


For Langmuir oscillations, œ? = œĝ = , and the kinetic energy of 


nve 


m ! 
the electrons is equal to the potential energy of the electric field E7/8z. 


Thus, in the presence of an isotropic distribution of Langmuir oscillations 
the effective temperature T,,, which defines the width of the average distribu- 
tion function, is given by 7,, = T,+T,, where T, is the true electron tem- 

dT a? 
perature, and J, 3 Ban 
the oscillations and arises from the adiabatic interaction of the particles with 
the waves. Note that because T, > 0, the effective temperature 7, cannot be 
smaller than Tș. 

The division of the interaction between particles and waves into reso- 
nance and adiabatic interaction becomes slightly more complex when going 
over to an inhomogeneous plasma. Let us for instance consider the simplest 
case of the Langmuir oscillations of an inhomogeneous plasma in the absence 
of a magnetic field.t We assume that the effect of the mean electric field E 
on the oscillations can be neglected. In the quasi-classical W.K.B. approxi- 
mation the kinetic equation for the perturbed distribution function f,, is 
then written in the following form 


. The value of T, determines the kinetic energy of 
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where f is the averaged distribution function. In this equation the function 
of 


Ap and its derivatives with respect to r and ¢ are considered independent of 


r, t, and the distortion of the wave packet due to the inhomogeneity is in- 
cluded through the second and third terms in the curly brackets. 


t In an inhomogeneous plasma the longitudinal waves are decoupled from the other 
modes only when the electric field of the wave acts parallel to the gradient of the un- 
perturbed density. For any other direction the longitudinal oscillations are coupled with 
the transverse ones and the Langmuir oscillations must therefore be accompanied by 
electromagnetic radiation. However, the intensity of this radiation is small and the oscil- 
lations can still be considered approximately longitudinal. 
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The kinetic equation for the averaged distribution function has the 
following form 


af e af 
apt Nd na Oe (I. 55) 


where the term S,,, which represents the collisions between electrons and 
waves, is given by the relation 


ð 
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The expression underneath the integral sign in (1.56) can be written in the 


symmetrical form 4# EÉ ofro + Exo feo and we can then substitute for the 
function fko its expression in terms of Eko from (1.54). The contribution 
due to the first term in the right hand side of (1.54) is proportional to 
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Bearing in mind that the correlation of the electric field approximates to a 
delta function, we can replace the real part in the round brackets in (1.57) by 
2x0 (w—kv), and the imaginary part can be represented in the form 
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The last two terms in the curly brackets in (1.54) lead to terms which are 


k 
proportional to the derivatives of Jke with respect to k and œ. Since Ex = Ey 


we finally obtain 
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where the singularity (o—kv)~* is integrated in the sense of its principal 
value. We have omitted a small term proportional to Iko (i.e. not containing 
a derivative of J.) because within the framework of our adopted zero order 
W.K.B. approximation we also neglected similar additions in the transfer 
eqn. (1.43). In (1.59) the first integral corresponds to the resonance interaction, 
the other two to the adiabatic interaction between the particles and the waves. 
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Multiplying (1.59) by mv and integrating the result with respect to v, we 
obtain the force F acting on the plasma due to the interaction with the 
waves. Since for longitudinal oscillations € = 0 we obtain 
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where e is the dielectric permeability of the plasma given by 
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According to (1.43) the second integral in (1.60) vanishes and consequently 
the force F reduces to the divergence of the Maxwellian stresses. The second 
part of the force, related to the kinetic energy of the oscillations, is obtained 
from the averaged function and, according to (J.52), can be represented in the 
form of a gradient of the pressure tensor 


(I. 61) 
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where we have replaced w* by wo = 
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The total force exerted on the plasma by the Langmuir oscillations is 
thus simply the negative gradient of the energy density of the electric field. 

Consideration of the adiabatic interaction of the particles with the 
waves 1S important in other cases when the energy and momentum balance 
in a plasma supporting oscillations is discussed; in particular this interaction 
must be included when discussing the “ anomalous ” diffusion of a plasma. 

It must be remembered that the division of the total interaction into 
resonant and an adiabatic parts is only possible for sufficiently small growth 
rates when one can refer to near-periodic oscillations. The condition 
y < max (w, kv;) must be fulfilled and this is also a necessary condition for 
the applicability of the quasi-linear approximation (v+ is the thermal velocity). 


(e) Enhanced Diffusion of a Plasma 


The expression for the collision term (1.59) can be generalised without 
difficulty to include longitudinal oscillations of a plasma in a magnetic field. 
For a two-component plasma an expression of the type of (1.59) must be 
written down for each component. The force F, acting on the component j 

Fo 
can then be represented in the form F; = ae F; where Fo is the total 
force acting on the, plasma as a whole and F; is the “ frictional ”? force 
between the electrons and ions which is transmitted through the oscillations. 
By using a relation of the form of (1.60), it is not difficult to show that the 
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total force Fp is the divergence of the Maxwell stress tensor. This force 
can only lead to a macroscopic motion of the total plasma, including the 
frozen-in magnetic field. The diffusion across the magnetic field is determined 
by the frictional force between the electrons and ions, F;,, or more precisely 
its component across the density gradient. This conclusion follows immedi- 
ately from the hydrodynamic equations of motion for each of the components 
which, as we know, describe fairly accurately the slow mean motion of the 
particles across a magnetic field. 

This result 1s very important and must always be borne in mind when 
investigating the possibility that some given oscillations affect the diffusion 
of a plasma. In particular, it follows that high frequency oscillations, in 
which the ion motion can be neglected, cannot lead to diffusion of the plasma. 
At first sight this assertion might seem improbable, since each separate 
electron in such a field performs random motions which may be considered 
as diffusion. The corresponding coefficient of diffusion for a separate 
particle can be determined either by calculating the mean-square of the 
displacement (Taylor (42)) or from quasi-linear theory (43); for resonance 
particles this diffusion is described by the first term in (1.59). But if together 
with the diffusion of the resonance particles we also consider the displace- 
ment of the remaining electrons due to the adiabatic interaction, we obtain an 
expression of the type of (1.60), according to which the total current vanishes 
identically. Thus, the interaction of electrons with waves, if the ions are 
stationary, leads only to a diffusion of separate particles, i.e. to effects such as 
an intensified thermal conduction, rather than to a diffusion of the plasma 
as a Whole. The interactions of these oscillations are therefore equivalent to 
electron-electron collisions. 

The neglect of this important conclusion in papers dealing with enhanced 
diffusion has often led to erroneous statements. For instance, in some experi- 
mental papers attempts are made to relate the anomalous diffusion directly to 
the high frequency oscillations in which it is known that ions cannot partici- 
pate. Theoretical considerations of enhanced diffusion are sometimes 
limited to only one component, for example, the electrons; in this case the 
problem of the mechanism making the diffusion ambipolar needs further 
consideration. This problem is automatically resolved when we determine 
the diffusive loss directly from the frictional force between the electrons 
and the ions 

A similar error is incurred by not taking into account the “ dragging 
along ’’ of the waves by the particles, as a result of their interaction. For 
instance, in Taylor’s paper (42) quoted above where Langevin’s equation 
is used for the investigation of the random motion of the particles, a rather 
vague assumption was made that the oscillations are isotropic in the labora- 
tory system of co-ordinates. This assumption, which is equivalent to assum- 
ing a strong coupling between the oscillations and the walls surrounding the 
plasma, led that author to the erroneous conclusion that the coefficient of 
enhanced diffusion cannot exceed Bohm’s value. A similar error was also 
committed in ref. (44) where the coefficient of diffusion of the electrons was 
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calculated on the basis of the thermal cyclotron oscillations, which were 
considered isotropic in the laboratory system of co-ordinates, rather than 
in the system of co-ordinates moving with the electrons (this error has been 
corrected in a later paper by the same authors (45)). 

In a complete discussion of the diffusion problem all these paradoxes are 
automatically resolved, and complete clarity, as we have seen, can only be 
achieved when the adiabatic interaction between the particles and the waves 
is considered. In those cases in which resonance interaction is impossible, 
the enhanced diffusion can be determined from the adiabatic interaction 
alone. As an example, we may quote the problem of the diffusion of a plasma 
due to drift waves excited by an external source, treated by Petviashvili (46), 
in which the electrons diffuse due to the resonant interaction with the waves, 
and the ions due to the adiabatic interaction. 


Il. INTERACTION BETWEEN WAVES IN WEAK TURBULENCE 


THE quasi-linear method investigated in the preceding chapter has a very 
limited field of application, being suitable only for the description of states 
which are so weakly excited that strictly speaking they ought not to be called 
turbulent, since the most important property of turbulence—the non-linear 
interaction between the oscillations—is not yet apparent. In real turbulent 
processes in a plasma, the interaction between the oscillations generally 
plays an important part. The consideration of this interaction is the subject 
of this and the following chapters. 


1. KINETIC WAVE EQUATION 


(a) Derivation of the Kinetic Wave Equation 


We are considering here only states where the interaction between the 
waves can be considered weak. This situation can be expected when the 
growth rate of the unstable perturbations is sufficiently small, i.e. y/w < 1. 
This condition as it stands, however, may be misleading because the value of 
w depends on the choice of the co-ordinates system. For instance, in the 
case of a set of sound waves travelling in the same direction, say parallel to 
the z axis, the frequencies of all the waves vanish when we transform to a 
co-ordinate system moving with the sound velocity. Even for a very small 
growth rate this must be considered a case of strong turbulence. However, 
for the same sound waves distributed isotropically, the weak interaction 
approximation can be used. Thus, although in what follows we shall use the 
condition y/@ < 1 as the weak coupling criterion, we must remember that the 
condition would be expressed more accurately by the statement that the 
growth time of the perturbations must be considerably larger than the 
characteristic time of conservation of the relative phases of the different 
waves. 

Let us first consider oscillations in which the resonance interaction with 
particles plays no part. This would occur, for example, in the case of a cold 
plasma where the thermal motion of the particles can be neglected, and the 
magneto-hydrodynamic equations can be used to describe the motion of the 
plasma. Now a whole series of such problems lead to equations whose basic 
structure is similar, and can be represented in Fourier space by the following 
general scalar equation 


(œ — Wy, — in) Ce = | Vko, k'o’ Cko Ck -—k’,@-o’ dk’ da’ (dl è 1) 


We shall use (II.1) as a model equation to explore the character of the 
31 
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interaction between the oscillations. The expression on the left hand side of 
eqn. (II.1) represents the linear part of the equation of motion and the 
expression on the right represents the Fourier transform of the quadratic 
terms. These quadratic terms describe the interaction between the various 
harmonics and the quantity of Vko, wa represents the matrix element of this 
interaction. 

Note that since Cko represents the Fourier transform of a real function, 
Cx, = C-k,- Using this relation in the complex conjugate of eqn. (II.1), 
changing the signs of k, œ and k’, w’ and noting that @-, = —@,, Y-k = Yks 
we obtain an equation analogous to (II.1), but with a different matrix element. 
Comparing these equations we obtain 


Viet = Vise —-w; —k’, -o’ (II. 2) 


If the growth rate yą > 0, small perturbations will increase with time 
until the non-linear interaction comes into play. In conditions approximating 
to equilibrium, the right hand side of (II.1) must be considered small, since 
in zero approximation we have 


(o—w,)C? = (II. 3) 


If the oscillations developed from random thermal motions, the individual 
modes can be considered completely independent. Multiplying equation 
(II.3) by Ci?2* and averaging the result with respect to the statistical ensemble, 
i.e. With respect to the random phases of the separate oscillations, we obtain 
for a stationary and spatially homogeneous system of oscillations 


(Cko Cro = Delo — olk- klw- w’) (II. 4) 


If we replace Cko by C? on the right hand side of (II.1), the non-linear 
term will play the part of an inducing force. The amplitude of the induced 
oscillations will be denoted by C) The presence of induced oscillations 
does not lead directly to damping of the waves, which becomes evident only 
in higher approximations. It follows, therefore, that the quantity CV must 


be of order a y/w, and therefore in the equation for CÉ} the growth rate y, can 
be omitted. Moreover, since the beats must attenuate more rapidly than the 
main oscillations, it is necessary to add to the right hand side of equation 
(II.1) for the beats in place of — iyẹ a small term iv, which represents the damp- 
ing due to the higher correlations. We thus obtain the following equation for 
the first approximation 


Ch) = (@-O tiv) | Wowo C2 Cro- dk’ do” (II 5) 


By a similar method we could determine the higher order corrections to 
the amplitudes. However, we are not so much interested in the correction 
to the amplitude as in the effect of the wave interaction on the characteristic 
frequencies of the oscillations, and more precisely on their imaginary part, 
i.e. the damping. To determine the magnitude of the additional damping 
we multiply eqn. (II.1) by Cx. and then average the result with respect to the 
Statistical ensemble, assuming, as before, that in zero approximation the 
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oscillations are uncorrelated (the random phase approximation). The right 
hand side, containing the product of three random variables, vanishes in zero 
approximation; the next approximation is to include successively in each of 
the three factors the first order correction C, which leads to 
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where k” = k—k’, w” = w— w’. 

The first two terms on the right hand side of this equation are proportional 
to Iko and contribute therefore to altering the characteristic frequencies. 
For our purposes it is sufficient to consider only the imaginary part of these 
terms, since it is this part which describes the damping of the waves due to 
the non-linear interaction. Their real parts, giving the frequency, can be 
neglected on the strength of the condition y/w < 1. Transferring the terms 
proportional to Jo to the left, they can be combined together with 7, 4,4 and 
we obtain the total growth rate ), (more precisely speaking the damping 
rate, since in the steady state the value of ), must be negative). Since the 
non-linear terms are small we can, within the framework of our approxima- 
tion, replace the small imaginary part v > 0 by —}), > 0 in the last term, 
and the relation (II.6) can then be rewritten in the following form: 


{(w 7” wy)“ a Fk ko = 4 Í liis, vol lra Ikro" dk’ dw’ (II. 7) 
where Uxe, ka” = Vio, wo + Vko, ko” and the right hand side of (II.7) is now in a 


symmetrical form. 
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as previously Iko = Ikô(@— wy). Substituting this last relation in (1.7) and 
considering that Im (w — wy +i»)! = —nô(w — wy), we obtain the kinetic 
equation for sustained oscillations : 
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where the integral with respect to k’ in the last term on the left hand side is 
the principal value, and k” = k—k’, œw” = w— ow’. 7 
In a transparent medium vx,” is real. For simplicity we shall only 
2° 
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consider this case here, leaving the general case to Section (2) where we shall 
consider the thermal motion of the particles. 

Note here that the addition to the right hand side of eqn. (II.1) of a cubic 
term of the form 


f Vko, k’o’, k”œo” Cko’ Ckro” Ci -—k’-k’,@a-w’-o” dk’ dk dw dœ” (II. 9) 


would lead to the appearance on the right hand side of eqn. (II.6) of an 
additional term 


Í Uko Ko Teco ty’ AK’ doo" (II. 10) 


where 
Uko, k’o’ — Vio; k’o’, ~k- w’ T Vko, kw, k’o’ + Vko, k’w’, ko (II. Í 1) 


Equation (II.8) defines the spectrum of the oscillations in the steady state, 
where the growth of the waves due to the instability is exactly compensated 
by their damping due to the non-linear interaction. In the absence of equili- 
brium the oscillation amplitude will vary with time with a growth rate equal 
to the difference between the linear growth rate and the non-linear damping. 
The kinetic equation for the waves, including non-linear effects, can therefore 
be obtained from (1.43) simply by replacing Yx by p, to give 
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This equation, together with the equation for the averaged function, 

constitutes the basis for describing weakly turbulent systems. In the case of 
; 00, : ‘ : 

a homogeneous medium, where a 0, this equation is analogous to the 

well known kinetic equation for phonons, the sound quanta in a solid (see 

ref. 11). The kinetic equation for waves has only very recently been applied 

to a plasma (refs. 10 and 13). 

In this wave eqn. (II.12), only “‘ three-wave ’’ processes are considered; 
these are the decay of the wave k into k’, k” and the inverse process of the 
merging of the two waves k’, k” into one. The 0-function of the frequency 
difference appearing in (II.12) requires @, = @,+@ ,, and this condition, 
which can be shown to represent the conservation of energy, considerably 
limits the permissible region of interaction in k-space, and indeed for many 
forms for the dispersion relation œ = w(k) three-wave processes are com- 
pletely forbidden. It is natural therefore to divide the possible dispersion 
relations into two groups, decay and non-decay relations, according to 
whether decay of one wave into two is or is not permitted. 

In an isotropic medium, for instance, a relation of the form 1 (Fig. 10) 
where the phase velocity decreases with k, is a non-decay relation, while for 
a spectrum of type 2, where the phase velocity increases with k, the decay 
conditions can be satisfied. 
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For a non-decay dispersion relation, we must go to the next approxima- 
tion to obtain a non-vanishing interaction, and we then consider “ four- 
wave ” processes in which the scattering of two waves on one another gives 
two further waves. In the kinetic wave equation such processes lead to 
terms which are cubic in the spectral function J. But as we shall show below, 


Dy 
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Fic. 10. Dispersion relations of decay (2) and non-decay (1) types 


when the thermal motion of the particles is taken into account, additional 
terms quadratic in J appear, and these will in general dominate the cubic 
terms. We shall not, therefore, consider non-decay dispersion relations in any 
greater detail here. 

Let us now discuss the two simplest examples of decay interactions in a 
plasma. 


(b) Interaction between Langmuir Waves and Ion-sound Waves 


Let us consider the simplest possible case, the longitudinal oscillations 
of a plasma in the absence of a magnetic field. We shall assume that the ion 
temperature T; is considerably smaller than the electron temperature T,, so 
that we can put approximately 7; = 0. Under these conditions both Lang- 
muir and ion acoustic waves can propagate in the plasma. Each of these 
oscillations taken separately has a dispersion relation of the non-decay type. 
The dispersion relation for ion acoustic waves is of the form of curve 1 in 
Fig. 10, and the frequency of Langmuir oscillations approximates to @ 9, so 
that the sum of three frequencies cannot vanish. Therefore the only possible 
three-wave processes are those in which Langmuir waves are scattered by 
ion-acoustic waves, and we shall consider these in detail. 

In the hydrodynamic approximation the oscillations are described by 
continuity equations 


On; 
Ot 
for the electron density n, and the ion density n;, the equations of motion for 


+ div(n,v,) = 0 (II. 13) 
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each of the components 


dv; e; 
mnj + VaT) =— nVo (iI. 14) 
J 
and the equation for the electric potential 
Ag = —4ne(n;—n,). (II. 15) 


We shall use the index s for the ion-sound oscillations, the index / for the 
Langmuir oscillations. Let us consider first the Langmuir oscillations. In 
the case of the Langmuir oscillations the ions can be considered stationary, 
i.e. n} = 0, so that eqn. (II.15) takes the form 


Ag! = 4ren! (II. 16) 


Further, since decays within the Langmuir oscillation branch are pro- 
hibited, the non-linear term in the electron continuity equation and the 
non-linear term quadratic in velocity in the equation of motion for the 
electrons can be omitted. The non-linear coupling between the Langmuir 
and acoustic oscillations is given by a term on the right hand side of the 
electron equation of motion (II.14). Transforming to the Fourier representa- 
tion we express n, in terms of g from (II.16). Taking the divergence of the 
electron equation of motion and expressing div v, in terms of n, (II.13), we 
obtain the equation for the potential of the Langmuir oscillations 
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where œ; is the square of the frequency of the Langmuir oscillations given by 


T, ; , e 
w? = w + ~ kt, For k, k’ small compared with the inverse Debye length 


the second term in the curly brackets in (11.17) is negligible, and recalling that 
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where k” = k-k’, œ = w— w’. 

For the acoustic oscillations, the ion continuity equation and the equation 
of motion can be linearised, since decays within the ion branch are again 
prohibited. Expressing n; in terms of g* from the linearised equations, we 
obtain from (II.15) 


Q2 
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where 04 = an oe the square of the ion Langmuir frequency. At small k 
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the frequency of the ion-sound oscillations is considerably less than Qo, 
so that in (II.19) the unit in the brackets can be neglected. 

To obtain the coupling between the ion-acoustic and Langmuir oscilla- 
tions, we must take into account non-linear electron terms which would 
give low frequency beats. Taking the divergence of the electron equation of 
motion and noting that the term in V@g can now be linearised, we obtain in 
Fourier representation 


nm J | bie) (Kvi ro” — fees +) dk de +k?T ni, = kneg? (II. 20) 


Since the frequency of the ion oscillations w, is small, we have approxi- 
nol 
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in terms of @* we obtain 
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where œw; = ci k* = s . This equation describes the interaction between 
l 
the Langmuir and ion-acoustic oscillations. 

In the linear approximation, to each wave vector k correspond two fast 
waves, @x = +a@,, and two slow waves, wy = tœ, The ee function 
of the potential for these oscillations will be denoted by /,7' and I‘ respec- 
tively so that, for instance, 


s o = IE O(w—kc,) +1, slw + ke,) (II. 22) 
In addition, since I, = I-k, -œ, we shall assume that w-, = —@, and 
consequently the transformation œw, k — — w, —k changes neither I, nor 


Ô (w — x). 

In principle it is possible to choose another method of differentiating 
between the two waves with the same wave vector, but propagating in oppo- 
site directions. Often, for instance, the frequency is defined to be positive 
and then the direction k defines the direction of propagation of the wave. 
In this system, however, a more complicated notation for the collision term 
between the waves is required. The transition from one representation to 
the other is simple and we shall use only the representation introduced above. 

Near each of the characteristic frequencies, the difference w?—@, can 
be approximated by 2a, (w — wx), and the eqns. (II.18) and (II.21) assume the 
form of the model eqn. (II.1). The kinetic equations can therefore be written 
down by analogy; we obtain the following system for the longitudinal 
oscillations of a homogeneous plasma 


oli. ôl new 
ôt kõr 4 T? 


è f cos’ al R. — IL Aelio — og) dk’ (II. 23) 
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It is evident from these equations that the main part in the interaction 
between the waves is played by the electron oscillations: beats between these 
oscillations excite ion-acoustic waves which then scatter the electron 
oscillations strongly. Let us for instance consider the case where a single wave 
with wave vector kọ propagates in the initial state in the plasma, so that 

t = Po Ôl — wk, )ðlk— ko) (Oraevskii and Sagdeev (47)). Such a wave may 
decay into a Langmuir wave with k’ % — kọ, propagating in the opposite 
direction and having a frequency œw} close to œw}, and an ion-sound wave with 
wave vector k” % 2k,. In this case the condition w, — œ} — @y = 0 can be 
satisfied. 

To investigate this process it is convenient to revert to the eqns. (11.18), 
(11.21) for the amplitudes. For each of the excited waves we write w = @, tiy 
and linearise the equations for the amplitude relative to the perturbations 
o'd(k +k )5(@—w}) and 0°6(k—2k,)6(w — oi), we obtain 


e 
2iyw' p! = —we F p? (II. 25) 
2 
, 5 COS 
ipo, p = — kg Qo (II. 26) 
m,W9 


whence, assuming œw’ % wọ we obtain 
2 2 
> Ik o, E 


T, 


e 


(II. 27) 


Thus y can be > 0 and for a single Langmuir wave a “ decay instability ” 
occurs, in which the amplitude of a small perturbation representing a super- 
position of Langmuir and lon-acoustic waves, related to the main wave by the 
decay conditions, increases exponentially with time. 

We could describe this process with the aid of the kinetic equations for 
the intensities J’ and J°. However, in this case a difficulty arises because 
the substitution in the kinetic equation of functions of the form /,6(k—ko), 
I'6(k +k) and I°6(k + 2k), leads to a divergence of the quadratic term. But 
if we recall that the function 6(@, — ww — wy) Was in origin an approximation 


to a function of the form [(w — wy)? + y7]7', we see that this function must 


be replaced at the resonance point by 1/zy and we again obtain (II.27). 
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This example shows again that the kinetic equation for the waves is 
strictly correct only for wave packets which are sufficiently broad in wave 
number space, with a smoothly varying function J. 

According to the kinetic eqns. (11.23), (11.24), the Langmuir oscillations 
alone cannot interact with one another and they must first excite ion oscilla- 
tions. This process is a relatively slow one. But as soon as the amplitude 75 
increases appreciably, the main process is the rapid scattering of the Lang- 
muir waves at the ion inhomogeneities. Since the interaction between the 
Langmuir waves depends explicitly on the electron temperature, we should 
expect that an important part in this process may be played by processes 
involving the thermal motion of the electrons, which will be considered in 
Section 2. 


(c) Interaction between the Alfvén Waves and Magneto-acoustic Waves 


As a second example, we consider the oscillations of an ideally-conducting 
plasma in a homogeneous magnetic field (49). For simplicity we assume that 
the plasma pressure is much less than the magnetic pressure, so that in the 
equation of motion the plasma pressure can be neglected. Then the plasma 
oscillations will be described by the following system of magneto-hydro- 
dynamic equations 


dv 1 
myn = 7 [curl H, H] (II. 28) 
as curl [vH] (II. 29) 
ôt 
ôn 
Moi II. 30 
F iV ny ( ) 


Suppose that in the stationary state the plasma is quiescent (v = 0) 
and its density is constant. Linearising and transforming to a Fourier 
representation the first two equations become 


1 
ZM NOko = 7> [[kH;o]H] (II. 31) 


— OH, = (kL v2 HJ] (II. 32) 


The component of the velocity parallel to the mean magnetic field H has 
disappeared from these equations, so that it is sufficient to consider only 
oscillations transverse to the field. Expressing Hy. in terms of Vko with 
the aid of (11.32) and substituting the result in (II.31), we obtain the equation 
for the velocity 

Vip = CAL KZ Vo + K(KV 0) | (II. 33) 


2 


where c4 = is the square of the Alfvén velocity and k, the component 


t 


of the wave vector parallel to the unperturbed magnetic field H. From this 


40 PLASMA TURBULENCE 


we obtain for the velocity component Vko in the plane (kH), 


w? = œw = c4 k? (II. 34) 
and for the component Vko, perpendicular to both k and H we have 
w? = œw? = cå k? (l1. 35) 


Thus, for a low pressure plasma, i.e. 6 = 8np/H? <& 1, two types of 
waves can propagate, namely Alfvén waves with a frequency w,, and magneto- 
acoustic waves with a frequency œ, In addition stationary perturbations 
of the density (w = 0) are possible; these correspond to the slow magneto- 
acoustic waves which are obtained when the plasma pressure is not neglected, 
but these will not be considered here. 

The phase velocity of the magnetosonic waves is a constant c4. It 
follows, therefore, that the non-linear interaction is very important for waves 
propagating in the same direction. This interaction leads in particular to the 
steepening of the wave fronts of an initially sinusoidal wave by the generation 
of higher harmonics. But for waves propagating in different directions, the 
dispersion relation does not permit decay processes. 

The Alfvén oscillations propagate only parallel to the magnetic field H, and 
this with the same phase velocity c4. We might expect, therefore, that these 
oscillations would interact strongly with one another. However, the matrix 
element of the interaction of two Alfvén waves propagating in the same 
direction can be shown to vanish, and as a result a single Alfvén wave can 
propagate with any finite amplitude. 

Thus, neither the Alfvén nor the magneto-acoustic waves by themselves 
are decay waves. However, each of them may decay into a pair of Alfvén 
and magneto-acoustic waves. The matrix elements of this interaction can be 
obtained from the equations of motion (11.28)-(11.30). Because of their 
cumbersomeness we shall not quote the general form of the expressions 
(these can be found in ref. (49)), but shall confine ourselves to the limiting 
case of almost transverse propagation where k, < k. Neglecting in the matrix 
elements the quantity k, compared with k, it is easy to transform the equations 
for vê and vf into the following form 


kk’ ; 
(w —k*c4)vi,, = le T (ko — Vko Vern) SIn & COS & — 
— tko Verge Sin? a} dk’ dw’ (II. 36) 


kk’ 
(w* — k2c4)ve,, =| 7 (SIN & COS 0 Vio Ung” F Dire Ving» SIN? + 


+ Uf rey! Vro» COS? a) dk’ dan’ + | Fesinacosax 


, oki —a'k2\ , 
x a Uk ew’ Veg” AK’ da’ (II. 37) 


[kk]... 
kk’ 


where sina = 
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In the equation for v* we have deliberately conserved the second integral 
term although in fact it vanishes for oscillations of such low intensity that 
the frequency shifts are negligible, so that we can write 

w? = cf k, œ? ack 
In actual fact, at very small k, and kj, and a sufficiently large amplitude of 
the oscillations, we may find that the spectral function {v{,.|* is broadened to 
cover a frequency range which may be comparable to or even larger than 
k.c,. In this case we can no longer neglect the Alfvén—Alfvén interaction 
because a transition takes place to strong turbulence where the restriction 
on “ three-wave ” processes due to the 6-function of frequency disappears. 

We shall not here write out the kinetic equations for the Alfvén and 
magneto-acoustic waves, since they can be obtained without difficulty by 
analogy with the case of the Langmuir and ion acoustic waves considered 
above. It is sufficient to note that according to eqns. (11.36) and (11.37), the 
matrix element of the interaction is of the order of unity for vf, vu’ ~ c4. It 
follows, therefore, that the lifetime t for decay in the case of a single wave 
with amplitude vo is determined by the quantity c4/@,, svo, and the character- 
istic energy exchange time between different modes in a diffuse wave packet 


is of the order of c4/@,. Vo. 


2. INTERACTION OF WAVES IN A PLASMA WITH CONSIDERATION OF THE THERMAL 
MOTION OF THE PARTICLES 


In the derivation of eqn. (II.8) we did not take into account the specific 
properties of a plasma, i.e. a system of charged particles interacting with an 
electromagnetic field, as the medium supporting the oscillations. However, as 
we know from linear theory, the discreteness of the medium gives rise to a 
specific damping of the waves, namely Landau damping, related to the reso- 
nant interaction of the waves with the particles. We should naturally expect 
that some corresponding effect should appear in the non-linear case. 

As we have seen above, the principal effect of the non-linearity 1s to lead 
to the appearance of beat oscillations with combination frequencies w— w’ 
and wave vectors k—k’. The resonant interaction between the particles 
and these beats gives rise to an additional damping of the waves which we shall 
call non-linear Landau damping. In the terminology of waves and particles 
this process corresponds to the scattering of waves by particles, i.e. in 
absorption of the wave k and re-emission of the wave k’, whereas the linear 
Landau damping corresponds simply to an absorption of the wave. This 
non-linear damping effect is of the same order of magnitude as the decay 
processes, and when the dispersion relation is of the non-decay type it is the 
most important non-linear effect. To consider the non-linear interaction of 
the waves with the particles we must set up a kinetic equation for the waves 
on the basis of the non-linear kinetic equation for the particles. We then 
obtain simultaneously the effect of the waves on the averaged distribution 
function. 


42 PLASMA TURBULENCE 


(a) Kinetic Wave Equation with consideration of the Thermal Motion of the 
Particles 


In order to keep the discussion as simple as possible we consider, to 
start with, the case of longitudinal Langmuir oscillations, and we then show 
how the equations can be extended to the general case. Suppose that in the 
equilibrium state the plasma is homogeneous and there are no electric and 
magnetic fields. We divide the electron distribution function into two parts, 
one averaged with respect to time f and one oscillating part f” which we 
expand as a Fourier integral. Separating the kinetic equation into two by 
taking suitable averages, we obtain 


0 0 
D Wisi — | k(o% 4 feo» dk dw dk'da’ (II. 38) 
Ct . m Ov 
e . of 
—i(o—kv) fko = — z K Oio” 
eð f i 
a m ðv k (Pro Sk-k', o-o (Oxo Sk-k’, 0-0") dk dw (II. 39) 


where Øko 1s the Fourier component of the potential of the electric field, given 
by 


4ne 
Pko = TE | fo dv (I. 40) 


Neglecting the quadratic terms in the eqn. (11.39), we can express fx. 
linearly in terms of Øko: 


e ôf 
= — —k ] -lik = ° 4 
Sivoo en (w v+iv) ay Pko (II 1) 


where we have introduced the small positive quantity v > 0, allowing for a 
small damping, in order to pass round the pole correctly. Substituting this 
expression into (11.40), we obtain 


e(k, @)Pko = 0 (II. 42) 
which gives the dispersion relation e(k, w) = 0. The quantity 
ô 
2 k 
osip -n II. 43 
,%) = — | —— ôv ; 
mk? J (w—kv-+ iv) ( ) 


represents the dielectric permeability of the plasma. 

Now let us consider the solution of the non-linear equation assuming 
that the amplitude of the oscillations is small. We write eqn. (11.39) in the more 
compact form 

Feo = (Bko kK) FP xo + J (Pico k’) (Pro Sk =k, 0-0 — (Pra fk-r, o-o) Ak’ doo’ 
(II. 44) 
where g represents the following operator 


e ô 
g= (o=kv+ i)! = (11. 45) 
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In eqn. (II.44) the quadratic terms can be considered small. Making 
use of this, we integrate this equation and represent /,., in the following form 


fro = (Siok) Ox. + f (Gre k’) (8r-r(k-KD Pro PỌk-k',o-o' 
— [Pko Px-k',a-0$ dk! doo’ +| (Skok )(8k-r' K”) X 
X (8k-w-k(K-— kK — KYS Pwo Pko” Pk-k'-k”, wœ- - o" 
— Pro Pura” Pk-r' =k”, w- w’ -07 — (Po Pra” Pk -k'- k”, 0- o’ ae x 
x dk' dw’ dk” dw” +... (II. 46) 


. Within our approximation we can neglect the terms of the fourth and 
higher orders in the oscillation amplitude in (II.46). Substituting this equa- 
tion into (II.40) we obtain the following non-linear equation for the potential 
Pko- 


e(k, WO) Prey a | Vag: ko (Pk! Py, - k, o-o (Pko Px —-k’,@- >, dk’ doy’ T 
Ez Í Vko, k'o’, kor Pko Pko” Pk- k'o o 
= Pko Pure” Ox -k’- k”, œ- o -— or? T Pko’ Ox" wa” Pk-k'- k”,@o-@’- o>} x 
dk’ dw’ dk" dw" (IL. 47) 
where the matrix elements are given by the following relations 


4ne f 


Moone! =~ x | Grok’) Be-x,0-o(k—k Sav (UI. 48) 


4ne 
eee =| (Grok) (By -wk” (Gy ye -a(k—k’—k") dv (I. 49) 


Equation (11.47) is similar in structure to our model eqn. (II.1), supple- 
mented by a cubic term of the form (II.9). We can therefore immediately 
use the results of the preceding section and, by analogy with (II.6) and (II.7), 
obtain the following kinetic equation for the waves with consideration of 
the thermal motion of the particles 


Are 
a(k, uw = Iho po f (2,6) een E Bye -w ft 
tket N hadak doo! + 


Uko, k” o” UK a", ko 
+I of PakO Tr dK deo! + 
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fi Vko ro'l Two lo” dk'da’ (I. 50) 
w) 


T 36 A 
where 
k” = k-k’, w” = w— w’, e+(k, œ) = e(k, œw + iv) 
Ukoko = Vko, ro T Vko, kor = Uko, k'o” 
The quadratic terms on the right hand side of eqn. (II.50) describe the 
additional damping and the shift of the characteristic frequencies due to the 
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non-linear interaction between the waves. For y/œ <% l, i.e. e’/e <1 
where ¢’ = Ree and e” = Ime, eqn. (II.50) can be solved by the method of 
successive approximations. The zero approximation gives £€fkoe = 0 which 
has the solution o = 4, 6(@—@,), where w, is the characteristic frequency. 


In the next approximation it is sufficient to consider only the imaginary part 
de’ \_ 

of eqn. (11.50). Denoting by y, the linear growth rate yk = — €” TF ! and 

adding a term including the derivatives with respect to time, we write down 


the imaginary part of egn. (11.50) in the following form 


1 él Ge'\ T! 47e 
oh = hht h(E ) 1 2 f Bu) geo Wei 
2 ot OW 


+ (k’ e i Mw dv dk’ + 
de'\~* Uko, ko" Vk”o”, ko Lk’ 
Li — l ea dk! + 
u Com | elk”, @, — ©, tiv) 


E N) orcs 5(@, — Oy — Oe ye Tye dk’ (I. 51) 
2\00/ kon 

The last term and the subtracted part of the penultimate term describe 
the wave decay processes and the remaining non-linear terms describe the 
scattering of the waves by particles; they include the linear Landau damping. 

We have written down this equation for the electrons only. To include 
the ions it is sufficient to consider their contribution to e and to the matrix 
elements vko, k'o and also to sum the second term in the right hand side of 
(11.51) over the two particle types. No difficulties are encountered in genera- 
lising these equations to the cases in which the plasma is located in a magnetic 
field: the only change is in the form of the operators ko. The generalisation 
to the case of randomly polarized oscillations gives considerably more 
complicated expressions, because the full dispersion function given by eqn. 
(1.33) must be used in place of ¢ and the one eqn. (II.51) must be replaced by 
separate equations for each of the possible polarizations. 


(b) Thermal Fluctuations 


So far we have throughout understood by f the distribution function 
averaged over small macroscopic volumes in phase space and therefore a 
continuous function. Similarly we have omitted from consideration the 
thermal fluctuations related to the discreteness of matter. We shall now drop 
this assumption and discuss the effect of thermal fluctuations. 

The kinetic Vlasov equation without the collision term can be also written 
down for a microscopic distribution function 


Salt, V, t) = ), Or —r,)d(v—v,) 


where the summation is over all particles of the given type. Let us split up the 
function f,, into two parts: an averaged part f and a fluctuating part f”, which 
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vanishes when averaged over the macroscopic volume. Neglecting the inter- 
action between particles, which is admissible for a rarified plasma for which 
the Debye number Np = nD? > 1, we obtain the correlation coefficients 


<f”, v, DSE, v, t) = dv—v')d(r—r’ — v(t—t’)) f(v) (II. 52) 


in co-ordinate representation and 


Siro Skol) = Òk -klw — w')5(@ — ky) 


in the Fourier representation. 

The correlation function (11.52) corresponds to freely moving particles 
and satisfies the kinetic equation with no interaction. In the presence of 
oscillations the macroscopic function f(r, v, tf) also becomes random, i.e. 
macroscopic fluctuations arise, and only these have been considered earlier. 
It is easily seen that the thermal fluctuations were omitted when we expressed 
fko in terms of Ox. In fact the general solution of the linearised kinetic 
equation for the longitudinal oscillations has the form 


ari a SS eee IT. 54 
Suc m dvo—kv+iv + Seal) ( ) 


JO) 
(27)? 


(II. 53) 


where the function ffe satisfies the following equation 
(w—kv) fi, = 0 (II. 55) 


In the absence of external beams we can understand fiw as a fluctuating 
function satisfying relation (II.52). Previously we neglected this additional 
term, assuming that the amplitude of turbulent oscillations considerably 
exceeds the thermal level. With the inclusion of the term f¥o the equation for 
Oko Decomes 


4 
eko)Pro == 2 f fln dv (II. 56) 


that is, an additional noise source appears. Accordingly, in eqn. (11.50) 
additional terms appear taking into account the thermal fluctuations, and it 


assumes the form (see (213)) 


eka) = i E ZE E ; | 0o- kv) dv+ 
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The first term on the right hand side of this equation represents the 
thermal fluctuations. We can represent eqn. (II.51) symbolically in the form 


= = 2yI+q—al? (II. 58) 
where the first term describes the build-up of the oscillations with growth 
rate y, q represents the source due to thermal noise, and the non-linear term 
describes the interaction between the waves. Thus in the stationary turbulent 
state for a large growth rate, we may neglect q, and thus J = «/2y. On the 
other hand, when y is negative and not very small, we can neglect the non- 
linear term in eqn. (II.58) so that J = q/2|y|. In this case only thermal 
noise is present in the plasma. For y — 0 this noise level diverges and to 
determine J we must preserve the non-linear term. This is illustrated in 
Fig. 11 which shows the transition from thermal to turbulent fluctuations. 


y 
Fic. 11. Dependence of noise intensity on growth rate of small oscillations 


(c) Wave-particle Interaction 


Let us now return to the equation for the averaged function (11.38). The 
expression on the right hand side of this equation, which we shall again 
denote by S.f, represents the collisions between the particles and waves 
and can be written in the form 


eo 
Sp = —Im—— | kP 
ef moe fk ko Y) dk dw (II. 59) 


where the correlation function Pko is defined by 
Proel — w skk’) = (Pko fro (II. 60) 


Utilizing eqn. (11.46) for fke and adding to it the term f4,, related to the 
thermal fluctuations (see (II.53)), we obtain in the random phase approxima- 
tion 


Pv) = sinc 


Oo 6(@ — kY) + IkolKgko)f — 


— Iko f (K' Bio) (k8rro)(k'B-w oS +k Brera) KEko) f Heo Uk’ deo’ + 


Il. INTERACTION BETWEEN WAVES IN WEAK TURBULENCE 47 


Vk" wo" ko ' n ! t , 
+ ho | eE gua) (K'ero K gwo) Mire dk dow + 


i 
+ sg] Uko, k'olK Bko) CK" Sara dS iko kro” dk" da’ (I. 61) 


where k” = k-k’, w” = w— w’, e(k, œ) = e(k, œw + iv). 

In this expression the first term describes the slowing down of the particles 
due to the polarization of the medium and the radiation of longitudinal 
waves by the Cerenkov effect. The second term, which is linear in Ake, has 
already been considered earlier when treating the quasi-linear approximation, 
and the remaining terms in (II.61) arise from the non-linear interaction. 

Let us first consider the simple case of a stationary stable plasma where 
the quadratic terms in J in (11.57) and (II.61) can be neglected. Determining 
from the kinetic wave eqn. (II.57) the equilibrium intensity of the fluctuations 
iko and substituting the result into the second term in (II.61), we obtain 

mel 

4nefd(w —kv) ae Ov 2e? a ae 

(2n)*k?e%(k,@)  mo—kv+iv zle(k,w)|2k* | IIIA EV IGN 
(II. 62) 


Substituting this expression in (1I.59) and remembering that according to 
(11.43) 


Pye ¥) z 


e” = Ime = re [azo kv) d I 
=Ime=—~ 5 Ay (œw— kv) dv (II. 63) 


we obtain 


ef ~— 


2e* f d(kv—k’'v’) 6f(v') ~ ôf 
ae f kVa(k, kv)? {sk ay! f(v K) dk dy (II. 64) 

In the region KD > 1 (D being the Debye screening radius), the dielectric 
permeability can be set equal to unity, and the integral (11.64) may be shown 
to reduce to the collision term in the Landau form (where the integral has 
to be cut off at the upper limit for k ~ 1/po, where pọ is the minimum 
separation of the particles in the case of binary collisions). 

In the form (II.64) the collision term has been obtained in refs. (50), (51), 
(55) (see also (53, 66-69, 214, 215)). It represents both binary collisions 
with impact parameter less than D, and also the interaction through longi- 
tudinal waves. It has been shown by Davidov (52) that in a plasma approxi- 
mating to equilibrium, the contribution to the collision term due to the 
Langmuir waves is only about one order of magnitude smaller than the 
contribution due to the binary collisions. In a non-equilibrium plasma the 
contribution due to the oscillations may be considerably larger. For instance, 
Gorbunov and Silin (54) have shown that in a strongly non-isothermal 
plasma T,/T; < m,/m; the interaction due to ion acoustic waves predominates 
over the binary collisions. Another case where the interaction through the 
waves is larger than the collisional interaction has been considered in refs. 
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(44), (45); this is the interaction between the electrons and the cyclotron 
radiation in a strong magnetic field. 

As we approach instability the thermal noise level increases, and conse- 
quently in (II.61) the second term becomes more important, which in the 
quasi-linear approximation leads to the diffusion of the particles in velocity 
space. In this case, and also in the case of a weak instability of the plasma, 
it is more convenient to consider an initial-value problem, so that in the kinetic 
wave ean. (II.57), we have to take into account the term containing the time 
derivative of the oscillation intensity, and in the expression for the collision 
term S,, it is necessary to consider the adiabatic interaction of the particles 


oh 
with the waves, i.e. the term with the derivative z In other words, we 


arrive at the quasi-linear approximation considered in Section 3 of Chapter I. 

For large oscillation amplitudes, when the plasma goes over into a turbu- 
lent motion, it is necessary to consider the non-linear terms in the kinetic 
wave eqn. (11.57) and in the wave interaction term (11.59). Some examples 
of such processes will also be considered in Chapter IV. We shall see that in 
many cases low frequency oscillations can be excited in the plasma by a 
longitudinal or transverse electron current. When the phase velocity of these 
oscillations parallel to the magnetic field is smaller than the thermal velocity 
of the electrons, a resonant interaction may occur between the electrons and 
the waves, and for the electrons it is then sufficient to consider only the quasi- 
linear terms. For the ions, the thermal velocity of which may be considerably 
smaller than the longitudinal phase velocity of the waves, it is necessary to 
consider the non-linear terms. We then find that the energy and momentum 
of the electrons is transferred by a resonance mechanism to the waves, and 1s 
then absorbed by the ions due to the non-linear damping of the combination 
waves. Processes of this type may give rise to an anomalous resistivity of the 
plasma, which determines both the anomalous diffusion of the plasma across 
the magnetic field and its turbulent heating. 


Il. METHODS OF CONSIDERING STRONG TURBULENCE 


]. THE WEAK COUPLING APPROXIMATION 
(a) Wave Equations in Weak Coupling Conditions 


HITHERTO we have assumed any interaction between waves to be small and, 
strictly speaking, infinitely small. Let us now try to examine what occurs as 
the matrix element of interaction increases. We shall again consider the 
model eqn. (II.1) which we write in the form 


(w z Ok)Cko = Í Vko, ko Ck’ Ck-w', w= o’ dk" dw’ (ITI. 1) 


where the frequency œp is now complex. 

Suppose the matrix element of interaction Vko, ho, though remaining 
small, increases, approximating in order of magnitude to unity. The inter- 
action between the separate waves broadens out the frequency spectrum of the 
oscillations and in the limit of strong turbulence the dependence of Ixo on the 
frequency no longer bears any resemblance to a 6-function (see Fig. 12). 


Fic. 12. Spectral functions for weak (1) and strong (2) turbulence 


We cannot, therefore, use the kinetic equation in the form (11.12), but as 
long as the matrix element remains less than unity, we can use the weak 
coupling method as an approximation even for turbulence which is not weak. 

It should be noted that according to eqn. (III.1), a separate mode k, œw 
interacts only with two other modes k’ œ and k” œ”. Since the number of 
modes is large and the amplitude of the oscillations of an individual mode is 
determined by its interaction with all other modes, it might be thought that 
the interaction of the wave k, œw, with each separate wave k’, œw’, would be 
comparatively small even for V ~ 1. Moreover, as we have seen further 
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above, one of the principal non-linear effects is an additional damping, which 
is implicitly determined by the right hand side of (III.1). It is natural to 
isolate this part of the interaction earlier from the right hand side by writing 
equation (III.1) in the form 


(w — OK + Mo) Cro = Nko Cko +f Vko, k'o Cko’ Cirer osu dk’ dw’ (III. 2) 


The terms ko Cko Which we have added to both sides of this equation 
represent the part of the non-linear interaction proportional to Cko. In the 
right hand side of eqn. (III.2), from which we have thus removed the self- 
action (i.e. the damping) of each mode, only the input from the beat interac- 
tion of different modes remains. These inputs will be considered small, 
which is justified for V < 1. Accordingly we again put Cro = CK + CÈR 
where C) < Ci”. To determine the amplitude of the induced oscillations 
CŒ) we need only the non-linear term on the right hand side of eqn. (III.2) 
and we obtain 


CÈ = (0-01 + Mw) | owo COQdk'doo’ (I.3) 


We now multiply eqn. (III.2) by Cy, and average the result over the 
statistical ensemble, assuming that the amplitudes C{°? with different kw are 
statistically independent. We now substitute Cpo = C04 CH) into the non- 
linear term, and just as in the case of (II.6), this term reduces to a sum of 
three, two of which are proportional to ko, and a third to the integral of 
the product Iwe’ Ly”. Defining the quantity Nko so as to eliminate the terms 
proportional to Iko, we obtain the following system of two integral eqns. 
(111.4) and (II.5) 


lo — Oy + Mkeo| Iko ss 3 | [Pko wol Tro lke” dk’ doo’ (IIT. 4) 
Uko, k’w” Vk”o”, ko ; 
Neo = {> OS Tig dk den’ (III. 5) 
a" ‘= Wy FE Mra” 


where 
Uko, k'o’ = Vko, k’o’ F Vko, k”œo” 
We now introduce Sko = (W—-W_ +ko) +, and eliminating nko, obtain 


Iko = Skool” Í [Vko, wol lea Ikro” dk’ dw’ (IIT. 6) 
Sko = Sko — Sko Sko | Sera” Mee,n0" Pero" bo lye dk’ deo’ (IIL 7) 


where Sp, = (w— w) +. 

The quantities Sko and Sv, have a simple physical meaning which can be 
seen by adding a small external source fko to the right hand side of (III.2). 
By repeating the argument, it may be verified that Sko represents a Green’s 
function which describes the response of the turbulent medium to a small 
“ force ” fxo, While SRo represents the same Green’s function in the linear 
approximation. 
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In the theory of fluid turbulence, equations of the form of (III.6) and 
(111.7) were obtained by Kraichnan (18), and Wild (56) has shown that these 
equations can be obtained from partial summation of terms in perturbation 
theory. To obtain Kraichnan’s equation we write the Navier-Stokes equation 
in Fourier representation 


(œ + ivk?)Jüko — K Pro =| Geo K Wyo dk' daw’ = fro (III. 8) 


Where Pro 18 the pressure, fko tlie external force, v the coefficient of kinematic 
viscosity. Using the incompressibility condition ku,,, = 0, and assuming 
without loss of generality kf. = 0, we can eliminate the pressure from 
eqn. (II.8) and obtain 


, k 
(w + irk’õuko = fro + | (Uko K) furor -yz (kuy-)} dk'dw (III. 9) 


Repeating the preceding arguments and averaging the result over angles, 
assuming fko to be isotropic, we arrive at the following system of integral 
equations 


Tico = |Steol?Qeco +4| Stool” i ka(k, Kor Turor dk’ do’ (III. 10) 
k? b(k, k” yo 
a | TETTE 
where Jko is the spectral function defined by 
I 0(@ — w’)6(k—k’) = Cu, Uy a> 
Sko = (œ + ivk? +ga) t 
Fw ôl — 0')O(K—K’) = [fro fwo? 


and the matrix elements a(k, k”) and b(k, k”) are defined (see refs. (18, 56)) 
by the relations 


dk’ do’ (IIT. 11) 


nL (kk’)*(kk”)* — (kk’)(kk”) (k’k”) 
ak) = 1-2 Gast apapa | D 
(k na (k'k”) (kk”) 


(IIT. 13) 


b(k, k”) = ktk'2 seer? 

These equations represent a Fourier transform of Kraichnan’s equations. 
It has been shown by Kraichnan that they lead to an incorrect asymptotic 
behaviour of the spectral function for k > œ; instead of the well-known 
Kolomogorov law I,dk ~ k75’? dk we obtain the spectrum J,dk ~ k7~3!? dk. 
We shall show below that this occurs because in the weak coupling approxi- 
mation used here, the adiabatic interaction of the distant harmonics is not 
separated out. 


(b) Weak Coupling in Kinetics 
We shall show by the example of electron Langmuir oscillations how the 
weak coupling equations can be set up for strong turbulence ina collisionless 
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plasma. We take as the initial equations 


4 
Pro = — a | Su dv (III. 14) 


(OK fio =~ prok F ZL fiy E TE E ei ile 


m OV 


(Il. 15) 
We again represent fko and Eko in the form fx. +f, of + eV), where the 
superscript (1) denotes the part representing the induced oscillations. When 
these forms for fand @ are substituted in the non-linear terms, we shall just 
as before be able to extract “‘ self-action ” terms proportional respectively to 
fio and Eko, and it is again natural to separate out these terms. Because of 
the non-linearity these terms are not proportional to one another, or in 
other words their ratio is a random value. We therefore write the self-action 
terms as a linear combination of fko and Øko: NkoSko Eko Pkor Where Nko 1S 
an operator acting on the variable v, and čka is a function of velocity. As a 
result we obtain 


_4 
gi) = ed | a dy (111. 16) 


ay, SF _ e @ 


on df= E al k'o, A2, dk! do’ 
m OV moy 


(III. 17) 


From these equations we express fi,’ and {$} in terms of an integral of 
Qo, fir. and substitute the resulting expressions in the non-linear term of 
a (11.15). Then, multiplying eae (111.14), (111.15) by oy, and averaging 
over the random phases of A9, of, we obtain, dropping the superscript (0) 


4re 


ho =- 45 f Po dv (IIL. 18) 


0 
(w—kv) Pyro = e ko a x 
m Ov m OV 


f= f 0 
x f Gyr KP yoyo K’Praolyras) dk’ do't, EI x 
m OV 


k’ o „k” k” a) A kosko y + 1 r 
<| (9x. pien 7 o + ko rror t dk do’ + 


B -fk Gaver kN EEES P Tyee! do (I. 19) 


é*(ko) 
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where 
sce: 0 
Oko = (@— kv + nko) ae = 
m OV 
of 
4ne? or 4ne Exc 
êlk, œ) = 1+ —,; ee mL ee ae 
mk“ J w—kv+n,,, k w — kY + ko 


~ 4ne ~ / ~ M7 ar. 1 ~. 
Vko, k'o’ = k? al {(ðko K Grok f) T (ko K V(Gure’ k'f)} dy (TII. 20) 


It can be seen from (III.19) that in this approximation the non-linear term 
reduces to a sum of two components, of which the first is proportional to Po, 
and the second originates from averaging in the following form <o*o). 
Identifying these enn ae with Nko Pko and Cho Iko, gives 
5 | Kraal dk’ do (III. 21) 


Nko lko = 
m m ôv 


cko = ma nov] Gx a” KP -r-o dk’ dw’ + 


k” k” koy Vicon’, Kokoy y i 
+ Zf Grew” fy {aces vot a) vit vo} dk do’ (III. 21a) 


where Ko = Pikol Tuy: 
Considering these relations and substituting the expression for Pko from 
(IIT.19) in eqn. (III.18), we obtain the first integral equation in the form 


l&(k, o)|?I,.. = 4 | Fico, we Iwo Tvor dk’ deo" (III. 22) 


Equation (III.21) can be used as the second equation and the third is 
obtained by substituting the expression for čko (III.21a) in eqn. (III.20) 


£ 4ne f 4ne z z 
i(k, 0) = 14 | Ook fav + TS f EBra) Kuro) P-o Ak’ de! dv + 


gne [ Tko, wc” Vireo”, ko i 
2 ko Ixo dk’ dw’ (III. 23) 
If we write approximately py. = JkoKf, eqns. (III.21), (III.22), (111.23) 
constitute a complete system of equations for the three unknown quantities 
Nko, E(k@), and i... Since Nko is an Operator, these equations are symbolical, 
and to explain their meaning it is necessary to perform an expansion into a 
power series with respect to Iko. A series expansion in the amplitude of 
oscillations and the selective summing of this series is also necessary for the 
rigorous justification of the above equations, which we have obtained 
essentially from semi-intuitive considerations. (It has been shown by 
Mikhailovskii that in the case of a plasma in a strong magnetic field integral 
equations of the form (II1.21), (III.22) can in fact be obtained by selective 
summing, similar to Wild’s summation, which as a first approximation leads 
to Kraichnan’s equations.) 
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(c) Resonance and Adiabatic Wave Interaction 


In the arguments given above we implicitly allowed only resonant interac- 
tions between waves. However, the interaction of modes of very different 
wavelengths does not necessarily have to be of the resonance type. (Both 
here and below we are discussing only that part of the interaction which 
leads to secular terms, i.e. to a damping of the modes, and not the part whose 
effect is to shift their characteristic frequencies.) 

In order to analyse the character of the interaction rather more carefully, 
we shall again consider the model egn. (III.1). Suppose for simplicity that 
the frequency œp increases monotonically with k, and the ratio of the growth 
rate to the frequency, Yk/@p, is small. As we explained earlier, in a stationary 
turbulent state the damping of each mode is compensated by the input due 
to the beat interaction. The wave k, œ then exists for a time ~ y ‘ and 


: ; : : w ; ; az 
occupies a region in space of size ~ —-. During a period of the order y~}, 


yk 
an individual wave disappears completely and is replaced by another wave, 
which originates from the beat interaction and consequently cannot be cor- 
related with the first. Although the mode may initially be localised, it in 
general spreads out as time passes, to fill a region of space whose characteristic 
size L will be of the order of the distance through which the wave propagates 
do œw 
dk © yk 

Thus the state of turbulent motion of a continuous medium must be 
regarded as a system of many wave packets. For y/w < 1 these packets 
exist for a very long time and are very extended, so that one can describe them 
as waves which are almost completely unlocalised in space. But as y/œ 
increases, we must explicitly consider that the elements of the turbulent 
motion are not the Fourier components, but wave packets. In other words, 
for a finite y/w, nearby Fourier components can no longer be considered 
weakly correlated. 

We can put the matter in another way as follows. Consider the region 
of wave numbers of order k. Since the waves in this region exist for a period 
t ~ yk, when we average the quadratic terms over time it is sufficient to 
cover times of the order of yg *. All slower oscillations, with frequencies 
w' S Yk, can be considered fixed during this period of time and need not be 
averaged. These modes can be regarded as forming an inhomogeneity of the 
background. For short waves, this inhomogeneity may be considered in the 
quasi-classical approximation discussed earlier. As we have seen, the wave 
packets are deformed and the wavelength changed as they propagate through 
such an inhomogeneity. As a result of these changes, the wave packets 
move about in wave number space, which leads to a strong correlation of 
nearby Fourier components, which now describe essentially one and the same 
wave packet. 

Thus the interaction between the wave k, w, and the slow wave k’w’ leads 
to an adiabatic variation of the wave vector k and of the frequency w of the 


during its lifetime, i.e. L ~ y7! 


II. METHODS OF CONSIDERING STRONG TURBULENCE 55 


wave packet under consideration. This interaction cannot be considered 
as the resonant input to the wave k, w from the nearby wave k” = k—k’, 
w” = w—w’, because the Fourier components Cre and Cko” describe the 
same wave packet and cannot be considered independently of one another. 

Thus, for finite y/o we must make a distinction between resonant and 
adiabatic wave interactions. To bring in this distinction we must change 
from the Fourier representation to a representation in terms of the wave 
packets. For each packet with a mean frequency w and mean wave vector k 
the region of integration with respect to k’, œ in the non-linear term is 
conveniently broken up into three parts: (1) the principal region where 
k' ~ k, œw ~ œ, (2) the long wave region where w’ ~ y and (3) the short 
wave region where y’ ~ œ (y/@ is assumed small). 

For small y/œ the spreading of the wave packets in the principal region 
can be neglected. In this case we can use eqn. (III.1) for the Fourier com- 
ponents. In the regions (2) and (3) the quasi-classical approximation can be 
used, i.e. an expansion in k’/k, w'/@ and k/k’, w/w’ respectively. Limiting 
the treatment to the first (or more precisely to the zero) approximation, we 
obtain 


(O= 4)Cro Z | Veo wor Cro Cera dk! dor’ + 
(1) 
Cee | (Vico, wo + View, wo”) Crea AK’ do’ + 
(2) 
mn | Vices, walCwal? UK deo’ (III. 24) 
(3) 

where for simplicity we have omitted the exponent exp (— iw't +ik'r) in the 
integral over the region (2), assuming Cw to be a slowly varying function of 
r, f. 

In (III.24) the last integral term can be omitted, because when we divided 
the equations for C into two—for the averaged function and the oscillating 
part—we should have included it in the equation for the averaged function. 
The second integral on the right hand side of (III.24) is proportional to 
Cko and consequently it leads only to a shift of the characteristic frequency. 
In the region (1) this shift can be considered approximately constant. If the 
matrix element Vko, xwa depends only on the frequency difference w—a@’, 
or is independent of w, w’, which is often the case, the frequency shift can be 
eliminated by changing the variable appropriately. The final result is that 
we recover eqn. (III.1), with the difference that the range of integration in the 
non-linear term is restricted to cover the region (1). It follows that separate 
consideration of the adiabatic interaction 1s important only for wave numbers 
sufficiently large that the lower limit of integration of region (1) is greater than 
the main turbulence scale ko, @o. 

It is now obvious why the simple weak coupling approximation led 
Kraichnan to an incorrect spectrum: this approximation over-estimates the 
part played by the large-scale fluctuations, which is in fact no more than the 
convection of higher modes which are deformed adiabatically in the process. 
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If we make the corresponding change in the weak coupling eqns. (IH. 10), 
(III.11) by setting the lower limit of the integrations over k’, œ equal to 
ék, Ew, say, respectively, where č is a small constant number (say č ~ 4), 
these equations become self-modelling and lead to the correct spectral form 
Idk ~ k~* dk. However, it is obviously sufficient to use the simple weak 
coupling equations for modes in the region of the main scale. 


(d) Improved Weak Coupling Approximation 


To improve our understanding of the interaction between modes in a 
turbulent medium, we shall here give explicitly a series expansion of the second 
term on the right hand side of eqn. (III.24), to take into account the adiabatic 
distortion of the modes under consideration by the long wave perturbations. 
Assuming for simplicity the matrix element to be independent of the fre- 
quency, we obtain 


(O— Oy — O) Cyan X | Vie Creo Cr-r, o-o Tk’ do" — 


(1) 


ac ôC 
-k | kvi Cro Ak’ doo’ — 9 Í "Dyn? Cra AK! deo! + 
ak 0w 
(2) Sa 
F 2 | Virro Ckror?s dk’ do’ dk" da" (II. 25) 
(3) 
where 
Vew = Viw t Vig k-k” w, = Í Viw Cro AK’ do’ (III. 26) 


(2) 

The last term in (III.25) has been written as an average over high frequency 
oscillations, denoted by <...>,, which is permissible since only the secular 
effect of the oscillations in the high frequency region (3) can be important 
for the modes under consideration. 

The second and third terms in the right hand side of (III.25) are propor- 

; 0@, 0M, : ‘ y 

tional to E= and EE They describe the distortion of the wave packet k, œ, 
due to weak spatial inhomogeneity and its slow variation in time. As we 
have said, this effect is the same as the distortion of the wave packets in an 
inhomogeneous medium which we considered earlier (see Section I.3b). 

In eqn. (III.25) it is convenient to change the variable to the relative 
frequencies v = œw— w, with the corresponding amplitudes denoted by C,,. 
It is easily seen with this transformation that the third integral on the right 


l , „ÔC ô 
hand side of (111.25) vanishes, since 1 KY = p-s Cyy. The re- 
Ot Ot dv 


maining terms retain their previous form, apart from the change of the 
frequencies œw, œw’ to y, v’. 

Using the argument of the previous section, we again separate out from 
the right hand side of (III.25) the quantity nky C,,, proportional to C,,, and 
write the amplitude in the form C$? + CŒ, where C{? describes the induced 
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oscillations. According to (III.25) the amplitude C{}) can be written as a sum 
of two components, Cii2+C where the first is due to resonance inter- 
action in the region (1) and the second given by 

(0) 


ace? 
Css ee | Ky CO, dk’ deo’ (L27 


ôk 
(2) 
originates from the distortion of the wave packets at the long wavelength 
inhomogeneities. 

Multiplying by Cf, and averaging over the random phases of the ampli- 
tudes C{°’, we obtain, in addition to the contribution due to the resonance 
interaction, terms of two new types. Two terms of the first type originate 
from the second on the right hand side of eqn. (III.25). Their sum is pro- 


. (1)* ac (0)* Cin ; . . 
portional to Cku “ak + Ci, k: These terms describe the diffusion 


of the wave packet k, vin momentum space due to the distortion introduced 
by the long wave length modes. The value of the corresponding diffusion 
coefficient will be proportional to the integral over region (2) of the quantity 
k'*Iy. Since the wave number k’ in this region is very small, this diffusion 
can be neglected.t Additional terms of the second type originate from the 
last integral in eqn. (III.25), where instead of Cwe or Cyra We must put CHP. 
Since Cy, and Cro” enter perfectly symmetrically, it is sufficient to replace 
only one of these, say Cy--, and then to double the result. Now for region (3) 
the long wave region (a region (2’) as it were) is represented by our region (1), 
and using this the corresponding contribution may be obtained without 
difficulty. It is proportional to J», and should therefore be included in the 
equation for gy. We finally obtain 


Ay 


v ” n" I ry! r ? al ty! 
oe | Pan Pere Ke ae dy! + | ew Pek Kf er gy (LIT. 28) 
y 


v” — Dy FI ry Op tiyy ok’ 
(1) (3) 
i,, = 418,,|? Í lonl ever dk' dv’ (III. 29) 
(1) 
where k” = k-k’, wo” = wm =w; Shy = (v-o tik) t 


We have marked the intensity I,, and the Green function Sy With a 
tilde to emphasize again that these values are calculated in a system of 
co-ordinates moving with the long wave pulsations. The true functions Iko» 
Sko are defined by the relations 


Iko = EN Sko = OS (II. 30) 
which are averaged over the long wave pulsations (we recall that œ is still a 
random value). 


t This assertion is no longer valid where for the modes under consideration there is no 
resonance interaction: in that case the distortion of the wave packets due to the long 
wavelength modes becomes the dominant effect. Such a situation occurs for instance, 
in the short wave (viscous) region of ordinary turbulence (see (57) (58)). 


P.T. 3 
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According to eqn. (III.26), œw, can be represented approximately by 
W, = VC, Where C; = Í Co dkdæw is the amplitude of the whole set of 


(2) 
long wavelength modes represented by region (2). C, is a random quantity 
given by the sum of a large number of individually random and weakly 
correlated amplitudes; thus the distribution of C, and consequently also 
that of œ, is Gaussian. Thus the relations (III.30) can be put in the form 


-22 1 or 
dies = Ea (S POT pri dwi, Sio — —_ e wo? Spani da, (III. 31) 
NWO TO 
where 
aoe = 2?) = 2\dyu|? Í Iva dk' doo’ (I. 32) 


(2) 

The calculation of Iko, Sko is of course only necessary if we are in fact 
interested in the temporal correlation of the fluctuations. In a number of 
problems it is quite sufficient to know only , = f iodo = f Ikydv, and we 
need then only calculate the relative values /,, and Suy- 

Equations (III.28) and (III.29) together with the relations (II1.30) con- 
stitute the equations of an improved weak coupling approximation. The main 
contribution to the damping of the waves is supplied by the resonance region, 
but, according to eqn. (III.27) a small additional damping arises in the short 
wave region (3). According to eqn. (III.29), the intensity of the oscillations 
I, is determined by the resonance interaction, but there is in addition a small 
contribution which we omitted, of the type of a diffusion in velocity space, 
which arises from the interaction with Jong wave fluctuations in region (2) 
and also, finally, an altogether insignificant addition may also be supplied 
by small occasional impulses due to the short wave oscillations of region (3). 

If we now reduce the coupling between the oscillations, we see that region 
(1) expands and in the limit y/w — 0 it occupies the whole of k space apart 
from isolated regions where the damping is not weak, and the perturbations 
can reach equilibrium with the remaining fluctuations fairly rapidly. In a 
plasma, there is a region of this type at short wavelengths KD > 1, where all 
oscillations are damped. This region in particular defines the binary collision 
term. Because of the weak (logarithmic) dependence of the collisional 
term (II.64) on the upper cut-off parameter, it is possible to write it in- 
advertently in the investigation of any given set of oscillations. In fact short 
wave fluctuations follow the oscillations adiabatically, although from the 
point of view of the quasi-linear approach (see Section II.2(b)) they must be 
considered as equivalent to all other oscillations. 


2. PHENOMENOLOGICAL APPROACH TO THE DESCRIPTION OF A STRONG 
‘TURBULENCE 


As we have shown, a quantitative theory only exists for the limit of weak 
turbulence where a kinetic equation for the waves can be set up. For this 
theory to be applicable the interaction between the waves must satisfy 
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several important limitations: the growth rate of small perturbations must be 
small compared with the frequency, the dispersion relation must be such 
that the laws of conservation of momentum and energy can be satisfied simul- 
taneously at discrete surfaces in momentum space, the interaction between 
the waves must provide additional dissipation in each region in momentum 
space, and so forth. In practice these limitations must fairly often be in- 
fringed, and then strong turbulence develops. 

Starting with the kinetic wave equation, we have shown above that the 
transition to strong turbulence leads to integral equations in the weak 
coupling approximation, in which the resonant and adiabatic interaction 
must be separated. At present we have no rigorous method of performing 
this separation and of reducing the integral equations, and the description 
given by this theory is inevitably only approximate. We have as yet no indica- 
tion of the accuracy of this approximation. 

However, in describing strong turbulence in a plasma we can use an 
analogy with ordinary turbulence. Here the principal results have been 
obtained from a purely phenomenological approach. It is natural, therefore, 
to use this approach in plasma turbulence theory. 

In the next section we summarize some of the results of the phenomeno- 
logical description of ordinary turbulence, and will show later how similar 
concepts may be used when considering plasma turbulence. 


(a) The Turbulent Jet 


Using the phenomenological approach, the mean velocity profile across 
a turbulent flow may be obtained without discussing the spectral function. 
We use the mixing length concept introduced by Prandtl (59), and consider 
one of the simplest problems, the velocity profile in a submerged turbulent 
jet. Such a jet is formed when a liquid or gas flows out of an orifice into a 
medium of the same state. When it emerges from the orifice the velocity of 
the fluid is nearly constant over the cross-section; as the distance from the 
orifice increases, the velocity profile becomes more and more deformed until, 
at a fairly large distance the jet becomes self-modelling, i.e. the velocity 
profile retains the same form but its width increases in proportion to the 
distance from the orifice. 

For simplicity we consider a flat jet emerging from a slit. The mean 
velocity depends only on the co-ordinates x parallel to the jet and y transverse 
to it. The mean pressure in the free jet can be assumed constant, and the 
longitudinal component of the averaged equation of motion can be written 
in the form: 


ĝu ðu 1 Oty 


—+v-—= II. 33 

“ax Oy p oy ( ) 

where u and v represent the x- and y-components of the mean velocity 
respectively, p the density of the fluid, t,, = —p<u'v’> is the tangential 


Reynolds stress arising from the transfer of longitudinal fluctuating momen- 
tum pu’ by the transverse fluctuations v’. The value of v’ is of the same order 


60 PLASMA TURBULENCE 


as that of u’, since they are related by the continuity equation. To obtain an 
estimate of u’, Prandtl assumed that the velocity fluctuations arose essentially 
from the conservation of the mean longitudinal momentum of a fluid element 
during its transverse displacement. This gives the estimate 


Ou 
oy 
where / is the mean transverse displacement of the fluid elements, which is 
the so-called “ mixing length ”. In a free jet the mixing length can be assumed 
to be approximately constant over the cross-section and proportional to the 


thickness of the jet, i.e. Z = cx, where x is the distance from the origin of the 
jet, and c is a constant. Thus we obtain approximately 


S 
Ile 
em 


(III. 34) 


Ou |ðu 
Tay = + pc?x? ay Bs (III. 35) 
and eqn. (III.33) becomes 

dus du du 0*u 
pas — = +2c?x? -L II. 36 
ay Se a) ( ) 

This equation and the continuity equation 
ðu ðv 

—+—-=0 HI. 37 
Ox i oy l ) 


are a complete set of equations, and determine the mean velocity profile in a 
turbulent jet. This system contains only one unknown constant c, which can 


be determined by comparing the theoretically calculated velocity profile with 
that measured experimentally. 


0 
Fic. 13. Velocity profile in turbulent jet 


Equations (III.36) and (III.37) were solved by Tollmien (60). Figure 13 
shows a comparison between Tollmien’s profile and the experimentally 
measured profile (this figure has been taken from ref. (61)). With an appro- 
priate choice of the constant c the theoretical profile agrees very well with 
that determined experimentally. The mixing length turns out to be an order 
smaller than the half thickness of the jet. 

While other authors have proposed different hypotheses relating to the 
character of the momentum transfer across the jet (see, for instance, ref. (61)), 
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the general validity of Prandtl’s description of the process in terms of mixing 
length theory has remained unchallenged. 


(b) The Kolmogorov-Obukhov Law 


Let us now go on to the problem of the fluctuation spectrum in a turbulent 
fluid. It is natural to limit the discussion to scales which are considerably 
smaller than the principal scale, so that it is possible to refer to local proper- 
ties of the turbulent flow. The theory of such a locally isotropic turbulence 
was developed by Kolmogorov and Obukhov (63, 64) and later by Heisen- 
berg (65). We shall discuss here only the principal results of this theory, and 
this in considerably simplified form. 

When the Reynolds number is very large, the fluid viscosity is negligible 
over a wide range of scales of the turbulent motion. Over this range the energy 
cannot be dissipated, but can only diffuse in wave-number space. This 
principal region of the spectrum of the turbulent fluctuations is called the 
inertial range. On the other hand, in the region of lengths, less than some 
small scale Ap—the so-called internal scale of the turbulence—the viscosity is 
dominant. In this region the energy of the motion is dissipated into heat. 

According to Kolmogorov, in the inertial range of the spectrum a quasi- 
equilibrium is established, in which there is a constant energy flux e through 
the spectrum into the short wave region À < Ag, where dissipation occurs. 
The value of this flux € determines the local properties of the turbulence. 
This hypothesis is equivalent to the natural assumption that the energy 
transfer between modes is of a resonance character, in which, as we have seen, 
the energy of a mode of scale A can be transferred only to modes with nearly 
the same scale. Thus a portion of energy handed down from large to smaller 
scales must pass through the entire range of scales of motion almost to dp. 
We may assert, therefore, that for each scale 4 the value of € is determined 
only by the fluctuation level at this scale, that is by the value of the spectral 
function for k = 27/2. In other words, e must be expressible in terms of k 
and J, only. The only dimensionally correct combination is 


e~ VI, k’ pl, k (III. 38) 
since J,k? has the dimension t~*, and kJ, the dimension of a square of the 
velocity. From eqn. (III.38) the well known 5/3 law of Kolomogorov- 
Obukhov follows immediately: 


£ 2/3 
I,dk = C (=) k75! dk (III. 39) 
p 


where C is a universal dimensionless constant. 

This law agrees well with experimental data. It has recently been shown 
experimentally that ¢ is itself a random quantity undergoing large fluctuations, 
and in accordance with this Kolomogorov and Obukhov have made some 
changes in the theory, but we need not go into details of these in this review. 

Thus the 5/3 law can be obtained from very general considerations of 
dimensional theory and the hypothesis of the absence of strong interaction 
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between pulsations of very different scales. To obtain the spectrum over the 
whole range of wave-numbers, including the viscous region A—Ay, more 
precise ideas about the character of the interaction between the pulsations 
are required. So far a large number of various hypotheses of a phenomeno- 
logical character has been proposed to represent this interaction. The 
spectra which can be obtained from these hypotheses are very similar to each 
other (see ref. (62)), and in the inertial range all phenomenological theories, 
except those which are obviously wrong, lead to the 5/3 law. 


(c) Wind Waves 


As an example very similar to phenomena which occur in a plasma, let 
us discuss the development and interaction of waves on a high sea. The 
dispersion relation for gravitational waves is defined by the relation 


wo =~ gk (IIT. 40) 


where g is the gravitational acceleration. 

This dispersion relation is of non-decay type. This means that for small 
amplitude of the oscillations the principal part must be played by four wave 
processes involving the conversion of two waves into two other waves. 
Using perturbation theory we can obtain a kinetic wave equation for this 
type of interaction, but this equation is very complex and its solution is 
difficult. A much simpler method is the phenomenological approach in which 
the general form of the spectral function is established from physical con- 
siderations (refs. 70, 71). 

Suppose that at some initial moment a wind starts to blow with a constant 
velocity above the surface of the sea. At first, when the surface of the sea is 
almost unperturbed, the principal mechanism governing the development of 
the waves will be their build-up by resonance with turbulent pulsations of the 
atmospheric pressure. The energy of each mode increases linearly with time. 
But as soon as the slope of the wavy surface attains some finite value, an 
instability mechanism will begin to play the principal part. The physical 
origin of this instability can be explained by reference to Fig. 14. When the 


Fic. 14. Excitation of waves on water surface 


slope of the wave is large enough, the air pressure will be slightly larger on the 
side A facing the wind than on the side B screened from the wind. This 
difference in pressure, due in the last analysis to the eddying of the wind 
above the region B, leads to an effective resonant transfer of energy from the 
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wind to the wave, since where the surface is rising the pressure is slightly 
smaller, and where it is falling, it is slightly larger than the average pressure. 

The instability mechanism, which is in fact non-linear, leads to an 
exponential increase of the amplitude of the wave. But since it operates only 
for a considerable slope, which for a given amplitude will be larger the smaller 
the wavelength, only perturbations with wavelength smaller than some 
critical wavelength will increase exponentially at any given time. In going 
towards longer wavelength the spectral function must therefore decrease very 
rapidly with A, probably more or less exponentially. 

In the short wave region of the spectrum the amplitude in the steady 
state is determined by the non-linear interaction between the modes. The 
simplest sufficiently strong interaction between the waves is their collapse, as 
a result of which the ‘“‘ white crests ” appear at the surface of the sea. This 
interaction limits the amplitude of the oscillations to such an extent that 
pointed “ tips ” appear on the waves. These “tips” correspond to dis- 
continuities of the gradient so that the second derivative of the surface 
displacement € must contain 6-functions, and for large k values the Fourier 
transform of č must decrease as č ~ k~?. It follows, therefore, that the 
spectral function of the vertical displacement č for large k values in the 
equilibrium region must behave as I,dk ~ k~*dk where dk = dk,dk, ~ kdk. 
Transforming from the variable k to œ we obtain I,,dwa ~ dw/w°. This 
spectrum has been obtained by Phillips (70). It agrees well with experimental 
data relating to the short wave region of the spectrum. The total spectrum 
of the oscillations of the surface, measured experimentally, can be repre- 
sented satisfactorily by the function 


I, = ag? > e7200? (III. 41) 


where « and wo, depend only on the time. This spectrum agrees perfectly 
with the above physical argument. 


IV. TURBULENCE IN A PLASMA 


1. NON-LINEAR DAMPING OF LANGMUIR WAVES 


THE discussion of specific turbulent processes in a plasma is best begun 
with the simplest case, that of Langmuir oscillations. In particular we shall 
consider here in some detail the non-linear damping of waves due to the 
resonant interaction. Assuming that the mean wave number of the wave 
packet under consideration is appreciably smaller than the reciprocal Debye 
radius D~* we can neglect the linear Landau damping and retain only the 
second and third terms in the kinetic wave eqn. (II.51) (the last term dis- 
appears since Langmuir oscillations are of the non-decay type). 
In the matrix element 


4re 
Dowo =- Ge | (rok Grok + Cok Eek) f} av (IV. D 


we keep only the first term, since for a small frequency difference w” = w—@’ 

the quantity 8ko” > Zkos Zero. By integrating by parts we replace gk’ by 
e kk e kk’ 

tk Se eC aoe o” and the remaining integral can be expressed 


in terms of the dielectric permeability (11.43). Thus, in zero approximation 
(0) _ e (kk’) ki 


Vio, ka” = 
ako” = m w? kè 


[e(k”, o”)—1] (IV. 2) 


In this expression unity can be neglected compared with g, since for 
2 


w < kv, e(kq@) = re > 1. Thus e(k”w”) cancels in the third term in (11.43). 


e . 
Replacing k’g,. by its approximate value — Tuk? =~ mw iB the 


first integral, it is easily seen that this integral is proportional to Vko”, kw, 
and the two integral terms compensate one another exactly. This means 
that we must consider the thermal corrections which we initially neglected. 

To determine which small corrections we must consider in particular, 
let us examine the matrix element 


4ne 


Oko” ko = k"? | US" k’) (Eko k)f+ (Lure k) (g_y -0 kf} dv (IV. 3) 


The main contribution to this integral is supplied by the region where 
Zkorf—> 0, ie. where xo X wo Therefore the term with the second 
65 34 
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derivative of f with respect to the velocity, arising from the operation of gwo” 
on f, is very small. Neglecting this term we can assume that the operator 
urn” acts only on gy and g_,, -w Performing the differentiation and keep- 
ing only the largest term, we have 


4ne e kk’ a 
Porno E ps os flok = eo) AV. 4) 


This shows that it is aaa to consider d real correction terms in 
the matrix element Vko, 47a", Since the corresponding contribution will be 
proportional to 1—7! (k’w”), and its imaginary part of the order of 
|e(k”w”)|~* < 1 and thus extremely small. It is sufficient to consider only the 
small terms of the first and second order of smallness with respect to kv/œ 
in the operators Zxo in the first integral term in eqn. (II.43) and in the matrix 
element Vko, ca”. Collecting all these small terms and writing E y = we 


finally obtain 
Al, mo 
ôt  2mkên?o’ 
Since for Langmuir waves œ” ~ @,(kD)* < @, kv = k'v+o" x k'v 
and with an accuracy up to higher order than the second in kv/œw (kv)? can be 
replaced by (kv) (k’v). In this approximation the factor multiplying Jx 
inside the integral in the expression (IV.5) is anti-symmetric under the 


interchange of k and k’. It follows, therefore, that in this approximation 
Langmuir waves are scattered without energy loss (9) 


r, 
— | k*I,dk =0 
=| 


| (kk’)?(k)?(0"—kwyk’ dvI, dk’  (IV.5) 


According to (IV.5) the characteristic reciprocal diffusion time of the waves in 
k-space is of the order 


T 


k21, dk (IV. 6) 


where kọ is the mean value of re wave daa The damping of the waves 
is described by small terms, not considered here, and the corresponding 
damping decrement is of the order of (kD)*v. 

The energy damping rate can be determined from the following considera- 
tions. By integration by parts it can easily be shown that 


k"? 
Uko, k"o"” = T Vu" "ko (IV. 7) 


and we can see from the expression (IV.3) for the matrix element v4-4” ko that 
if the exponentially small remainders at the points œ = kv and w’ = k’v are 
neglected, and only the remainder at the point œ” = k’v retained, the ex- 
change of k and k’ is equivalent to taking the complex conjugate. In other 
words, the matrix element satisfies the symmetry condition 


Vy" @" ko = Un no", k'o’ (IV. 8) 
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From this relation and relation (IV.7) it can easily be seen that the imagi- 
nary part of the coefficient of Jw, in the expression below the integral in the 
third term on the right hand side of eqn. (11.51), is simply multiplied by 
~k?/k’? when k and k’ are interchanged. The second term in eqn. (II.51) 
has a similar symmetry. Making use of this symmetry we find that in the 


absence of decay the non-linear terms in eqn. (11.51) conserve the “ wave 
2 , 


number ” N, = é,/@, = g7 k zæ In other words, during the scattering 


of Langmuir waves on the electrons, the wave energy decreases in proportion 
to the frequency 


2.) (IV. 9) 


é,/@, = const, i.e. Ek =F, (1 F 5 
2MW5 


where the constant &,, represents the energy for k-0, i.e. t> 00. From 
eqn. (IV.9) the rate of damping of the energy can be determined from the 
rate of decrease of the wave number. 

The result so obtained admits of an explicit quantum mechanical inter- 
pretation. This process represents a scattering of plasmons at electrons. 
The condition œ” = k”v required by the 6-function represents in fact only 
conservation of energy: 


ha, — ha, + Apy = 0 
where Ap = — hk” is the momentum transferred to the electron. When the 


ð 
electron energy states are populated normally, with “a < 0, each plasmon 


loses energy on the average during scattering, so that its frequency and hence 
its wave number decrease. Since the frequency of a Langmuir plasmon 
depends only very weakly on the momentum k, the variation of its energy 
during scattering will be considerably smaller (by a factor k? D’) than the 
variation of momentum. For an inverted population of the electron states 
(i.e. for a “ negative ”? temperature) scattering is associated with an increase 
in the energy of the plasmon (74). 

For small amplitude of oscillation where higher wave interaction pro- 
cesses can be neglected, the scattering of the waves at the particles considered 
above is the main process. In contrast to conventional hydrodynamic 
turbulence where the interaction of the modes leads to a cascade of the energy 
towards large k, in the case of Langmuir oscillations the energy flux is 
directed towards the region of small k, where the linear damping is expo- 
nentially small. 

It can be seen from (IV.9) that for small values of kọD the energy of 
the Langmuir oscillations changes only slightly for £ + oo. In other words, 
single scattering processes of Langmuir waves on the electrons do not lead 
to their total relaxation in a homogeneous plasma. But when we consider an 

Ol, 0m, Ol, 
inhomogeneous plasma, we must add terms of the form Uz and =a 
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to the kinetic equation for the waves. The second term, which is the more 
important, describes the energy flux in k-space due to the distortion of the 
wave packets. At small k, the corresponding rate of change in the mean 


Wo 
Jy 7 K@o becomes comparable 
x 


to vky and kọ ceases to vary with time. pee -i (IV.6) ‘y can write the 
= (kD) =kD, From 
this we find that € = —v (ky D) E = ~ Kk D2? Since ~ equilibrium 


wave vector associated with this flux ko = 


condition for this 1 


T 
value kgD = G23) depends only weakly on & and x, the damping 


time of the Langmuir waves in an inhomogeneous plasma is to order of 
magnitude given by (Dkw) *. The time necessary for a deformation of 
the wave packet such that the wave number k attains the value D~* will be 
of the same order of magnitude; at this wave number Landau damping 
becomes important. In a homogeneous plasma, where this mechanism does 
not operate, the relaxation of the Langmuir waves is determined by the slower 
4-wave interaction processes and by the scattering of the waves at the ion 
fluctuations. 


2. EXCITATION OF ION OSCILLATIONS BY AN ELECTRON CURRENT 
(a) Excitation of Ion-sound Waves 


Let us consider the problem of the excitation of ion-acoustic waves by an 
electron current (72). Suppose that the ion temperature 7; is considerably 
smaller than the electron temperature T,. Such a situation arises for instance 
in a weakly ionised plasma where the ions lose energy in collisions with cold 
neutral gas atoms, or in non-stationary conditions where the ions are not 
heated by the electrons. The results of the stationary problem considered 
here can be used for the non-stationary case, since the time necessary for the 
establishment of the stationary oscillation spectrum is relatively small. 

For T; < T, ion oscillations are excited when the directed (current) 


velocity of the electrons u exceeds the velocity of sound c, -JZ — (see (75)- 


(77)). As u increases above c, the excitation of oblique waves becomes 
possible as well as of those propagating parallel to the current, and for 
u > c, the cone of permissible directions of the unstable waves opens up, so 
that practically all waves having a positive projection of the phase velocity in 
the direction of u are unstable. 

Two waves with wave numbers of about the same size and propagating 
at an angle to one another set up beats with small phase velocity since in this 
case œ” = w—w’ X 0 and k” #0. Since the decay of ion-acoustic waves 
is prohibited, the non-linear Landau damping at the beats makes the main 
contribution to the non-linear interaction of the waves. This damping limits 
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the amplitude of the waves, and in stationary conditions some oscillation 
spectrum independent of time must be established. 

The non-linear damping of ion-acoustic waves on the ions resembles 
the damping of Langmuir waves on the electrons. The considerations of 
the preceding paragraph can therefore be largely repeated. We can again 
omit the last term on the right hand side of the kinetic wave eqn. (II.51), 
because decay is prohibited. Moreover, the non-linear damping of the waves 
on the electrons can again be neglected because it will be described by the 
same terms as for Langmuir oscillations, except that m, is replaced by m,. 
As in the case of the electron oscillations, the non-linear terms in zero 


fee T; 
approximation with respect to the small parameter kv/œ ~ J T cancel each 


other out. The first order correction, which is in fact of the same order of 
magnitude, also vanishes and the next two terms of the expansion have to be 
considered. 

In zero approximation the matrix elements are determined only by the 
ions and are given by the following relation: 


ki"? 4ne?(kk’) 2kk’ aw" ôf: 
(0) (0 mar. n t n NNT 
Uko, ko” = k2 Wr = — ape? | + k? 2G) MES- av (w —k v) ! dy 
(IV. 10) 
i.e. they are proportional to the integral | k” To "—k"y)~! dy. 
The dielectric permeability also contains this aad 
Ane of, 4ne*n 
k”, M 1 k” = "n" k"! d : 
e(k”, œ”) T irf Oy "(œ@” —k”v) + dv + KT, (IV. 11) 


In this case the integral term in (IV.11) is considerably (of the order of 
T,/T; times) larger than unity, which may be neglected. Thus in the second 
integral term of the wave eqn. (II.50) it is sufficient to consider the corrections 
of the second and third orders of smallness in the imaginary parts of the 
matrix elements — 


kk’ Of; 
IM tko, ko” = TZ E in Veg” ko = a lee (@— ky? ròl" = k'y)k" 2 dy (IV. 12) 


These corrections can be obtained without difficulty and the kinetic 
equation for ion-acoustic waves reduces to the following form: 


1 de’ Gl 991 kko” 
LeOn a men = et f i-s a 0 + 


n a)" 3 12 
+ (3-4 rs =e 7} alo "yk" Dry dvdk’ (IV. 13) 


k?a@/] @ kK 


4ne?n l ; 
where 22 = ——., J, is the spectrum of the potential and g’ = Ree. 
m 


i 
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For T,/T, > 0 when the main contribution to the integral comes from 
the region wm” ~ 0, this equation takes the form 


Q2 61, OÉTI, 


m?w’ w k?n 


f (kk’)?{[kk’k”? — (kk")(k’k")J5'(o") Ly dk’ 


(IV. 14) 


where 6’ denotes the derivative of the 6-function and œ” = c,|k—k’|. 

For T, < T, the linear Landau damping on the ions can be neglected 
and there remains in Im e only the contribution due to the electrons which 
leads to the build-up of the oscillations. Including the wave damping due to 
collisions between ions and neutral gas atoms, and supposing that the electrons 
have a Maxwellian velocity distribution displaced by a velocity u along the z 
axis, we obtain the TERS en for Ime: 


p% Q2 1 
—k*I ay 0 — kv) dv — — 
me = (w—kv) dv de? oz, 
Qf 1 
= 2 (u cos 8@—c,) — — (IV. 15) 
E 4c? wr; 
x... Ale 5 4Ane*n 1 | T: 
where v2 = E wg = oe oe the mean collision frequency of 
i 


e 

the ions with the neutral gas atoms, and cos 0 = k,/k. Thus the collisional 
damping decreases as œw increases. For a sufficiently large q; it is appreciable 
only at very small w and its effect may be included by simply cutting off 
the spectral function towards low frequencies at the frequency at which 
Ime changes its sign, so that the waves damp. In this approximation the 
quantity k*Ime in the principal region of wave numbers can be considered 
independent of k. It follows immediately from (IV.14) that J, ~ 1/k?. We 
can demonstrate this in more detail. Suppose for simplicity that the waves 
are built up only in a relatively narrow cone of directions within an angle 6. 
(Such a situation arises, for instance, when the plasma is contained in a 
cylindrical tube which absorbs waves which propagate in the transverse 
direction.) In this case k” ~ Ook, and consequently the quantity kk” ~ k”? 
can be neglected compared with kk”. Transforming to k-space and using 
spherical co-ordinates, and integrating the 6’ term in (IV.14) by parts, we 
obtain 


Q2 ðI Qs m2 T, al ee 
Soop aoe, hye tf k20? A do +3k | Ohrada | (IV. 16) 


where 6, is the angle between k and k’ and dQ’ is the solid angle element. 
Integrating this equation with respect to Q assuming <0,>? = 67, and 
denoting by /, the integral of J, over Q, we obtain 


Ol, Ol, 
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where y = — 302 (Im e} is the mean growth rate of the small perturbations, 
QZ T, 02 
and A = 20.T?n This equation can be written in more explicit form: 
ol, 1 @ 
me ay (Ak Ti) = = 2y; I — AktI? (IV. 18) 


Noting that the energy of the waves is proportional to J, and is given by 


6, = mnv? = I, (IV. 19) 
we can interpret eqn. ([V.17) as the wave energy transfer equation. The 
second term on the left hand side describes the diffusive energy transfer 
towards the smaller wave numbers due to the scattering of waves on the ions, 
and the second term on the right hand side describes the non-linear damping 
of the waves at the ions. As we have established earlier, in the principal 
region of wave numbers k*Imé = const, so that y, is proportional to k, 
Y = ak say. Thus, from eqn. (IV.17) in steady state conditions, when 


a Aa a el The value of ky can be d 

` z= | 6 ace = $ = 
Ot k FAR? k where Ko const e value of Ky can be deter 
mined from the condition that the spectral function vanishes at k ~ D7}, 
since at larger values of k the growth rate becomes negative. From this we 


; = ] 
obtain kọ ~ D~*. Moreover, the spectrum J, ~ k`?’ In yp must be cut off 


at small k, where Im e changes its sign. 

We note that the non-linear terms in eqn. (1V.18) do not have a “ wave 
number” N, = &,|@,|~*. In other words, according to this equation 
f N,dk = const for y = 0. This result follows immediately from eqn. (IV.14). 
Thus the non-linear interaction of the ion-sound waves constitutes a coherent 
scattering at the ions, during which the “ wave number ” is conserved and the 
energy decreases in proportion to the frequency. 

We obtain the frequency dependence from the relation Adk = I,dw, 


which gives I,, ~ @~* In a Thus the energy of ion oscillations in the 
steady state is concentrated in the region of low frequencies, that is near the 
lower limit of stability. On the other hand, the spectral function of the 
electric field E$ = k?1,, which determines the effect of the oscillations on the 
electron distribution function, reaches a maximum for k ~ D™}. 
Substituting for the growth rate in eqn. (IV.17) the value given by linear 
theory, we obtain the following expression for the spectral function J, for 
u > C: 
pa” TE . 7 1 
~~ — — LILY 
K v, T, 702 4ne?k? kD 
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This is a very large value, and it is obvious that in the presence of oscillations 
of this magnitude the distribution function of the longitudinal electron 
velocities cannot be Maxwellian. In other words, we need to consider the 
feed-back of the oscillations on the averaged distribution function. For 
this purpose it is sufficient to use the quasi-linear approximation. 

In the presence of a sufficiently strong longitudinal magnetic field, such 


_ et 
that Q, oe > Qp, the motion of the electrons across the magnetic field 


can be neglected, and the oscillations lead to the development of a plateau 
on the longitudinal velocity distribution function in the interval from c, to u. 


e? 
Ex, k 
As a result the growth rate is reduced by the factor 2" 2 Zo 7v7 where 
e S 


ky ~ D™* and qt, is the mean collision time between an electron and the 
neutral gas atoms. Taking this effect into account the spectral function 
becomes to order of magnitude 


1 J m,T, 1 1 T? 
I~ ——=-+ In ——= IV. 20 
E 100. N c, m; T; Ooto kD ek? ( ) 


so that even when the formation of the plateau is taken into account the 
oscillation amplitude is still fairly large. 

In the absence of a magnetic field the effect of the oscillations on the 
electron distribution function becomes even more important. In this case 
each separate wave with wave vector k sets up a plateau in the k direction 
Over a small region of velocity space of the form of a thin layer. This layer is 
perpendicular to k and situated at a distance c, from the co-ordinate origin 
(see Fig. 15). Since c, < v,, the aggregate of all oscillations sets up a plateau 


Fic. 15. Isotropisation of distribution function by ion acoustic oscillations 


in the angular dependence of the distribution function in the resonance 
interaction region. Thus the ion acoustic oscillations developing in some 
cone of directions with aperture angle 0), must tend to make the electron 
distribution function isotropic in the external volume of a cone with aperture 
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angle n/2—@05. For u > c, so that in an unbounded plasma practically all 
waves with a positive projection of the phase velocity on the z axis are un- 
stable, the angle 09—>7/2, and the oscillations should lead to a complete 
“ blocking ” of the current, so that the directed velocity of the electrons 
could not exceed c,. In a bounded plasma, let us say in the plasma of the 
positive column of a glow discharge, together with the effects just described, 
we should also take into account the energy transfer of the waves to the walls 
Ol; 
oY 
propagating across the discharge should be damped, and the angle 6, may be 
fairly small even for u > c.. 

Let us note yet another circumstance. Above, we estimated the value 
of the spectral function averaged over the angles. If we wished to refine 
this calculation by determining the angular dependence, we would have to 
solve an integral equation with a degenerate core, and such an equation by 
no means always has a solution. In other words, the angular dependence of 
Im £ may turn out to be different from the angular dependence of the integral 
term in (IV. 16). The integral term then cannot completely compensate the 
growth rate, and there must be directions in k space in which the oscillation 
amplitude will continue to increase. We might expect, therefore, in the case 
of ion sound oscillations in an unbounded plasma, a tendency towards the 
formation of separate wave packets with fairly sharply defined angular 
dependence. 

Moreover, from (IV.16) we see that the interaction of the waves increases 
as the angle 0, between their wave vectors increases. Thus an individual 
mode may strongly interact with its ‘‘ distant neighbours ” but only weakly 
with waves travelling in the same direction as itself. Therefore the non-linear 
interaction not only does not resist, but even strengthens the tendency 
towards a localisation of the waves in a few separate directions. In an 
unbounded plasma this process will obviously develop until four-wave 
processes become important. In a bounded plasma we must add the term 
oe A into the kinetic equation for the short wave oscillations and this term 
will also affect their angular dependence. However, for the oscillations with a 
long wavelength of the order of the dimensions of the inhomogeneity, this 
approach is not suitable and in those cases where the condition for their build- 
up is fulfilled, we should expect a predominant development of one of these 
characteristic oscillations. 


(b) Build-up of Cyclotron Oscillations by the Electron Current 

The ion-wave instability considered above only occurs in a non-isothermal 
plasma where T, > T, For T; ~ T, ion acoustic waves do not occur be- 
cause of the strong damping on the ions, and the growth of oscillations is 
possible only for u 2 v, (77). However, in the presence of a magnetic field, 
which increases the “‘ elasticity ” of the plasma, oscillations can be excited 
at smaller values of the directed velocity u. 


by adding the term U-— to the kinetic wave equation. Accordingly waves 
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Let us first consider the case of a strong magnetic field where the magnetic 
8p 
field pressure is much larger than the plasma pressure, i.e. p = H2 < l. 


Drummond and Rosenbluth (78) have shown that in this case the longitudinal 
electron flux may excite ion cyclotron oscillations. For P < 1 these oscilla- 
tions are purely longitudinal, so that we can write E = —Vọ. Their frequency 
is given by the relation (75): 


p -=s -s -9; 
eH , i i 
where Q; = me B the ion cyclotron frequency, s = k? p? = 7.02 T,(s) 
; i ; 2T, 
is the Bessel function of imaginary argument, and v, = ae 


Equation (IV.21) refers only to waves with small k,, since the ion damping 
loa—-Q; 
omitted in this equation can be neglected only for k, < at Since the 
Te T, 
x —e *I,(s) attains a maximum of ~ 0-2— for s = 1'5, 


i 


i 


„ O- 
quantıty 


Q, TT, T, 


y” 


Fic. 16. Dependence of frequency and growth rate of cyclotron oscillations 
on wave number 


T; Q; 


we have k, S$ — i - Tv Tv and then according to eqn. (IV.21) the instability 
develops only for 
T, TN Jm, 
u >u. T. — 0; = 1s(2) a Ve (IV. 22) 


TN? 
For a hydrogen plasma this condition becomes u 2 ole) . The 
dependence of the frequency w and the growth rate y on the wave number ky 
is shown qualitatively in Fig. 16. This relationship between y and k, applies 
when the instability condition (IV.22) is satisfied by a considerable margin. 


As u decreases the band of unstable wave numbers becomes narrower. 
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Within the framework of the quasi-linear approximation, the problem of 
the sustained oscillations has been studied in the paper quoted above (78). 
These authors evaluated the oscillation amplitude at which a plateau appears 
in the electron distribution function. However, in the presence of a longi- 
tudinal electric field the development of the oscillations does not stop at this 
stage, since the electric field tends to destroy this plateau, and the distribution 
function will always tend to relax towards a displaced Maxwellian distribu- 
tion, and the oscillation amplitude will continue to grow until the non-linear 
processes enter into the picture. 

Oscillations with a frequency slightly larger than the cyclotron frequency 
set up beats at frequencies close to zero and to double the cyclotron frequency. 
The oscillations with very low frequency transfer energy to the ions very 
effectively and it will be these in particular which will lead to the limitation 
of the oscillation amplitude. The problem of sustained cyclotron oscillations 
has been considered by Petviashvili (79) and Karpman (80). 

Since the phase velocity of cyclotron waves parallel to the magnetic field 
w/k, ~ u is considerably larger than the thermal velocity of the ions, the 
non-linear damping of the waves, as in the preceding paragraph, will be 
described by a term containing the function 6(@—’) beneath the integral. 


When we integrate the ô function with respect to k, we obtain a factor of the 


form , which vanishes at the point where w reaches its maximum 


da 
dk, 
as a function of k,. Consequently, at this point the amplitude of oscillation 
in the steady state must vanish, even though according to Fig. 16 the growth 
rate is greatest at this point. 

When we consider the formation of the plateau in the electron distribution 
function, we obtain from the non-linear damping the following estimate of the 
equilibrium spectral function 


an a a 


TT; u? ł dk, 
a IV. 23 
k T ea 


dk, | Idk, =A =F A 


, eH 1. 
where A is a numerical factor of the order unity, 2, =~, z Í$ the 
e e 


electron collision frequency and Z, the spectral function of the electric 
potential. This expression refers to the case where u is much greater than 
the critical value u,. It gives a fluctuation level considerably exceeding the 
thermal fluctuation density which in this region of wavelengths is of the 


T- _ 
order of magnitude of | Ie dk ~ zz (nopi) - 


(c) Magnetic Sound Build-up 


The phase velocity of the oscillations under consideration is of the 
order of magnitude of u, and as we increase 8 = 8zp/H7? so that the phase 
velocity approaches the Alfvén speed, they can no longer be considered 
longitudinal. In other words, in this case curl E # 0, the oscillations 
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become much more complex, and the growth of longer wave magnetosonic 
oscillations becomes possible. 

The dispersion relation for magnetosonic oscillations including Landau 
damping has been determined in (81, 82). Extending these results to include 
the case of a longitudinal current, we obtain for œ < Q;: 


n E | 
or ak irii Eg u—ke,)—* i} (IV. 24) 


vek, * m;c? k? 
This shows that instability occurs = u > C4, but the growth rate is com- 
paratively small, given by y/w ~ p- SE 


Equation (IV.24) can be wee ‘approximately almost up to a ~ Q,, 


corresponding to kp; ~ VE P. When the wave number increases above this 
value, the magnetosonic oscillations change fairly rapidly into purely elec- 
tronic oscillations corresponding to the so-called “‘ atmospheric Whistlers ”’ 
The ions do not participate in these oscillations which cannot therefore be 
excited by a longitudinal current. Also for shorter wavelengths we have the 
cyclotron branch, which we considered earlier, which for P < | connects 
with the ion acoustic or slow magnetosonic wave. 


Thus the magnetosonic oscillations can be excited only for kp; < aif p. 
Their dispersion relation is of the non-decay type (see Chapter II.1(b)). 
It is true that they can excite Alfvén waves by resonant transfer, but since the 
latter are practically not absorbed by particles, this effect does not lead to a 
limitation of the oscillation amplitude. Consequently even in this case the 
amplitude of the steady state oscillations is determined by the non-linear 
Landau damping at the beat frequencies. The mean square of the amplitude 
of the steady state oscillations E can be estimated from considerations 
of the energy balance. We compare the rate of energy transfer from the 
electrons to the ions with the rate of absorption of energy by the ions, and 
obtain approximately 


~ 72 tate Ms 2 4 


where Ej is some mean value of the spectral function, k ~ p; ' JB Byw~ op — 


is the linear growth rate and ĵ the growth rate as reduced by the maio 
of the plateau in the electron distribution function. The right hand side of 
(IV.25) represents the non-linear damping of the waves. To estimate the value 
of B we observe that the amplitude of the low frequency beat oscillations is 
directly related to the non-linear terms in the equation of motion for the ions 


m,k c*E? 
and can be estimated as E, S Pi The energy of these oscillations 
61,2 
is in order of magnitude equal to m;n—, ~ EP = m?n——-E; and its 


= e? H° 
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21,2 2 
. e C 
absorption rate can be estimated as Bc,k——,; TA 


factor p takes into account that the interaction leads to damping only for 

oscillations with frequencies sufficiently close that the phase velocity of the 

beats is of the order of the ion thermal velocity (cf. the factor 7,/T, in eqn. 

(IV.16)). Comparing this result with eqn. (IV.25) and noting that the energy 
2 


where the 


e e . ° e C ° 
of the main oscillations is of the order of magnitude ~ m,nE; m we obtain 
2e? 
Q?H?” 
(IV.25) and setting k ~ p; * q B, we finally obtain an approximate value 
of the intensity of the oscillations of the electric field 

H?v7 u 1 
c* V cave M Qt, 

It is apparent that the kinetic energy of the steady state oscillations is 
broadly speaking Fy Qt, times smaller than the thermal energy, and is 
concentrated in oscillations with frequencies near the cyclotron frequency Q,. 

Both magnetosonic and cyclotron oscillations lead to a plateau of width u 
in the electron distribution function. As a result the electron current decreases 
by a small fraction of order u/v,. This decrease of the current for a given 
electric field can be represented as an additional “ anomalous ” resistivity. 
Consequently the effective conductivity of the plasma ø is defined by the rela- 


tion 
1 1 
Sx (1 7 =) (IV. 27) 


O Oo r 


approximately B ~ u;kB Substituting this value in the relation 


where g, is the conductivity due to binary collisions. 

It can be shown that one-third of the energy dissipated by the additional 
resistivity is transferred to the oscillations and so to the ions. As a result 
additional heating of the ions occurs. In stationary conditions, when there 
are no energy losses from the ions, and the electron temperature is maintained 


at a given level, the anomalous heating of the ions must increase their 
"2 


, u 
temperature to such a value that the energy transferred to the ions a 
e~o 


m, n A 
is equal to the energy — — (ETT), transferred by the ions to the 


nN; Te 


electrons. Writing co = e*nt,/m,, j = enu we have 


1 m; u? 
TET. (1 Fa ot) (IV. 28) 
3m, v 


e 


Thus for a sufficiently large value of u/v, this turbulent heating of the ions 
may lead to a considerable “ breakaway ” of the ion temperature from that of 
the electrons. 
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In addition to this anomalous heating these oscillations must lead to 
an enhanced diffusion of the plasma across the magnetic field, since both 
the electrons and the ions participate in the oscillations. The diffusion is 
determined by the magnitude of the mean drift along the density gradient, 
parallel to the x axis, and the diffusion flow can be represented by 


d fe 


Substituting here the value given by eqn. (IV.26) for Ej and considering that 
the only contribution to the integral comes from the resonance region, we 
obtain for the diffusion coefficient the approximate value 
2 
D 
Diy Go re (IV. 30) 
Ve VeC a JQ. 

where Dg ~ pw; is the Bohm value. 

The cyclotron oscillations lead to a diffusion coefficient of the same order 
of magnitude. 


3. DRIFT INSTABILITY OF A PLASMA 


In this section we shall deal with the drift instability of a plasma which 
gives rise to anomalous diffusion of a rarified plasma in a homogeneous 
magnetic field. The drift instability was described by Tserkovnikov (19) 
who showed that in the presence of a temperature gradient in an inhomo- 
geneous plasma oscillations may be excited with a phase velocity across the 
magnetic field of the order of the drift velocity of the particles. Such oscilla- 
tions are naturally called drift waves. 

In Tserkovnikov’s paper (19) only waves which propagate across the 
magnetic field were considered. Later Rudakov and Sagdeev (20) considered 
the more general case of oblique drift waves which transform into ion acoustic 
waves as the angle between the wave vector and the magnetic field decreases. 
It has been shown in (19) that in the presence of a temperature gradient these 
may be growing waves. 

A considerable influence on the further development of the theory of 
drift instability was the paper by Rosenbluth et al. (21) on the finite Larmor 
radius stabilization of the flute instability of a plasma. This paper stimulated 
a whole series of further investigations of the drift instability of the plasma 
with respect to short wave perturbations with transverse wavelength of the 
order of the mean Larmor radius of the ions. 

These papers were concerned with the investigation of the stability of 
a plasma of such low density that particle collisions could be neglected. 
On the other hand, drift waves are also obtained in the investigation of the 
stability of a weakly ionised plasma where collisions between charged particles 
and neutral gas atoms are dominant. Timofeev (23) (see also (91, 92)) has 
shown that as the neutral gas pressure in a weakly ionised plasma is reduced, 
a peculiar instability appears which arises from a combination of the drift 
motion of the charged particles across the magnetic field and the diffusion 
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along the field. We shall call such an instability a drift-dissipative instability. 
It has been shown by Moiseev and Sagdeev (24) that such an instability also 
occurs in a dense fully-ionised plasma. It has also been shown in (85) and 
(92) that thermal conductivity and viscosity may appear in addition to dif- 
fusion as dissipative processes leading to instability. 

In this section a brief survey will be given of the drift instabilities in an 
inhomogeneous plasma (a more comprehensive survey is due to Mikhailovskii 
(22)). We shall assume that the plasma pressure is much smaller than the 


l A 8 
pressure of the magnetic field, i.e. 8 = = <1. The magnetic field will be 


assumed homogeneous. 


(a) Drift Waves in an Inhomogeneous Plasma 


Consider a low pressure plasma in a homogeneous magnetic field H 
parallel to the z axis, and suppose that the plasma density varies along the 
x axis. For simplicity we shall assume that the plasma temperature is con- 
stant, and that the ions are cold. Then in equilibrium the ions are at rest, 

T, dn eH 


— where 2, = —. 
m,.Q,n dx mE 


e 


and the electrons drift with a velocity v9 = — 


. d 
This is the so-called Larmor drift. We assume that n < 0 so that vo > 0. 


We consider perturbations with transverse wavelength considerably 


ld 
l1 where k = ==" such 


n dx 
perturbations can be represented as plane waves of the form exp (— iwt +ikr). 
We limit the discussion to start with to the case in which œ < k,c,, and we 
assume the electric field curl-free, so that we can write E = —Vg. Thus 
we are discussing first the longitudinal oscillations of the plasma. Provided 
k, is not too small, so that the phase velocity of the wave is considerably 
smaller than the thermal velocity of the electrons, the latter follow a Boltz- 
mann distribution, so that the perturbation of the electron density n, can be 
expressed in terms of Q: 


smaller than the characteristic dimension a = K` 


ne e 
epee a (IV. 31) 
n 

where n is the unperturbed density and T, the electron temperature. The 

expression for the perturbation of the ion density n; can be obtained from 

the continuity equation 


On: 
~ + div(nv, =0 (IV. 32) 


where v, is the macroscopic ion velocity, which can be determined from the 
equation of motion 
Ov; e e 
—=-— --Ve+—-f[v,H IV. 33 
Ot M; 4 ae ] l ) 
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For œ < Q; we then obtain the transverse and longitudinal components 
of the velocity 


Va =i— = [hk] pook 19 (IV. 34) 


H Q; 


and 


v;,, =——@ (IV. 35) 


For œw < Q; the second term in (IV.34) can be neglected and then the velocity 
of the ions is determined entirely by the drift in the electric field. In this 
approximation, using relations (IV.33), (IV.34) and (IV.35), we obtain 


l ET 
Ta & a) st (IV. 36) 
n wo m,w*/ T, 
where 
TKK 
Ox = k, vo = 0 


In a dense plasma, where the Debye radius is negligibly small, all oscilla- 
tions can be considered quasi-neutral, i.e. n; = ne Substituting in this 
relation the values for the perturbations of the electron and ion densities 
obtained above, we arrive at the dispersion equation 


w* —ww, —k2c? = 0 (IV. 37) 


l- l a 
where c, = J r is the velocity of sound. This is in fact the required 
i 


dispersion relation for drift waves. The dependence of œw on k, determined 
by this equation is shown in Fig. 77. Curve (1) refers to a wave propagating 


1 
„KzĆs 
7 


Kz 


-KzCs* 


Fic. 17. Dispersion relation between œ and k, for longitudinal 
oscillations of an inhomogeneous plasma 
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in the direction of the electron drift, for which w/k, > 0, and Curve (2) to the 
wave propagating in the opposite direction, 1.e. with olk, < 0. As we can 
see, the absolute values of the frequencies of these two waves differ slightly 
from one another. In particular, for k,— 0 the frequency of the first wave 
tends towards w,, while the frequency ‘of the second tends towards zero. 
The transverse phase velocity of the first wave for k, = 0 is equal to the 
drift velocity of the electrons. In order of magnitude it is equal to c,ry/a, 


I; e ° . 
where r,= J zo ! the Larmor radius of the ions at the electron 
Séi 


ee 


temperature (for the usual Larmor radius of the ions J fi z We retain the 
i? ti 

notation p;). For k,/k, > rg/a both drift waves go over into ion acoustic 

waves. Therefore it is natural to call the first the accelerated, and the second 

the decelerated ion acoustic wave. In a strong magnetic field ry < a, and in 

this case the transition to a simple sound wave takes place for almost trans- 

verse propagation where k,/k, < 1. 

Returning once more to the derivation of the dispersion equation, it is 
easy to see that the effect of the inhomogeneity on the oscillation frequency 
arises from the transverse drift of the ions. By itself this motion is incompres- 
sible, i.e. div v;, = 0, and in a homogeneous plasma it does not lead to a 
change in density. However, in an inhomogeneous plasma, even an in- 
compressible displacement of the plasma ¢ along the density gradient leads to 


d ; 
a perturbation of the density n; = —¢ = This perturbation of the density 


leads to the change of the dispersion equation. 

These considerations remain valid when the ion temperature is different 
from zero. For if the transverse wavelength is considerably longer than the 
Larmor radius of the ions p;, and the phase velocity parallel to the z axis 
is considerably larger than the thermal velocity v;, we can still use the hydro- 
dynamic equation of motion (IV.33) only including the ion pressure gradient. 


c 
This term leads to the Larmor drift y LbVp:], which for H = const is 


incompressible and therefore supplies no contribution to the change in 
density. Thus, for T; Æ 0 the dispersion relation ([V.37) is unchanged 
provided w/k, > c,. It follows, therefore, that even for 7; = T,, so that 
c, ~ v, and ion acoustic waves do not occur in a homogeneous plasma, a 
drift wave (the accelerated sound wave) may propagate in an inhomogeneous 
plasma with a frequency œ ~ w,. The phase velocity of this wave parallel 
to the magnetic field may considerably exceed the thermal velocity of the 
ions, and consequently it is not subject to strong ion Landau damping. 

As k, decreases the phase velocity of the accelerated wave increases 
and when k, is sufficiently small it may reach either c4 or v,, whichever is 
smaller. We limit our consideration here to the case c4 < ve 1e. B= 


82p/H? > m,/m,;. 
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For = ~ c, the electric field can no longer be considered curl-free. For 


a ~ c, the lines of force of the magnetic field may be said to be no longer 
completely rigid and they become slightly bent. In this case it is possible 
purely formally to introduce a longitudinal potential y, determined by the 


l ô E i ee 
relation E, = so, while since in the case of the slow oscillations under 


consideration the transverse component of the electric field can be assumed 
curl-free, we can retain the transverse potential ọ; E, = —V,q. By using 
these two potentials we take account of the bending of the lines of force, but 
continue to neglect any changes in field strength due to compression of the 
field. For Y = of course even the curvature of the field disappears. 

Since according to our assumption the thermal velocity of the electrons 
considerably exceeds the phase velocity of the wave w/k, ~ c4, the electrons 
reach equilibrium along the lines of force leading to the following relation 

H? dn |. 
H dx ik eny (IV. 38) 
where H, is the x-component of the perturbation in the magnetic field given 


c 
by H; =— 


T,ik,n,.+ T. 


c 
z [kE]; =i z okle- Thus we obtain for the perturbation 


of the electron density 


EEE a fee as 
en 4 Sy (1 P) (IV. 39) 


We need to relate o and w. For this purpose we use the z component of 
° 2 


ð c 
the equation 7 ae curl curl E, where jis the current density. For these 


oscillations we obtain 
2 


c 2 
Ta k,ki(o—-w) (IV. 40) 


To determine j,, we use div j = 0, from which we obtain 


J: = 


j, = ik; + divj, = ik; ‘edivn(v,, —Vv,,) 
Thus the longitudinal current arises from the small difference in the trans- 
verse velocities of the electrons and ions, so that in the expression (IV.34) for 
the ion velocity we must consider the second (inertial) term. In addition there 
is a small difference between the ion velocity and the electric drift velocity 
arising from the effect of the finite Larmor radius (21). In a strong magnetic 
field, the mean Larmor radius of the ions, although small, is nevertheless not 
negligible and the electric drift of the ions is not determined by the electric 
field at the centre of the Larmor circle, but by some value averaged over the 
Larmor circle. It is easily verified that when we average the electric field of 
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a plane wave over a circle of radius v,/Q; and subsequently over the Max- 
wellian distribution in v,, the effective electric field is less than the value 


kif 
at the centre of the circle by a factor e~* J,(s), where s = =r and J, is the 
mì; 


Bessel function of imaginary argument. For small s this factor is approxi- 
mately equal to 1— s, and in this case a more accurate expression for the ion 


velocity may be written down in the aaa form 
cC @ 
s T e+ Fa ki (IV. 41) 


Substituting this expression for the ion velocity in div j, and noting that 


va = [hk] — [nk] = 


v1. = = [hk]o we have 


nm;c? k? 
z= =o p Ote (IV. 42) 
whence, using (IV.40), we obtain the following relation between @g and w 
basg ll (IV. 43) 
ke, 


Using this relation and (IV.31), we obtain the perturbation of the electron 
density: 


2. 9 
He ¢ a =) - (IV. 44) 


Finally, comparing this equation with the perturbation of the ion density 
(IV.36), neglecting the second term in brackets, we obtain the dispersion 
relation 


(w — o) (0? +@,@—c*k?) = 0 (IV. 45) 
This splits up into two: @ = w,,, which describes the drift wave which we 
already know, and 
w +o0,0—cik? = 0 (IV. 46) 
which shows that for k, > 0 the Alfvén waves in an inhomogeneous plasma 
also go over into drift waves. To differentiate these waves from those con- 
sidered earlier, we shall call one the accelerated Alfvén wave (w/k, > c4), 
the other the decelerated Alfvén wave (w/k, < c4). 
Retaining small terms of the order k?p? previously neglected, we obtain, 


for p > ae instead of an actual intersection of the branches, the picture 


shown in Fig. 18. There, as previously, waves with positive frequency 
(branches 1 and 3) propagate in the direction of the electron drift, and waves 
with negative frequency (branches 2 and 4) in the direction of the ion drift. 
In an inhomogeneous plasma, the usual sound and Alfvén waves are replaced 
by four different waves. 
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Fic. 18. Dependence of frequency of drift waves on wave number 


In the following sections we investigate the stability of an inhomogeneous 
plasma to drift waves. In a collisionless plasma the instability arises from the 
interaction between the drift waves and resonant particles represented by 
Landau damping, which must be discussed on the basis of kinetic considera- 
tions. This calculation is carried out in Sections (b) and (c) where we consider 
separately long wave (k,p; < 1) and short wave (kip; 2 1) perturbations. 
In Section (d) we consider the stabilization of the flute instability due to 
finite Larmor radius. In the following Section (e) the instability of a plasma 
with cold ions (T; = 0) is considered in a weak magnetic field, when the 
Larmor radius of the ions calculated from the electron temperature exceeds 
the characteristic dimension of the inhomogeneity. In Section (f) the growth 
of cyclotron oscillations in an inhomogeneous plasma is considered. The fol- 
lowing three sections are then devoted to the consideration of the drift- 
dissipative instability, related to collisions between particles, and in the last 
section non-linear drift oscillations will be discussed. 


(b) Drift Instability for kp; <1 


To describe the oscillations of a collisionless plasma with a wavelength 
much greater than the mean Larmor radius of the ions, the drift kinetic 
equation can be used (the same approximation has been used in ref. (20) 
which we follow in the present section). We shall limit the discussion to the 
case of longitudinal oscillations, and for these the linearised drift kinetic 
equation for particles with mass m and charge e has the following form 

e Of c Of 
—wot+k ‘'—-—k,o-- —--k,o— =0 IV. 47 
(—O+k wf — koa o (IV. 47) 
where f is the unperturbed distribution function. With the aid of (IV.47) we 
obtain the density perturbation 
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n’ =—e MEE ky =) IV. 48 
j m ôv, mQ x) @w—k, v, +iv aa 


where we have added the small positive quantity v -> 0 to pass the pole 
correctly. In the special case that the distribution function fis Maxwellian 
and the temperature constant, the density perturbation n’ can be expressed 
in terms of the function 
+o z 
— x? dx 


Y — =? =<? t2 4 anz? 
(z) = hi eas e fe TEn (IV. 49) 
0 


pe e | ee eo btain f 
SO ra Tay TK so that we obtain from (IV.48) 
n’ O+ Wy w ep 
—= Y -1 IV. 50 
n { kz vr a E T l ' 


where vr is the thermal velocity. Using the quasi-neutrality condition we 
obtain the following dispersion equation 


O + Wy w wW — Wy w 
Y Y =2 IV. 5] 
k, vi G -) k, Ve 6 = i ) 
where 
2T 2T 
a or =f ae eae Oe Wy = Ky vo 
M; m 
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Let us consider the drift waves with phase velocity along the z-axis in the 


3 
range vu; < kL <v,. In this range ion Landau damping can be neglected. 


For generality we allow the temperature to vary in the x-direction and the 
electron function to be shifted relative to the ion function by a velocity u so 
that there is a longitudinal electron current. In this case the dispersion 
relation assumes the form 


Ox ke (1 Wy DE in 
o 


+ — — k,u—@yz+4@,4)=90 (IV. 52) 


d\nT. . 
where 4 = T For 7; = 0 this equation differs from (1V.37) by a small 


imaginary term which arises from the interaction of the wave with the 
resonant electrons, the longitudinal velocity of which coincides with the phase 
velocity of the wave. For w/k, > c, the third term in (IV.52) can be neglected 
so that we have approximately 


— u lo | 
O = Wy, y= Jy (= - eee n) (IV. 53) 
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For — — <x | and —* * < v, the growth rate is much less than the frequency. 


7 k, 
Moreover, in the absence ofa Kona current the P occurs only 
f dinT, 
orn = Jinn 
quently long wave parte for k,p; < 1 are unstable ai in the pre- 
sence of a longitudinal current, although for sufficiently large values of y, to 


be precise for n > 2, perturbations with k, ~ = become unstable (20). 


e 

The instability of the drift waves for u # 0 is an extension of the ion wave 
instability (see (IV.2)) to the case T; ~ T,, when there is no normal ion 
sound. On the other hand, the instability for u # 0 can also be regarded as 
similar to the current-convective instability (see I.2). In effect, as in the latter 
case, the perturbation increases due to the drift of the particles in the electric 
field of an oblique wave. The transverse component of the field, which leads 
to the drift, arises simply as a consequence of the perturbation of the longi- 
tudinal electric field due to the interaction of the resonance electrons with the 
wave. The only difference between the instability considered here and the 
current-convective instability is that in this case the effects of collisions are 
produced by Landau damping. 


(c) Drift Instability for kıp; Z | 


It is evident from (IV.52) that for n = 0 and u = 0 the growth rate 
vanishes only because of what would appear to be an accidental circumstance, 
namely that the frequency of the oscillations exactly coincides with w,. Any 
effect which shifts the oscillation frequency from the value w, leads to growth 
or damping of the drift waves. Ina plasma with cold ions (7; = 0) it may be 
necessary not to use the drift approximation for the ions but to include the 
inertia term in the equation of motion, i.e. the second term in (IV.34). When 
the corresponding change in the expression for the ion density is made, it can 
be shown that eqn. (IV.52) takes the following form 


KiTe |; _@« Keck | ivn 
mQ? w o? ue 


The transverse inertia term becomes important for k,ry ~ 1, where ry 
is the Larmor radius of the ions at the electron temperature. There is then an 
instability even in a currentless plasma for 7 = 0. 

When the ions are hot a similar effect occurs for k,p,; ~~ 1. For 
Te = T; = T and 4 = 0 the dispersion equation for longitudinal oscillations 
(w < k,c,4) can be obtained by solving the kinetic equations for both the 
electrons and the ions and has the form 


ens w+w u) w — w, w 
Io(s ) PP 2Y r(e) ea Pa = (IV. 55) 


z“i 


——(w—k,u-—w,y+4o0yn) =0 (IV. 54) 
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2T 2T k? T 
i ne -JZ = kip? = +, i 
where i Ve S S 1 Pj mQ? lọ is the Bessel 


function of imaginary argument, and Y is the function defined by eqn. 
(IV.49). Equation (1V.55) differs from (IV.51) only in that in the first term an 
additional factor e~ *Jo(s), appears, which represents the decrease of the ion 
drift velocity due to the averaging of the wave field over the Larmor orbit 
(see Section (d)). 

This equation was investigated in ref. (83); the results were that only 
the accelerated sound waves, which propagate in the direction of the electron 


drift (Curve 1 on Fig. 17), are unstable. For v; < = < v. in (IV.55) we may 


k 
kv; - 
put n = = — y(— ~ —i\/n and we then obtain 
k, v; W k, Ve 
Ps ae Wy B1— B) 
as) = e IV. 56 
osoz y VNT n OBY 


where p, = e~h (s), and for k,v, > w, the growth rate is small compared 
to the frequency. However, if we decrease k, until the phase velocity of the 
wave parallel to the z axis approximates to v,, the growth rate y becomes of 


-4 [Mi 
Kzas 24m, 


Fic. 19, Dependence of growth rate y on longitudinal wave number kz 
for the drift instability 


the order w. The results of the numerical calculation of the growth rate as a 
function of k, near the maximum are shown in Fig. 19 for three values of s 
and for k, = 0. 

For small s the unstable oscillations can be divided into two branches. 
The frequency of one of the branches is given by the relation (IV.56) and the 
corresponding growth rate decreases as s? for s > 0. The frequency of the 
second branch can be obtained by using the asymptotic expansion 
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1 1 
Y(Z) = F -+ 553 for z > 1. Equation (IV.55) then becomes the cubic 


2 
==, (IV. 57) 


kz Wy PRT ; 
im For a > 0:4 complex roots appear in this equation 
L Z 


and the corresponding growth rate is of the order of w,. On Fig. 19 section 1 
corresponds to this instability. 
Figure 20 shows the boundary of instability for k, = Oandm,/m, = 1840. 


where œ? = 0,;Q, 


Tn (1+s) 


Fic. 20. Region of drift instability in a low pressure plasma 
(unstable below the solid curve). For dashed curve see text 


For small s it is defined by the condition œ < œ, andis given approximately by 
ak, = s, whilst for s 2 1 the limiting value of k,a is determined by the ion 
damping and is almost independent of s. The dashed curve on this figure 
is the locus of maximum growth rate as a function of k,a. 

The above relations refer to the case œ < k,c4, i.e. they are correct for 
the whole region of k, only if v, & cy, and consequently B < m,/m;. In the 
opposite case $ z m,/m;, we must consider the curvature of the lines of 
force when k, is small. Mikhailovskii and Rudakov (84) have shown that in 
this case the instability develops only on the accelerated sound branch 
(Curve 1 on Fig. 18). Forn = Oands 2 1 the frequency of these oscillations 
is given approximately by 


pes Wx B,C? 


ORT ee (IV. 58) 
and the growth rate y by 
Vro? 
yx (IV. 59) 


|kz|Vebs 
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where € = 2k aj b. The growth rate y as a function of k, reaches a maximum 


wo [me : 
~ BN mB for č x 1. As k,p; increases, the frequency of the oscillations 
S 


. . . mi 
remains approximately constant, the growth rate increases and for s ~ — B 
m 


e 


it becomes of the order of the frequency. 

These results refer to the case of constant temperature. However, when 
we include a temperature gradient the instability still exists at practically 
any given ratio between the density and temperature gradients (85, 86). 

As p increases, the ion thermal velocity approaches the Alfvén speed 
and the ion Landau damping becomes more and more important. At suf- 
ficiently large £ values, as has been shown in ref. (93), this effect completely 
stabilises the drift wave. Figure 21 taken from ref. (93) shows the critical £ 


O-| 


0:05 


Kı Pi 


Fic. 21. Region of drift instability as a function of $ 
(unstable below the curve) 


at which the oscillations become stable as a function of k, p; (stability occurs 
above the curve). For $ > 0:13 the drift instability is completely suppressed 
by the ion damping. 

In addition the drift instability must be absent in systems which are not 
very long. For, as is apparent from Fig. 20, the drift instability develops only 
for perturbations with small k,; more precisely, k,a must be of the order of 
107+. In devices of limited length the wave number k, cannot be small and 
the plasma must be stable. 

A further effect contributing to the stabilisation of the drift instability 
is the shear of the lines of force, refs. (88, 89). We shall characterize the 
magnitude of the shear by the parameter 0 = a/L, where a is the transverse 
dimension of a plasma column and Lis a length along the column such that 
the rotational transform angle referred to this length (94) varies across the 
column by a value of order unity. Broadly speaking 0 is the angle between 

P.T. 4 
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two lines of force separated from each other by a distance a. The projection 
of the wave vector on to the line of force varies across the column by about 
k0. If we assume k, ~ p; ', then even for the very small value 0 ~ p;/a we 
have k0 ~ a™~!, so that at some point in the column we must have a region 
where k, ~ a™ +, i.e. a region of strong absorption at the ions. Thus, for 
L ~ a?/p; we expect stabilisation due to shear. For some perturbations the 
value of 0 necessary for stability may be slightly larger than p;/a because of 
the reflection of the wave at “‘ potential barriers ”, i.e. turning points where 
for the frequency œw, k, becomes imaginary (see refs. 88, 89). But even in this 
case the values of @ required for stability are very small. 


(d) The Stabilisation of the Convective (Flute) Instability 


As we showed above, the inhomogeneity of the plasma leads to a consider- 
able change in the dispersion relation w(k) in the region of small frequencies 
of the order of w*. In other words, for oscillations with phase velocity of the 
order of the drift velocity the properties of the plasma are appreciably dif- 
ferent from those of a conducting fluid described by the magnetohydro- 
dynamic equations. Rosenbluth et al. (21) have shown that the drift effects 
may, in particular, considerably influence the magnetohydrodynamic flute 
instability of a plasma and under certain conditions complete stabilisation 
of this instability appears possible. 

Qualitatively the possibility of this stabilisation may be seen directly from 
Fig. 18. From the magneto-hydrodynamic point of view, the frequency of 
Alfvén oscillations is given by œ = k,e, and for k, = 0 it becomes zero. 
Perturbations with k, = 0 remain constant along the lines of force of the 
magnetic field and it is these perturbations which are referred to as flute 
perturbations. The fact that the frequency of the oscillations vanishes means 
that for these perturbations the plasma shows no “elasticity ” whatever, 
and they can therefore be destabilised by any curvature of the lines of force 
or any gravitational force acting in the direction of decreasing density, no 
matter how small. When we consider the drift effects, however, we can see 
from Fig. 18 that the relationship between the oscillation frequency and k, 
becomes more complex and the transverse motion cannot be considered 
completely inelastic. 

Let us consider this effect in some detail. Suppose that a cold plasma 
(T; = T, = 0)1s situated in a strong magnetic field acting parallel to the z 
axis, and subject to the effect of a gravitational force with acceleration g. 
We shall assume that g is parallel to the x axis and in the direction of de- 
creasing density. In the equilibrium state, in a co-ordinate system where the 
mean electric field vanishes, the electrons are at rest and the ions drift parallel 
to the y axis with a velocity vio = g/Q;. In the perturbed plasma the velocity 


of the electrons is determined by the electric drift [hvo], and the perturba- 


tion of the ion velocity v; can be determined using the equation of motion 
which for perturbations of the form exp (— iwt + ikr) takes the following form 
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e 
—iw'm,V; = —iekọ + zLvH] (IV. 60) 
where w = w+k,g/Q;. For w’ < Q; we obtain approximately 


= ~ [hk]o ie (IV. 61) 


? k 
H Q, 1P? 
and substituting this expression into the continuity equation for the ions, we 
obtain 


o ck, dn Gi cw’ k? 
w'n,;—-i— — IH) 
H dx” HQ, 


and the perturbation of the ion density is given by 


o=0 (IV. 62) 


k 
nwe oa (IV. 63) 
w 


For the electron density perturbation we obtain instead of (IV.63) the much 
simpler expression 


ckyK 
y 
n/n = IV. 64 
because we can neglect the electron inertia term. 
Substituting the values for n; and n, so obtained into the equation k*9 = 
4ne (n;—n,) and assuming for simplicity k, = 0, we obtain 


Qk 1 Í c2 

7, fren =) 1 2 (IV. 65) 
The terms in brackets are due to the difference between the drift velocities 
of the electrons and the ions arising from the gravitational field, the unit 
represents the inertia term in (IV.62) and the right hand side represents 
the effect of possible departures from quasi-neutrality. In a dense plasma, 

where c, < c this right hand side can be neglected. 
There is a large factor of approximately Q; in front of the brackets in 
(IV.65) and therefore the oscillation frequency is considerably larger than 


k . 
k,g/Q,;. Expanding (o+ wy to first order in this quantity, we obtain 


from (I1V.65) the much simpler dispersion equation w*+gx = O and thus 


oO = etd, gk so that the plasma is convectively unstable. 

According to (IV.65), the instability originates from the small difference 
between the drift velocities of the electrons and the ions. Naturally any 
additional factors which may influence this velocity difference must also 
considerably influence the stability of the plasma. One of these factors is the 
effect of the finite Larmor radius of the ions, and including this, using (IV.50) 
for small s, we obtain, instead of ([V.63) 

n; 


- (Sa-9-s) 2 (IV. 66) 


n 
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This shows that the effect of the finite Larmor radius of the ions can be 
included in the dispersion eqn. (IV.65) by simply introducing the factor 1—s 
in the first term in the round brackets. For c, < c the equation can then be 
represented in the form 


w?+o0,0+ 9K =0 (IV. 67) 
whence we obtain 
2 
Ox. Ox 
mE oe en , 
w 5 J 4 Kg (IV. 68) 


2 
w 
This shows that for gk < T the plasma is convectively stable. If we introduce 


the effective radius of curvature instead of g according to the equivalence 


T dg os 
g= 7» the stability condition takes the form 


4 
k?p? > = (IV. 69) 


mQ? anda =K 7 
practice, even perturbations the transverse length of which considerably 
exceeds p; will be stable. 

For s z 1, Mikhailovskii (95) has shown that there is also a slight change 
in the first term in the dispersion equation for the flute perturbations, and the 


dispersion equation valid for all s has the following form 
Q; / e "Iols 1 l—e SI,(s) c? 
( o(s) “ie K) c —— 

k \o+kyg/Q; œ S c 


The result of the investigation of this equation carried out in (90) is 
shown on Fig. 22. As this figure shows, in a dense plasma (c4 < c) all 


where p? = Thus for a < R which is often the case in 
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Fic. 22. Region of flute instability (hatched) 
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perturbations with sufficiently short wavelengths are stabilised, and the 
Stability condition is that given by the relation (I[V.69). As the density 
decreases the plasma becomes unstable, but for a very small density the stab- 
ility is again restored, since ultimately a transition will take place to separate 
non-interacting particles which can no longer be described as a plasma. In 
the intermediate region there is a gap where the plasma is unstable. This gap 
extends into the high density region, but only for extremely short wavelength 
perturbations, which cannot make any appreciable contribution to the 
transverse diffusion. 

We considered earlier only short wave perturbations for which the 
quasi-classical approximation is valid. Qualitatively the conclusion that the 
flute instability may be stabilised also refers to long wave perturbations 
except the so-called m = 1 mode, corresponding to the transverse displace- 
ment of the plasma as a whole. With this type of displacement the perturbed 
electric field in the filament is almost constant over the cross-section so that 
the drift velocities of the electrons and ions coincide closely and the stabilisa- 
tion due to the finite Larmor radius is absent. 

An instability of the m = 1 mode is more conveniently described by the 
momentum equation. For c, < c the momentum of the electromagnetic 
field and the electrostatic forces can be neglected, and the equations for the 
momentum for the displaced plasma filament can be written in the following 
form: 

dP 
7 F (IV. 71) 


where P = m,N oe is the momentum of the filament, N the total number of 


ions per unit length, r, the radius vector of the mass centre and F the gravita- 
tional force. Suppose that g increases linearly with distance from the sym- 
metry axis: g = br. Then the force F = m,Nbr, and we obtain from ({V.71) 
rf. = br or w? = —b, and consequently the plasma filament is unstable 
relative to this displacement. In an actual case of plasma in a trap with 
inhomogeneous magnetic field, the force F is determined by the integral 
. (L;+T,)n l , 
over the volume of the quantıty -n where n is the electron density 
and R the mean radius of curvature of the lines of force. When the mean 
curvature of the lines of force is proportional to the distance from the 


1 r 
symmetry axis, 1.€. RIR? where a and Rọ are constants, the force F 
ano 


Po 
can be represented in the form F = 7R TT where Po = fpdr x const, 
0 


c 


and ry = P,~ ‘frp dr is the radius vector of the “ pressure centre ” of the 
plasma. If the temperature of the plasma is constant over the cross-section, 
ry is equal to r, We then again have ř. = br., so that the plasma 1s 


unstable. 
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(e) lon-sound Instability of an Inhomogeneous Plasma 


We have considered so far only oscillations with frequencies much less 
than both the electron and the ion cyclotron frequencies, so that both the 
electrons and ions were ‘“‘ magnetised’’. Let us now consider the other 
limiting case where the magnetic field is so weak that its effect on the ions can 
be neglected (90). In other words, we shall consider oscillations with a fre- 
quency considerably larger than Q,, but we shall still assume that the frequency 
is less than Q,, so that for the electrons we may use the drift approximation 
(we assume kp, < 1). For simplicity the ion temperature 1s assumed to be 
zero. We Shall limit the discussion to the case of a low pressure plasma in a 
homogeneous magnetic field and to “ potential” oscillations for which 


E = —Vqo. With these assumptions the equation of motion for the ions, 
which we shall suppose at rest in the equilibrium state, takes the form 
wm,V; = eko (IV. 72) 
from which, together with the continuity equation 
on;—kv;n = 0 (IV. 73) 
we obtain the perturbation of the ion density 
k? 
nj/n = S (IV. 74) 
w* m; 


The density of the electrons is determined by relation (1V.50). Comparing 
n; and n, we obtain the dispersion relation for these oscillations 


eT, ao 
oe ety \ at (IV. 75) 


mœ?  k,v, kv. 


For oscillations with phase velocity w/k, < v, this assumes the simpler form 


22 
kcg pW Wy 
7a iN 


=1 (IV. 76) 


ze 


From this we obtain the frequency of the oscillations œ = kc, and the growth 


- R: ( ) l 
} l m, O w ( v. I /) 


From this relation we see that for œ < w, the plasma becomes unstable. 
We have already established that the same condition applies to the growth of 
drift waves. In the case considered here, the condition means that the drift 
velocity of the electrons must be higher than the velocity of sound. It can 
also be put in the form 


rk > 1 (IV. 78) 


Te l 
where ry = J 0 is the Larmor radius of the ions calculated from the 


l E> 


electron temperature. 
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According to (IV.76), the growth rate increases as k, decreases. This 
increase continues as the longitudinal phase velocity approaches the thermal 


m, 
velocity of the electrons, i.e. almost to k,/k ~ J —, when the growth 
M; 


rate attains its maximum value y ~ œ ~ kc, For a further decrease of k, 
both the frequency of the oscillations and the growth rate decrease. 


(f) Cyclotron Instability of an Inhomogeneous Plasma 


When the magnetic field increases condition (IV.78) is infringed, and 
the 10n-sound instability goes over into the drift instability considered earlier, 
which develops for oscillations with a frequency considerably smaller than Q,. 
Mikhailovskii and Timofeev (96), (97) have shown also that oscillations with 
frequencies near multiples of the cyclotron frequency, œ = nQ,, may also be 
unstable in an inhomogeneous plasma. From this point of view, the drift 
instability can be regarded as a cyclotron instability with n = 0, and the ion 
sound instability as an instability at very high harmonics of the cyclotron 
frequency. 

For oscillations with a frequency close to nQ,, when the contribution of 
the other harmonics to the density perturbation can be neglected, the disper- 
sion equation for T; = T, takes the form 


Ota, . fo—nQ\ _ W — Wx w ) 
Y | —— Je ‘I Y =2 IV.9 
k, v; ( k, v; Je n(s) + kv C v ( ) 


where /, is the Bessel function of imaginary argument of order n, and the 
other symbols are the same as those introduced earlier (see (IV.55)). For 
kp, 2 1, we must introduce the additional factor e~**p21,(k*p2) in the 


~ 


second term on the left hand side of eqn. (IV.79). In addition, for very 
large values of the wave number k the quasi-neutrality condition may be 
infringed, and it is necessary to add the term k? D? to the right hand side of 
eqn. (IV.79). 

Equation (IV.79) can be considerably simplified for small kz, where 


E Q; k, i act . ° 
Y sa ) X e The condition for the growth of such oscillations 
k, U; a—nQ, 


can easily be shown to have the form obtained previously, w < @,, 1.€. 
kkp? >n (IV. 80) 


when waves propagating in both the electron drift and the ion drift directions 
are unstable. 

It is interesting to note that cyclotron oscillations can also be built up 
for purely transverse propagation (k, = 0) if the following condition is 
fulfilled 


2\ 4 
Kp, > 2n - $ 4 (IV. 81) 


where n is the number of the harmonic considered. In the case of a rarefied 
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plasma, the quantity = underneath the root sign in (IV.81) can be neglected 
and for n = 1 this condition takes the following form 

Kp; = 2c,/c (IV. 82) 
The growth rate of the cyclotron instability may reach values of the order 


4 
(=) Q,. If the condition (IV.82) is not fulfilled, an exponentially small 


M; 
‘residual ” instability may develop (98). 


£ 


(g) Drift Dissipative Instability 


In (b) and (c) we considered the growth of drift waves in a collisionless 
plasma. It has also been found that the drift instability may develop in a 
dense plasma in which an important part is played by collisions between the 
particles (23), (24), (90), (92), and the collisions are in fact responsible for the 
instability. To distinguish between this instability and the collisionless 
instability, we shall call the former the drift-dissipative instability. 

We assume that the electron mean free path A, is considerably smaller 
than the longitudinal wavelength, so that the diffusion approximation can be 
used for the longitudinal motion. For simplicity we shall limit the discussion 
to the case of a plane layer of plasma, although in fact we have in mind a 
cylindrical plasma column in a tube of radius a ~ k™t. The transverse 
length of the perturbations will be assumed to be considerably smaller than a. 

Let us start with the case of a weak magnetic field, assuming that the 
ions are cold and the Larmor radius calculated from the electron temperature 
ry > a. We assume furthermore that in equilibrium there 1s no mean electric 
field and the ions are at rest. These conditions may occur in practice in a 
positive column at a very low neutral gas pressure. If the mean free path 
length of the ions A, is of the same order as 4., and is larger than a, then the 
ion collisions can be neglected and for oscillations with frequency œw > Q; 


k 
we have from (IV.74) n,/n = ao The perturbation of the electron density 
i 


is obtained from the continuity equation 


ck, dn 
on, + —— o—k,v,,n =0 IV. 83 
H dx? ( ) 
where for v, we substitute the electric drift velocity, assuming that 9,7, > 1, 
so that the transverse diffusion and mobility of the electrons can be neglected. 
The longitudinal velocity of the electrons in diffusion conditions is determined 
by the relation 


0p On, et, /0@ T, ôn 
eR PE a ane Sec ear Quer Mean co 
re = Deas e" Gz mM, p en > ee 


where b, is the mobility, D, the diffusion coefficient and t, the mean collision 
time between electrons and neutrals or ions (for collisions with ions it is 
necessary to substitute the relative velocity for vz, in (IV.84), but since the 
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longitudinal motion of the ions in these oscillations can be neglected, this 
effect is not important). 


Substituting the expression for v,, into the continuity equation for the 
electrons, (IV.83), we obtain 


ne Q oTe eT, (IV. 85) 
n D,k -io m,* l 


For w < D,kz this expression can be simplified to 


n ik,k \ ep 
< = | 1- -— }— IV. 86 
Fao mn 
Comparing n; and n, we arrive at the dispersion equation 
ikk \~? 
* = ket (1 - 2 ) IV. 87 
w Cs Q Ttek? ( ) 


The frequency of the oscillations is complex and one of the roots has a posi- 
tive imaginary part, so that the plasma is unstable. Since œ ~ kc, this in- 
stability may also be classified as ion-sound instability. It is evident from 
(IV.87) that for Q.t, > 1 perturbations will be built up in the first place 
which are strongly elongated parallel to the magnetic field, i.e. k, < ky. 

The growth rate y attains a maximum y ~ œ ~ ck for k} = k,x/Q,1,. 
Then k2D,w~1 = kry > 1, so that the assumption œ < k?D, is fulfilled 
throughout except for the region of very small k, which is not of great 
interest. 

Forr,yk S 1 the longitudinal diffusion of the electrons must be considered 
in (IV.85). It can be shown (90) that in this case the plasma is stable to 
acoustic oscillations with frequency œ > Q; 

A plasma with these parameters can be set up, for instance, in the posi- 
tive column of a glow discharge. In stationary conditions such a plasma will 
diffuse to the wall, and if the diffusion speed is of the order of c,, these oscilla- 
tions cannot develop before reaching the wall, since their group velocity is 
also of the order c,. Since the transverse diffusion coefficient D, 1s of the 
order of 2v. (Q,1,)~*, the growth condition D,« < c, can be represented 
approximately in the following form 

2 
ECR (IV. 88) 
aN m; 


l 


T, , 
where A, is the electron mean free path, r= ee o? a=xKxK!. Thus this 
j 22; 


instability occurs, in the absence of a longitudinal current, only in a sufficiently 
rarefied plasma where the electron mean free path is not too small. 
For A; S a we must also consider collisions between ions and neutral 


bi 
gas atoms. This effect (85) leads to the instability condition 22;t; > J p 


4* 
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where b, is the electron and b; the ion mobility. When the magnetic field 
decreases condition (IV.88) is replaced by the inequality pẹ < a, where pe 
is the mean Larmor radius of the electrons, since for pẹ > a the effect of 
the magnetic field on the electrons can be neglected, and there will be no 
instability. 

Let us now consider waves with frequency œ < Q;. The perturbation 
of the electron density is again given by (IV.85), and including the effect of 
friction between ions and neutral gas, the perturbation of the ion density is 


given by 
y 2 : 
Ta i a (1 cee je (IV. 89) 


wt,/ M; 


for Q;t; > 1, where 1, is the mean collision time of ions with neutral gas atoms 
(see (IV.63)). Comparing n, and n; we obtain the dispersion equation which 
can be conveniently written 


1 D, k? 

w + iw (oD. + =) — i0,0% — -=0 (IV. 90) 
where 
cT, dn k? 
ERR p 726.40) 
ce ”eHn dx = kz Meet 
. 2 kiT, 22 

It is readily seen that D,kz/w, = O27 kirą. Therefore in a strong 


magnetic field, when xr, <1, for perturbations with kr, < 1 we may neglect 
D.k2 compared with w,. If in addition the collision frequency v; = t, * is 


sufficiently small, then (IV.90) can be simplified to 
w*+iww,—iw, o, =0 (IV.91) 


This equation is also valid for a fully-ionised plasma (24), if by te` * 


we understand the mean collision frequency of electrons with ions. It follows 
from eqn. (IV.91) that for w, > @, 
2 
Ww 
w; = Obi“, WM, =- iw, (IV. 92) 


and for œ, < Ox 


O,2= +Jio,o, (IV. 93) 


In either of these limiting cases one of the roots has a positive imaginary part, 
giving an instability. For œ, ~ @, the growth rate attains its maximum 
value y ~ MW ~ Wy. 

In a weakly ionised plasma the collisions between ions and neutrals give 
rise to an additional damping, and for ke,t; < 1 all perturbations are damped 
(90). 

In the other limiting case kry > 1 drift waves with frequency œ <Q, 
may grow as well as ion-sound waves, but their growth rate, which is by 
definition smaller than Q,, is therefore smaller than the growth rate of the 
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ion-sound oscillations. Drift waves can grow only for 2,1; > 1, when the 


[hVo] 


drift velocity : H 


The entire region of the drift-dissipative instability in a weakly ionised 
plasma is represented schematically in Fig. 23. The abscissa is y = rp' K”! 


is larger than b;VYọ. 
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Fic. 23. Region of ion acoustic (vertical hatching) and drift (horizontal hatching) 
instabilities in an inhomogeneous weakly-ionised plasma. For diagonal hatching see 
text. (1) results of Geller (131), (2) results of Golant (111) 


which is proportional to the magnetic field, and the ordinate the quantity 


Te T 
č = tik which is proportional to l/ap, where p is the neutral gas 


pressure and a the tube radius. For y less than unity but greater than 


m i 
nigh ae where the Larmor radius of the electrons is of the order a, 
li 


an ion-sound instability will develop. The corresponding region is marked in 
Fig. 23 by the vertical hatching. The region of the drift instability proper 
(œ <Q,) is marked by horizontal hatching. For yn > 1 it has a lower limit 
č = |, but fory < 1 the instability develops only for Q;t; > 1. 

Thus, for sufficiently low neutral gas pressure, a weakly ionised plasma 
in a homogeneous magnetic field is unstable even in the absence of a longi- 
tudinal current. In the presence of a longitudinal current, the instability 
region expands towards smaller ¢ values, since even for 4; < a the current- 
convective (spiral) instability previously discussed may develop in a plasma 
with a longitudinal current (see I.2.b). Moreover, the longitudinal current 
may also have an effect on the drift-dissipative instability. In the presence of a 
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longitudinal current the expression for the perturbation of the electron 
density changes to 
n, k2+ik,«(Q,1.)' et 
I Balai a y (IV. 94) 
n D,ke-—iwt+ik,u me 
where u is the longitudinal (current) velocity of the electrons. 

Accordingly, in expression (IV.87) for the square of the frequency of 
the ion-sound oscillations, we obtain an extra factor 1 +iu/D,k?, and in this 
case the maximum value of the growth rate again occurs for kz ~ ky K/Q.T.. 
For a sufficiently large Q,t, the maximum growth rate may increase by a 


u Qr 


4 
factor (STE) compared with the maximum value y ~ kc, obtained 


previously. 

The longitudinal current also strengthens the intrinsic drift instability. 
Using eqn. (1V.94) the dispersion equation becomes more complex than 
eqn. (IV.91), taking the form 


k 
w + iw(@,+ik,u)— iw, (osi a ) =0 (IV. 95) 


kiti 
2 
where, just as before w, = pz Sete Wy = C, kyKry. 

As we established earlier, in the absence of a current the increment 
attains its maximum for œ, ~ wy, ie. for k ~ ki xry(Q,t,)*. The effect of 
the current becomes appreciable if at the maximum the value of k,u/kirj, is 
K 
ky 
exceeds this value. Then the point of maximum growth rate shifts towards 
large k,, since the value of uk, increases with k,, and the actual maximum 
value slightly increases compared with the currentless plasma. 

For sufficiently large k, values, when œw, > k,u we have approximately 


ik 
ox a,(1+ ee ) = a(143 : ). (IV. 96) 


wskiry Vek Ae 
where A, is the electron mean free path. We must assume k,/, < 1 because 
otherwise the diffusion approximation could not be used for the longitudinal 
motion of the electrons. Equation (1V.96) shows that as k, decreases the 
growth rate increases. For small ryk, this increase continues until œ? attains 
values of the order w,k,u/k{rj, and then the growth rate begins to decrease 
with k, The maximum value of the growth rate, of the order of 


2 \4 
Yı = uk (- a (IV. 97) 


; kp; (uwa.e\* 
is reached for k, ~ E ew : 


comparable with œ, i.e. for u > ue = ¢,rH k: J Q.t,. We assume that u 


ve 


e 
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For not very small values of k,r, it may be found that the value of k,u 
becomes larger than w, before we reach this value of k,. In this case we can 
neglect w, compared with k,u in the dispersion equation (IV.95), and then the 
maximum value is 


y. = 0,/kir}, (LV. 98) 


for k,u ~ w,/kiry. Thus, in the general case y = min (y1, Y2) > @,. For 
k, ~ « eqn. (IV.97) can be used if 


Me ? 1 e $ 
FyK < (z (Ae KY a (IV. 99) 
M; T 


(h) Transition to Collisionless Dissipation 


To conclude our discussion of the drift-dissipative instability, we shall 
briefly consider the transition from dissipation due to collisions to collisionless 
dissipation (90). Using the simplified expression for the collision term, 
assuming that after scattering at the neutrals the electrons have a Maxwellian 
distribution, the linearised kinetic equation for the electrons can be written 
down in the following form 

Of c. of 


e n 
Ok) i kog mi eo = = | 
(= 10+ il) fay +5 ikea 5 Te Egy Pro = —Yehuot vez S (IV. 100) 


where v, = Te ', and from this we obtain the electron density perturbation 


: W-W,y,  (M+iv, ive . (mtiv.\]~') eg 
stan arent re | hier ced ow) I 
where the function Y(z) is defined by eqn. (IV.49). 

For v, > 0 we recover the result previously obtained for a collisionless 
plasma. In the other limiting case t,k,v, = 4,k, < 1l and wt, < 1 we obtain 
the hydrodynamic result ([V.85). The transition from collisional to collision- 
less dissipation occurs broadly speaking at 2,k, = 1. A more exact treatment 
of this transition (90) shows that oscillations with œ ~ v, lead to some 
broadening of the region of the ion-sound instability. This additional 
region is marked on Fig. 23 by the oblique hatching. 

At very small k, values the function Y can be expanded in inverse powers 
of the argument. For wt, > 1, retaining the first two terms of this expansion 
leads to the expression 


He w — Wy dO oe l ep 
—= =] — 1 + ————_ ] >= IV .102 
n í @ | t alot ivm, E : i 


e 


Comparing this expression with the ion density perturbation which, for 
drift waves with T, = T, and œ/k, > v; 1s given by 
n; & ae ep 
= | —- —5S ---- — 


n w w T 
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according to eqn. (1V.66), we obtain the dispersion relation 


W—-Wy, lw + ive) 


= 2 
W+ Wy w; 


(IV. 103) 
2 


where %2 = Ta Ues ve =t, '. We have found earlier (see eqn. (IV.57)) 
L 


w 
that in this equation complex roots appear for ve = 0 and P > 0-4. Inthe 


presence of collisions the region of instability expands towards large k, 
values, when œ, > w,, and we have approximately 
2 
© © ostiv S (IV. 104) 
W, 

This expression is similar to eqn. (I[V.92) which was obtained in the 
hydrodynamic approximation, and for 7; = 0 it can be shown that the factor 
2 in the second term of (IV.104) disappears and the two expressions are then 
identical. 


(i) Current-convective (Spiral) Instability 
We shall now show that the instability of the positive column described in 
(1.2.b) transforms into the drift instability as the collision frequency decreases. 
Consider the simplest case Q,t; > 1. We can then use eqn. (IV.89) for 
the perturbation of the ion density and (IV.94) for that of the electron 
density. Comparing n; and ną we obtain the dispersion equation 


ik,u 
w? +iw(v;+ D,k2+0,+ik,u)—v(D,k?2 +ik,u)— iw, Q ae i 4 = 
i'n 


(IV. 105) 


For u = 0 this equation goes over into (IV.90) and for v; = 0 it coincides 
with (IV.95). For ryk < 1 the term k?.D, can be neglected compared with o, 
in the second term of (IV.105) as well as in ([V.90). In addition, for kry < 1 
we can neglect k,u compared with œw, in the second term, and in the last term 
w, compared with k,u/k{rj. In this approximation for œ < œ, we obtain 
from (1V.105) 

_ kto”; jOxk.ulkirn—ViDe ki (IV. 106) 

w tyi w, t Vi 

For sufficiently small k, values the first term in the expression for the imagi- 
nary part of the frequency is larger than the second, and the perturbations 
increase with time. This constitutes the current-convective instability which 
originates from the gradient of conductivity across the magnetic field. 
According to (IV.106) the expression for the growth rate can be represented in 
the form 


b,dinnk, œ b k o 
seep E a es a E E 
y i b, dx ky w+ y; e b, Q?7? O+ y: (IV. 107) 
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where b; is the ion and b, the electron mobility, E the longitudinal electric 
field, related to u by u = b,E. For perturbations with not very large k, 
values, the second term in (IV.107) can be neglected, and the first attains a 
maximum at œ, = v; given by 


dinn 
=U 
y Tx (IV. 108) 
where 
b. 
U = 4b,E J 5 (IV. 109) 


(j) Non-linear Drift Flows 

If the conditions for the drift instability are fulfilled, small perturbations 
will increase with time until the non-linear interactions come into play. To 
help to visualise at least qualitatively the character of these non-linear oscil- 
lations, we can reduce the equations for the non-linear motion to a single 
equation of the hydrodynamic type. We shall discuss here the derivation of 
equations of this type for non-linear flows corresponding to the drift-dis- 
sipative, current-convective and flute instabilities. 

Let us start with the drift-dissipative instability in a currentless plasma. 
For simplicity we neglect the collisions between the ions and the neutral gas 
atoms, and assume that T; = 0. In these oscillations with a characteristic 
scale parallel to the magnetic field considerably larger than that across the 
field, the longitudinal motion of the ions can be neglected. The transverse 


motion is principally governed by the electric drift v = T [hV¢), with a small 


correction due to inertial effects. In the continuity equation this correction 
can be neglected, and using the incompressibility of the drift velocity, div 
v = 0, we obtain 


a +vVn =0 (IV. 110) 


Furthermore, summing the electron and ion equations of motion, we 
obtain the hydrodynamic equation 


d 1 
mn +Vp =—[jH] (IV. 111) 
dt C 
where j is the electric current density and p the plasma pressure. 
Applying the operation (curl), to eqn. (IV.111), we obtain on the right hand 
s RR , a H oj, 
side the expression =. divj,, which because of div j = 0 equals z 
The value of j, can be determined from the electron equation of motion, 
neglecting the longitudinal ion current, giving 


2 al 
i= (2 : S (IV. 112) 
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We thus obtain 


d 0 0 
curl, (mn) = -gz eTenz_ (eo— Tin n) (IV. 113) 
, | | c[hV¢ | 
Equations (IV.110), (IV.113) together with the relation v = ——__— are 


H 

then the required equations of non-linear motion. From the hydrodynamic 
point of view, the potential @ can be considered as a function of the current 
of two-dimensional motion (v, = 0). Equation (IV.110) describes the 
variation of density in each plane z = const, and the right hand side of eqn. 
(IV.113) establishes the connection between the flows in the different planes. 
It is easily verified that in the linear approximation these equations lead to the 
dispersion eqn. (IV.91). 

In the special case of a helical flow where all quantities are functions of 
only two variables, namely the distance r from the symmetry axis and 
€ = 0— kz where 8 is the aximuthal angle, eqn. (I[V.113) can be integrated 
once, giving 

dv Ho x35 cT ð 
m;n Ti +Vp,= 3 aa eo oH? ae” n) (IV. 114) 


is the conductivity of the plasma, e, the unit vector directed 


NT, 


e 
where o = 
e 


along the radius and p an arbitrary function of the co-ordinates. Equation 
(IV.114) together with the continuity eqn. (IV.110) and the incompressibility 
condition 


divy = 0 (IV. 115) 


describes a hydrodynamic helical flow during which the fluid is subjected to 
an additional force determined by the right hand side of eqn. (IV.114). The 
second term in the expression for this force leads to the instability. 

We can similarly obtain equations for the non-linear flow corresponding 
to the current-convective instability. To this end we need only replace the 


Sexi ioi ; Ô 
inertia term by the frictional force m,nvv; in eqn. ([V.111) and add to — = 
the external electric field E in the expression for 7. The result is to replace 


(IV.113) by 


0 ĝ ô 
curl (m; nvyv;) = — gz eten (2 — eE— sir In n) (IV. 116) 
and for the helical flow we then obtain 
H?’ c cE Tc @ 
m,nV;V+Vpy= — ——9r7?k? — —— — — —] 
nv.v+ Vp. x r e, (v, Hrk  eHr ov an) (IV. 117) 


For small k the dominant term is the second term in brackets, which for 
Ejk > 0 corresponds to a force acting in the radial direction. Neglecting 
the other terms, we obtain the equation of motion for an inhomogeneous 
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incompressible fluid in a porous medium in the presence of a radial gravita- 
tional force. 

Finally let us consider the case of a flow corresponding to the flute 
instability including the effect of the finite Larmor radius. For simplicity we 
again replace the diamagnetic expulsive force by an effective gravitational 
force with acceleration g. Since the velocity of the electrons and of the ions is 


again determined mainly by the electric drift v = ees we again have 
a two-dimensional incompressible flow. This flow can be described by 
the magnetohydrodynamic equations including a collisionless ‘‘ oblique ” 
viscosity (see for instance (99)). It has been shown in (100), (101), that in the 
linear approximation these equations lead to exactly the same results as are 
obtained from the kinetic equation for k,p; < 1. 

For Q;t; > 1 the equation of motion assumes the following form 


dv 
m;n— T +Vp, = m;ng+E (IV. 118) 


where the force F is given by the expression 


cT T 
F = TAA Ag—(VnV) Vo) = ra {Vn cur], v + (VnV)[hv],} (IV. 119) 


be ues H ] 
and the effective pressure p, is given by P= 2nT + gaT nT n A 
ja i 


For the incompressible flow considered here, p, can be taken to be an 
arbitrary function of the co-ordinates r, 0. 

Now it is the force F which leads to the stabilisation of small flute per- 
turbations. However, there need not be such a stabilisation for finite per- 
turbations. For it is obvious that an isolated plasma tube in a vacuum, at 
a point where g # 0, must be ejected in the direction of g. Such a motion 
is similar to that occurring in the m = 1 mode; the electric field is constant 
over the tube cross-section and there is no difference between the motion 
of the electrons and of the ions. One might think that such a tube cannot be 
in equilibrium when it is immersed in a plasma of considerably smaller (or 
considerably larger) density. In fact we assume the opposite, that a steady 
state flow is set up in the plasma which does not lead to any flow towards the 
walls. Such a flow is stationary in some moving system of co-ordinates. In 
this co-ordinate system the equipotential surfaces œ = const must coincide 
with surfaces of constant density n = const, since (vVn) = 0. Considering an 
individual plasma tube bounded by a surface n = const, we can show that 
the integral of the force F over the volume of this tube vanishes. For instance, 
for the integral of F, we obtain by integration by parts 


cT a7@ 67 
o <£ Rs id IV. 120 
| Feat Ha,” (=? 7 dx BE pay) ( ) 


where the integration on the right is taken over the surface n = const. 
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Therefore we can remove n from the integral, and obtain the integral over a 
0 
oy 
we can show that the integral of F, also vanishes. Thus the total additional 
force acting on a separate tube in a plasma bounded by a surface disappears. 
It follows, therefore, that such a tube can only maintain itself in equilibrium 
for certain special flows when the pressure p, is distributed over its surface 
so that it exactly balances the frictional force. 

From these considerations it would seem that for sufficiently strong per- 
turbations, when plasma tubes become separated from the bulk plasma, and 
the surface n = const becomes multiply connected, stabilisation due to 
finite Larmor radius may be absent. 


closed path of the complete differential a ), which must vanish. Similarly 


4. TURBULENT DIFFUSION OF A PLASMA 


Since the only reason for the drift instability is the inhomogeneity of the 
plasma, oscillations developing in consequence of the instability cannot die 
away until the inhomogeneity is completely destroyed. In other words, such 
an instability must lead to turbulent diffusion. 

The magnitude of the “‘ diffusion ” flux brought about by the oscillations 
so developed can be determined from the following considerations (103, 102, 
212). As we established earlier, in the case of a drift instability (raK < 1) the 
ions move across the magnetic field mainly due to the electric drift 


v- c[hVo] 
= H 


In a homogeneous magnetic field the drift motion is incompressible, div v = 0, 
and in the presence of drift oscillations the plasma moves across the magnetic 
field without additional compression or dilation, i.e. turbulent convection 
takes place. In this case the displacement of the plasma by a distance č 
leads to a perturbation of the density n’ = ¢xn. In other words the density 
fluctuations are determined by oscillations of the displacement č. If these 
oscillations vary harmonically, they do not lead to any net flow of the plasma 
when averaged over time. On the other hand, if the amplitude of the oscil- 
lations increases with time, then each succeeding half-period of the oscilla- 
tions leads to a slightly greater displacement of the plasma than the previous 
one, and as a result a mean plasma flux q = <én'> = yxn<é?> occurs. 
To order of magnitude this relationship between the flux g and the dis- 
placement of the plasma € or the density perturbation n’ can be retained 
when the interaction between waves comes into play, since even in the pre- 
sence of an interaction the characteristic rate of growth of individual wave 
packets and of their transformation into other packets is of the order of the 
linear growth rate y. 

Consequently, to determine q, we need to know <é€?7> or <n’?>. 
For y~ œw, when strong turbulence develops, the oscillation ampli- 
tude increases to such an extent that the perturbation of the density 
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gradient k,n’ becomes of the order of the mean gradient, i.e. k,n’ ~ xn. 
In this case, therefore, the steady state intensity of the density oscillations 


with wave number k is of the order of magnitude ny |? ~ 


k? 
2 dn 
qq= — = <n >= 3) Kn = 
dx 


For weak E when y/@ < the or wave equation can be 
used to evaluate the oscillation intensity. If the dispersion relation is of the 
decay type, that is if the condition w, = wo, —@,_, can be satisfied, not 
identically but at some surface in wave number space, we have to order of 
magnitude y < n’*> = V*<n’'*>?. In this relation the left hand side 
describes the development of the waves due to the instability, and the right 


hand side represents the effect of the interaction between the waves. The 
2 


; k 
square of the matrix element can be estimated as V? ~ O73 3, since for 
a7 


n ~ = when the gradient of the density perturbation is of the order of the 
mean gradient, the interaction between the oscillations ought to lead to a 
damping rate of the order of the frequency. Thus for weak turbulence 


the intensity of the oscillations must decrease as the increment decreases, and 
2 


K 
to order of magnitude the quantity <n’? > ~ ta mr Approximately, 
y? \ dn 
wk?/ dx’ 
in this relation represents the turbulent diffusion coefficient. To estimate the 
value of the diffusion coefficient, we can use for k, the value for which the 
diffusion coefficient is a maximum, and write to order of magnitude 


Yk 
D~ IV. 121 
(3 a l ) 

A similar relationship holds good for the ion sound instability (rąk > 1), 
since analogous considerations can be applied to the electron drift motion 
across the magnetic field. Indeed the case ryk > 1 seems simpler since both 
in the absence and in the presence of collisions y ~ œ ~ kc, and consequently 


D ~ c¢,/k 
In the presence of a strong magnetic field, when xr, <1, the turbulent 
diffusion is more sensitive to collisions, and more detailed consideration 1s 
therefore necessary and is given below. 


therefore, gq = — The factor multiplying the density gradient 


(a) Rarefied Plasma in the Absence of a Longitudinal Current 

Let us consider a plasma of such a low density that collisions can be 
neglected. More precisely, we assume that collisions on the one hand are so 
rare that they have no effect on the frequency and growth rate of the 
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oscillations, whilst on the other hand they occur sufficiently frequently to 

maintain the Maxwellian velocity distribution of the particles. The coefficient 
of diffusion corresponding to these assumptions will be designated by D,. 

To evaluate D, we use results obtained in Section 2(c). Let us first 

Srp m 

consider a plasma column of unlimited length. For B= P a 

i 

only perturbations with kp; Z 1 are unstable. As a function of k, the growth 

rate reaches a maximum at k, ~ "N B. In this case œw ~ kxp,v;, and 


Cy mM, 
v/@ pepe SS 
Ve m; 


(we assume T; = T,). It can be shown (see below (f)) that the main contribu- 
tion to the diffusion comes from oscillations with k,p; ~ 1, and in this case 


m m 
D, ~ — p?v;k ~ — Kp; D IV. 122 
S mp i mp oe B ( ) 
where Dz ~ pw; is Bohm’s coefficient of diffusion. 
Thus in this range of $ the diffusion coefficient decreases as f increases, 
and for B = 107+ the plasma becomes stable and D, vanishes (leaving only 


the classical diffusion due to collisions). 
For $ < a we may have, in addition to the short wave instability with 


kip; ~ 1, a “ hydrodynamic ” instability for perturbations with large 
wavelengths for œ% = w? (see (IV.57)), and these perturbations make the 
greatest contribution to the diffusion coefficient. Since œ ~ @, < c,k,, the 
minimum permissible value of k, is ~ @o/c and we obtain, using (IV.121), 
the corresponding diffusion coefficient 


Dew D w p;D IV. 123 
= Kp; . 
S B m,B i*~B ( ) 


Wo 
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Fic. 24. Coefficient of turbulent diffusion in a tube of infinite length 
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In this region of magnetic field the coefficient of diffusion is inversely 
proportional to H, but in absolute value it is smaller than Bohm’s coefficient. 
As the plasma density decreases, when wy, becomes of the order ck, i.e. for 


P~ = x’ p?, the diffusion coefficient reaches its maximum possible value of 
i 


the order Dp, and does not change when the density further decreases. In 
this case y ~ w ~ w@,, and the maximum scale of the turbulent pulsations 
is determined by the tube radius, i.e. ky ~ K. 

The relationship between D, and $ for an infinitely long column is shown 
schematically in Fig. 24 by the full line (on a double logarithmic scale). 

Now let us examine the effect of shortening the plasma column on the 
diffusion. The value of the wave number x,,,, at which the growth rate 
reaches its maximum, increases with f, and is given by 


m m 
kim ~K Je k?p? for B<—x’p? 
m; mM; 


m 
kim ~ KB. /— for = K pi < B< — 
i 


= m 
and kim ~ KNB for po 


i 


Therefore the effect of finite length of the plasma column is more important 


Mmi 
for low density plasma. This effect starts at L < a J P (kp)? where L 


is the length of the column and a its radius. For smaller values of L the dif- 
fusion coefficient decreases from its value for infinite length, and is given by 


m,\? /« \3 a ko Jo 
Ds (= (=) oe MOn L k m; ( ) 


m.(K\? a J m 
> om, (=) par for L m; ( ) 


where ky = 2zn/L is the minimum longitudinal wave number. 

Thus, in a tube of infinite length, as $ decreases the diffusion coefficient 
is at first equal to the value for the infinitely long column, but when the value 
of k,,, reaches ko, the diffusion coefficient remains constant and does not 
increase as $ decreases further. This is shown in Fig. 24 by the “ dashed ” 
line. 


and 


a Mme , l 
According to (IV.124), for L < J Ks the coefficient depends relatively 


a Me 
~ J — and small f values the 


weakly on the length of the device: for L a 
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coefficient of diffusion is equal in order of magnitude to D, ~ xpjv,; for 


a me | | 
L > J — the coefficient D, decreases fairly rapidly as L decreases, and 
m; 


a 
L 
perturbations with k, = 0 which are sensitive to the boundary conditions 
at the ends of the tube). 


for — = 107! complete stabilisation occurs in general (we do not consider 


(b) Low Density Plasma carrying a Longitudinal Current 


a e 
Consider a plasma column of limited length, say =- ~ J m We 


i 
have shown earlier that in the absence of a longitudinal current only short 
wave oscillations with k,p; ~ 1 are excited in the plasma. The corresponding 
diffusion coefficient is comparatively small, namely D, ~ kp;D,. In the 
presence of a longitudinal current, as has been shown in Section 2(b), an 
additional instability appears at long wave perturbations with k,p; < 1. 
Since the diffusion coefficient ([V.121) is proportional to the square of the 
transverse wavelength, this instability increases its magnitude considerably. 

Let us once more consider a plane inhomogeneous plasma layer with 
density gradient along the x axis and assume that T; = T, = const. In 
equilibrium or, more precisely, in a state averaged over time, the electric 
field in laboratory co-ordinates, in which the ions are at rest, is determined 
by the density gradient 


T 
Vo = ——Vn (IV. 126) 
en 


The longitudinal electron current excites oscillations of the potential o’ 
and density n’. According to Section 3(b), waves with longitudinal phase 
velocity considerably smaller than the electron thermal velocity are excited 
in the plasma and for these waves all electrons except the resonance electrons 
can attain equilibrium parallel to the magnetic field. Therefore the perturba- 
tion of the potential can be represented in the form 


T n’ 
p’ =—In (1 + 4 +, (IV. 127) 
where the first term represents the Boltzmann distribution, and the small term 
@~, is related to the resonance electrons. 

If the Larmor radius of the ions is appreciably smaller than the charac- 
teristic wavelength of the perturbations considered, the hydrodynamic 
approximation can be used for the ions. In this approximation it is sufficient 
to consider only the electric drift in the continuity equation, since the Larmor 
current does not lead to a change of density 

Cn’ 


i + 7 [hV(9+¢’)] Vin+n') =0 (IV. 128) 
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Using relation (IV.127), this — can be written in the form 


On’ 2cT 
ees < [bvn] Vn +— = [hV] Vin +n )= (IV. 129) 


Ot 
For ọ, = 0 this equation is linear in n’. This circumstance, which indicates 
that there is no interaction between the modes for a Boltzmann distribution 
of the electrons along the lines of force, is very important for all that follows. 
In eqn. ([V.129) the second term describes the transfer of the perturbation 


| dn 
along the y axis. For sa const the velocity of this transfer 
I 


is also constant, and by transforming to a moving co-ordinate system, the 
second term in (IV.129) could be eliminated. However, for a general distribu- 
tion of the mean density n(x) this velocity is also a function of x, and because 
of the resulting differential motion the density perturbation is deformed 
with time. To demonstrate the character of this deformation, we shall go 


over to a spectral representation. We put n’= f Noe t dk do, 


where Mgko(r, t) is the slowly varying amplitude of the wave packet k, œ. 
For k > x, so that the wavelength of the perturbation is considerably smaller 


dn 


See, : 
than the transverse dimension a, the quantity P can be expanded in 


series near the point x = xọ under consideration, and we shall retain the 
first two terms of this series. Since in transforming to the Fourier representa- 


ri 
tion x— Xo > i=, we obtain 
° ôk 


x 


ON key dw. hko = 
— i(w—20.)ny + a +2 rp TA = Tae g, Oath — 
Vn u 
mD Al Oko k-kr' o-a dk’ dw’ =0 (IV. 130) 
n Ue 
cTk, dn 


—, and u is the longitudinal mean electron 
eHn dx 


where œ. = k, vo = — 
velocity. 

In eqn. (IV.130) the last two terms, the first of which describes the build- 
up of the oscillations by the resonance electrons, and the second the non- 
linear interaction between the waves, come from the Fourier transform of the 


last term in (IV.129). For gy, we insert its approximate value 


eee IV. 131 
Piko = i r- v, n ( ) 


which can be obtained from the expression for the density perturbation 
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(IV.50) for w/k, < ve Strictly speaking, the expression (TV.131) should con- 
tain an additional factor k,/|k,|, which we have omitted, assuming k, > 0. 
In eqn. (IV.130) the first term is considerably larger than the rest, and 
therefore as a first approximation we obtain œ = 2@,. In other words, the 
function n,e is close to m,6(@—2@,). The remaining small terms in ({V.130) 
describe the evolution of the wave packet in time. In particular, the term with 


ke shows that the deformation of the wave packet due to the differential 


ved e e e 
drift velocity vo leads to an increase of the wave number k, with time, Le. 


to a flow in wave number space. To order of magnitude this term is equal to 
K 
Ws n 
k, ko 
and for small k, in particular it will be the decisive term: all perturbations 
will “ drift ? into the region of large k, faster than their amplitude grows 
k ; = ; ; 
due to the instability. For — > 2e the instability begins to be dominant, 
K u 
and the distortion of the waves can be neglected. Neglecting the second 
term in (IV.130) and introducing the new variable v = w—2a,, we obtain 
for the region k/k, < ufve: 


VNky — ino. = ny,— i = =f Ngy Ngy Ak’ dv’ (IV. 132) 
Ve n Vve 
where k” = k—k’ and v” = v—v’. (We neglect here the non-linear Landau 
damping at the ions, assuming u > 1t;.) 

In the linear approximation, it is clear from eqn. (IV.132) that the charac- 
teristic frequency v is purely imaginary, so that we are concerned with an 
aperiodic instability which must lead to strong turbulence. This can also be 
seen directly from eqn. (1V.132) according to which the non-linear interaction 
can balance the linear growth of the oscillations only when the perturbation 
of the density becomes of the order of the mean density. This result is a direct 
consequence of the absence of interaction between modes, noted above, fora 
Boltzmann distribution of the electrons. 


Let us estimate the value of the diffusion flux q -7 n') We 


. 0 , 
must substitute — for E, since only the resonant electrons undergo 


diffusion (it is easily seen that the term containing <7 In (1 +=) 
ye n 


vanishes when averaged against n’). Using (IV.131) we obtain 


CE ft TI kN ia kd n 
= — V/t— 
eH Ve LA n AR GV. 133) 


q 
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Where Nko 1s the spectral function of the density: 
CNkoy Ne’) = Ngo d(@ —@ ') ô(k —k’ ). 
In the case of strong turbulence the density perturbation is of the order of 


cT u 

the mean density and consequently to order of magnitude q~ Ti kin, 
eH v, 

where k, is a typical value of the wave number in the region of strong 


turbulence. As we have shown above, this region begins at k, ~ K— 
v 


X 

Substituting this value of k, into the expression for q (since the effective 

value of k, is in any case not smaller than this value), we obtain to order of 
cT dn 


eH dx 
Thus, in the presence of a longitudinal current the diffusion coefficient 
of the plasma in a very strong magnetic field is of the order of the Bohm 


magnitude q ~ — 


value.t This result is valid only when the value of * is not excessively small. 


Ve 


In fact, for : ~ J A as well as the build-up of the long wave perturba- 
M; 


tions considered here, there is an instability at perturbations with kip; ~ 1 
which also leads to diffusion. Since the diffusion coefficient determined by 
the short wave perturbations is of the order of xp4v,, the long wave perturba- 
tions must be damped at a rate of the order of ~kixp?v;. Clearly the results 


given above are only valid if for perturbations with k, ~ K- the growth rate 
u 
yr 2 w, is larger than this damping rate, i.e. 
Ve 
u —— 
—>/xop, (IV. 134) 


Ve 


and for smaller values of u/v, the coefficient of diffusion must drop quite 
sharply to a value ~ «p;Dp, since the differential rotation of the plasma 


(F Æ 0) shifts the long wave perturbations into the region of large k, 
x 

where they are damped by diffusion arising from the short wave perturbations. 
In the presence of even a small shear of the magnetic field, of the order 
0 ~ kp; the short wave oscillations are stabilised, and then the turbulent 
diffusion coefficient may rise to ~ Dp even for small values of u/v,. 


r : , CE n’ 
+ In ref. (103) we assumed for the velocity fluctuation v ~ 7 ~ aa too low a value, 
K 


; u , 
namely v ~ yn’/kin, so that the result obtained there, D ~ 3 Dz, is also too low. 
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(c) Diffusion of a Very Low Density Plasma 


The results. given above are also invalid unless collisions can in fact 
maintain a Maxwellian distribution of the longitudinal velocity of the 
electrons. In a very tenuous plasma this may not be possible, and then 
conclusions based on the specific form of the distribution function may need 
considerable revision. A decrease in the collision frequency must have a 
considerable effect on the diffusion either in the presence of a longitudinal 
flux, or for B > m,/m,, where the oscillations are built up by a small number 
of resonant electrons. For the case p < m,/m, on the other hand, where the 
instability in a long tube is of a hydrodynamic type, and all electrons partici- 
pate in the growth of the oscillations, the results are not sensitive to reduction 
of the collision frequency. 

Consider first the case of diffusion due to longitudinal current. For 
the preceding results to remain valid, the distribution function must be 
similar to a displaced Maxwellian distribution with an accuracy up to u/v,. 
On the other hand, the diffusion loss leads to an appreciable distortion of the 
distribution function in a time of the order of the characteristic diffusion time 
k?D ~ x?Dz,, and this distortion can be destroyed by collisions only if 


u . ; i 
ve— > K*D,. Introducing the parameter S = 1.p;k?, this condition can be 


e 


written in the form: 
S<u/v, (IV. 135) 


In a very tenuous plasma this condition is infringed, and a plateau 
develops on the distribution function. The escape time is then determined 
by the rate at which collisions can restore the Maxwellian distribution, and in 


l . u 
order of magnitude this cannot exceed —v,. 
v 


e 
Turning to the case of a currentless plasma, for B > m,/m, we also find 
that only a small number of resonant electrons, with longitudinal velocity 


m nn : TE : 
U.<C,4 = 0, J —< participate in the diffusion. A decrease in the col- 
m;ß 


lision frequency then again leads to a decrease in the diffusion coefficient. 


(d) Difusion of a Dense Plasma 


The expressions given in Section 4(a) above for the turbulent diffusion 
coefficient D, refer to a fairly tenuous plasma where the collision frequency 
Ve is smaller than œ ~ w, for all waves including k, ~ x. This condition 
can be written in the form 


m; 


S> (IV. 136) 


Me 


For smaller values of the parameter S = A.p;K* the friction between the 
electrons and ions due to collisions becomes important, leading to the 
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drift-dissipative instability. The collisions lead to an increase in the optimum 
wave number k,, at which the growth rate reaches its maximum, and from 
the condition œ, = w, we obtain fork, =x 


m,\* 0; Kp; 
kız (| =) x*o,,/—=7=—— | 
(=) SPUN 2, NVL D 
The longitudinal phase velocity corresponding to this value of k, is of the 


+ + 
order of w/k, X v; (=) @ , and if this velocity is smaller than the Alfvén 
velocity c4, and k, given by (IV.137) is larger than the minimum possible 
value ky = 2n/L, then y ~ œ ~ w for ky ~k, and the corresponding 
coefficient of turbulent diffusion attains the value Dg ~ pw; According to 
(IV.137) the condition w/k, < c, can be written in the form B < (Qta) +, 
and if this condition is violated we can have y ~ œ only for larger values of 
k,. The diffusion coefficient then decreases to 


ee (IV. 138) 


where D, = p?/t, is the classical coefficient of diffusion. At the same time 
the optimum wave number increases: k, ~ x*p,$?Q,1,. 

In a device of limited length, when the minimum wave number ky = 22/L 
is larger than the optimum values given above, the coefficient of diffusion 
is determined by shorter wave perturbations, for which again y ~ w (the 
contribution from the long wave perturbations ~ y*/wk? varies as k? as k, 
decreases). The corresponding coefficient of diffusion can be estimated as 


ote IV. 139 
D x D,{—== 
(Janz) ana 


+ 

This expression in its turn is valid only for kọ/K < (=e) A,K, Since 
i 
otherwise the optimum wave number k, becomes so large that viscous 
damping of the oscillations ~k{p?/t; becomes important. Moreover, 
expression (IV.139) is valid only for kgA,<1, where the diffusion approxi- 
mation can be used, while if this expression gives a greater value than (IV.138) 
the latter must be used. 

We shall now consider the effect of a longitudinal current. We have 
shown above that in a tenuous plasma the diffusion coefficient in the presence 
of a longitudinal current may reach the value Dz. This conclusion requires 
k,4, > 1, so that collisions can be neglected. Since the maximum possible 


value of k, is determined by w/k, ~ = ~ v, the condition for the validity 


of the collisionless approximation 1s 
S>1 (IV. 140) 
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For smaller values of S = 1,x*p, the diffusion approximation can be used. 
According to Section 3(g), in a strong magnetic field («*p7 < 1) and with a 


sufficiently large longitudinal current (u/c, > x2 p? Q.T.) the growth rate of 
small perturbations (IV.97) is larger than the oscillation frequency ~@,. 
The corresponding coefficient of diffusion can be estimated as 


p~e g ) ~ Dal) (IV. 141) 
K ACU AK Ul» 


ue = C KFN QT, (IV.142) 


where 


Relation (IV.141) refers to a plasma of fairly low pressure (8 < (Q,1,)~*) 


UC, Kp? 


272 
e’*c 


greater than k,, the diffusion coefficient varies inversely as kọ and is given by 


4 
in a long device (o<k. = x( l As ky increases and becomes 


; 4 
b~p, (=) ae s (IV. 143) 


B 
ko \Us vk. 


This relation can be used as long as the condition kod, < 1 is satisfied, or 
provided that the value given by (IV.143) is less than that given by (IV.139). 
According to eqns. (IV.141) and (IV.143), the diffusion coefficient in a long 
tube carrying a current such that u > ux is independent of the magnetic field, 
while for a shorter tube it decreases as H™+. 

According to the results of Section 3(g), the longitudinal current may also 
increase the diffusion coefficient in the case of a weak magnetic field ryk > 1, 
where an 10n-sound instability occurs. 

All these estimates for the turbulent diffusion coefficient are very rough, 
and for systems of limited length they may not be applicable because we have 
not considered the possibility of perturbations uniform along the length of 
the system (k, = 0), which require special consideration for each specific 
case, since they must be sensitive to the boundary conditions at the end of 
the tube. 

Since turbulence developing from the drift instability is usually strong, 
to obtain a more accurate determination of the turbulent diffusion coefficient 
it would be necessary either to use the weak coupling approximation, or to 
introduce the free mixing length. We shall consider below two special cases 
for which a more detailed calculation has been carried out, one a case of 
strong and the other of weak turbulence. 


(e) Turbulent Positive Column (17) 


Consider an infinitely long positive column in an insulating tube in a 
longitudinal magnetic field, much higher than the critical field. We established 
in Section 3 that such a column is unstable to perturbations which are greatly 
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extended along the magnetic field, with a growth rate maximised over k, 
given by 
dlnn 


=U 
7 dx 


(IV. 144) 


E E E. 
for Q;t; > 1, where U =4b,E J b . This instability gives rise to a convective 
flow in the plasma during which the plasma escapes to the wall in the form of 
separate helical tubes, while outside the plasma isolated “ bubbles ” appear. 


Since at the boundary of these “ bubbles ” the relative gradient y is very 
n 


large, the instability causes them to move about fairly rapidly with the plasma 
and the motion assumes a chaotic and turbulent character. 

The magnitude of the diffusion flux in such a turbulent plasma can be 
determined by using the mixing length concept. Suppose / is the effective 
length or mixing length by which the plasma tubes are displaced by the 
convective motion, before the coherence of this motion is destroyed due to 
the interaction with other perturbations. Since the perturbation of the density 
arises primarily from the convective motion, the density fluctuation level is 


; d l ; l 
given byn = l = . The velocity fluctuation vu’, related to the density fluctua- 


, 


, , , dn : n 
tion by yn’ ~ v To can be evaluated as v = U —, and consequently the 
n 


diffusion flux q = <n’v’> is given by 


y 
h (IV.145) 
7 


q= -UI 


In a discharge in an insulating tube the walls have no stabilising effect 
on the pulsations, and by analogy with turbulent jets in an ordinary fluid the 
value of / may be assumed constant over the cross-section and proportional 
to the radius of the tube a. Using expression (IV.145) it is then straight- 
forward to set up a particle balance equation to determine the diffusion losses 
and the radical density distribution. The mixing length / can then be deter- 
mined by comparing the results with experiment. Such a comparison was 
carried out by the author (17), who obtained the following results for values 
of the magnetic field at which the turbulent loss rate is equal to the loss rate 
without a rnagnetic field: for discharges in helium Z/a = 0-15, in hydrogen 
lJa = 0:10, and in nitrogen //a = 0-12. Thus just as for the turbulent jet, we 
have l/a ~ 0:1. 

Figure 25 shows a comparison of the experimentally measured and 
theoretically calculated dependence on the magnetic field of 0 = E,/Eo, 
where E, is the electric field in the turbulent discharge, and Ep the electric 
field for H = 0. The good agreement between the experimental and theo- 
retical curves shows that the description of the turbulent plasma on the 
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basis of the mixing length concept is fully justified. Moreover, experi- 
mental measurement of the radial density distribution (104) confirms the 
validity of the assumption / = const. The theoretical and experimental 
density distributions are compared in Fig. 26. The boundary condition 


Fic. 25. Dependence of ratio of the longitudinal electric field E, in a turbulent discharge 
to the field E, for H = O on the neutral gas pressure, for discharges in helium 


Fic. 26. Comparison of theoretical and experimental radial distribution of plasma density 
(1) Discharge without magnetic field. (2) Turbulent discharge 


q = n,U was assumed in calculating the theoretical distribution where n, is 
the density at the wall. This boundary conditon is analogous to the introduc- 
tion of the extrapolated length in ordinary diffusion. Figure 26 shows that the 
actual density distribution is closely similar to the theoretical form and very 
different from the distribution in the absence of the magnetic field, shown on 
the same figure. 
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(f) Turbulent Diffusion of a Tenuous Plasma with B > ie 
ni; 


The second example to be considered is of interest in that it shows how 
peculiar may be the interaction between the oscillations in a turbulent 
plasma. 

Consider an inhomogeneous low density collisionless plasma in a strong 
magnetic field and suppose p;x < 1, 8 > m,/m,. It was shown in Section 3 
that such a plasma is unstable to short-wave perturbations with kp; ~ 1. 
For $ > m,/m, the growth rate of perturbations with kp; ~ 1 is considerably 
smaller than the frequency, so that the motion of the plasma is weakly 
turbulent and can be described by the kinetic wave equation. Such an equation 
has been set up and analysed in (105) (see also (106)). Only the results of these 
investigations will be given here. 

The dispersion relation of the excited drift waves is of the decay type, and 
the decay of the waves leads to diffusion of energy in momentum space. In 
addition to the decay interaction, the non-linear damping of the waves due 
to beat oscillations with phase velocity w"/k? S v; is important for these 
waves. The matrix element of the wave interaction is small for modes with 
wave vectors of similar absolute value, but is appreciable for modes with very 
different wave numbers. Because of this, perturbations with k,p; ~ 1, 
though not strongly interacting among themselves, can nevertheless com- 
pletely suppress oscillations with k,p; > 1, and the spectral function has a 
distinct maximum at k,p; ~ 1. The amplitude of the oscillations of the 
potential and density near the maximum are determined by the balance 
between the linear build-up of the oscillations due to resonant electrons 
and the non-linear damping at the beats, and the shape of the spectrum is 
determined by the decay interaction. As a result of the suppression of the 
short-wave perturbations the value of the turbulent diffusion coefficient is 
determined by (IV.121) not at the maximum of this expression, which is 
attained for kip? ~ m,/m,B, but at k,p; ~ 1, where the value is 


Daci prv,K (IV. 146) 
m,; B 


| Tk{m,\* 

The potential fluctuations are of the order of ọ'~ — ae and 

e \m, 

m, \* a 
those of the density of the order of n'/n ~xp; (z . This is to be 

M; 
expected since the perturbation of the density arises from the displacement 
of the plasma by about a Larmor radius. 


(g) Electrical Conductivity of a Turbulent Plasma 


The classical diffusion of a completely ionised plasma is determined by 


120 PLASMA TURBULENCE 
the frictional force between the electrons and ions. Including this frictional 
force the hydrodynamic equations take the form: 


1 e 
Vp, = — enE — ~ [vA] + - [jH] + = (IV. 147) 


dy 
dt 


min | + Vp, = enE + c [vH] — j (IV. 148) 


NT, 


A e l 
where p, = nT,, pi = nT;, j = en(v—V,), 6 = is the conductivity, and 


e 


v the ion velocity. 
From (IV.147), (1V.148) we find that in equilibrium j =p [bYp], 


where p = p,;+p, is the total plasma pressure. If the plasma is inhomo- 
geneous along the x axis and E, = 0, we can determine the diffusion velocity 
due to collisions from the y component of either of the above equations. 


C: 1 dp 


Eo TEE IV. 149 
Jy MeTeQ2 dx ( ) 


v, = — 

When oscillations are present in the plasma, we must include in equations 
(IV.147) and (IVY.148) additional terms of the form + <en’E’>, where n’ 
is the density fluctuation, and E’ that of the electric field. These terms 
describe turbulent diffusion, and the y component in this expression leads 
to diffusion along the x axis. Formally they might be treated as an additional 
“ frictional force ” representing a decrease in the effective conductivity of the 
plasma. However, in practice enhanced diffusion cannot generally be 
attributed to an enhanced transverse resistivity. This would be possible if the 
reason for the instability were only the difference between the drift velocities, 
but as we have seen, the instability really arises from the inhomogeneity of the 
plasma, and the addition of a further velocity difference between the electrons 
and ions by the imposition of an external field is not generally equivalent to an 
increase of the density gradient. 

Nevertheless, the oscillations may have a considerable influence on the 
conductivity of the plasma across the magnetic field. It is well known that 
the effective conductivity of the plasma across a magnetic field depends 
strongly on the boundary conditions on the Hall current. This effect can be 
observed directly from eqn. ([V.147) for Vp, = 0, v = 0. If there are no 
constraints on the component of the current perpendicular to the external 

u 
(Qe Te)" 
Hall current is prevented from flowing, a Hall electric field, perpendicular 
to the external field, is set up in the plasma and the normal conductivity 
j, = cE, is restored. Any non-uniformity of the plasma conductivity 
makes the passage of the Hall current difficult, so that even a small non- 
uniformity of a strongly magnetised plasma may have a considerable influence 


field E, (the so-called Hall current), then j, = E,, but if the 
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on the flow of current along E. This effect has been demonstrated by Yoshi- 
kawa and Rose (107). 

Consider first the simplest case, of high frequency oscillations where 
the ions may be considered at rest. For simplicity we put o/n = const, so 
that the oscillations of the conductivity are determined by those of the density 
(corresponding in practice to a weakly ionised plasma). Then the equations 
for the fluctuations n’, E’, j’ and the mean values n, E, j, can be written in the 
form: 


1 
T,Vn' = —en’E—enE’ + -[j/H] + Zy (IV. 150) 
C O 


1 
T,.Vn= —enE-—e<n’E’> + A] + i (IV. 151) 
From (IV. 150) we obtain 
2 
j, = eS [H, T. Vn’ +en’E-—enVo'| + az (T,Vn'+en'E—enVg’) (IV. 152) 
o 


where Vo’ = — F’. 

Transforming to a Fourier representation (assuming x <% kı) and 
using div j' = 0, we obtain a relationship between the perturbations of the 
potential @, and the density nų 


epy(k*? + Q2t2ke — irk, QT.) 
n, 
= (T,k? + T,k? Q21? — ik, E e0212 — ik, E, Dnd (IV. 153) 


(we assume that E, = 0). Substituting this value for g, into the averaged 
eqn. (IV.151), we obtain an expression for the additional current 


.  , GE, 
Oj a~ Q272 


which for Q,te > 1 becomes 
ec 
oj, = — 7 n Ey 


k T, dn k2(k? + Q2tekZ)N,n? 
E A re 55 a r Ve tS 
O Í [Etat k, E, F en = (k? +0212k?)? +k?k QZT? ( ) 


and similarly 
6j, =- <n'E;) 
n 


o kz y T,dn krk Q.t hk? +Rtek:) Nun? dk 
D | fz were “zt ondx] (k?-+0212k2)? +K KOT? 
(IV. 155) 
When the inhomogeneity of the mean distribution is small, so that 
KQ,t. < k?, we can neglect the last terms in the denominators of the ex- 


pressions below the integral sign, and the expressions for the additional 


P.T. 5 
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current are then somewhat simpler. Assuming in addition that the spectral 
function of the density N, is isotropic, we can perform the angular integration 
in the above expressions, obtaining the simpler expressions 


_ y o mI) i 
a Epon] 6/,=—-yoE, (IV. 156 
J=% meal xt e ndx J pa ( ) 


where y = | Nun7? dk SL 


Thus in the presence of these oscillations the transverse current increases 
by a factor of almost Q,1,. In the case of non-isotropic oscillations the effect 
may be considerably larger. For instance, in the presence of perturbations 
which are greatly extended along the magnetic field, so that k, ~ Q7't7’, 
the transverse conductivity may reach values of the order ø. ‘This effect has 
been demonstrated experimentally (107). 

The effect of the inhomogeneity on the longitudinal conductivity is 
considerably smaller; the conductivity is only slightly less than the collisional 
value. Moreover, when the spectral function N, depends on the direction of 
the wave vector, the transverse electric field may excite a longitudinal current. 
This effect is described by the second term in braces in eqn. (IV.155). It 
occurs only when the perturbations are inclined to the magnetic field on the 
average, so that <k,/k,> #0. According to (IV.154), the longitudinal 
field E, then excites a transverse current. This effect is analogous to the 
appearance of an averaged particle flux in a turbulent positive column where 
the mean drift originates from the longitudinal electric field. 

The above results refer to high frequency oscillations (œ > Q,) where 
the ions are at rest. For low frequency oscillations the effect may be consider- 
ably smaller. In this case we must take into account the perturbation of the 
ion velocity v’ in the equations of motion, which become: 


1 

T,Vn = —enE—e(n'E’) — (n'[WH]) + - [jH] + Z j (IV. 157) 
1 , 1 e ! 1 °? en, 

T,Vn' = —en'E—enE' — A A H] + “Li H] + pa (IV. 158) 


0 = enE+e(n'E’) + <n’ [vA] = = j (IV. 159) 


d ’ , € , en e; 
mn, = enE +en'E + -n[vH]- —j (IV. 160) 


We can determine j, from the sum of eqns. (IV.158), (IV.160) and then 


use the condition div j = 0 together with the longitudinal component of eqn. 
(IV.158) to obtain 


i af „C œ À 
divj) = =i 3. (enk?¢,—ien, k E.) 


I 


o 
= 5 (T ny —nek2@, —ik,eE,1,) (IV. 161) 
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Expressing x in terms of m and using the linearised ion eqn. (I1V.160) it is 
straightforward to obtain the additional “ frictional force ” 


eH 
F, = <en'E,> + p n'v?’ 


, P c 
and from this the additional transverse current 6j, = Ar j 


n 
OF Ny 
Ok Fok O -2 dk (IV. 162) 


(we assume that the frequency œ is real). 


(kkZE,, I ky k, k*E,) 
Oj x = 


+ 
For isotropic oscillations this gives dj, ~ yoE, (a S ) , Which is 
idele 
considerably smaller than the value obtained above for high frequency 
oscillations. Even perturbations greatly extended along the magnetic field do 
not lead to a marked change in conductivity. At the same time oscillations 
of this type may lead to a considerable flux of charged particles 


f C / 
ds = Cogn’) & = Esn’) 
that is to turbulent diffusion. The smallness of 6j, and F, immediately 
shows that the diffusion flux is determined by the electric drift, i.e. 
cE, 
H 


From relation (1V.161), we can also determine the change in longitudinal 
current due to the inhomogeneous conductivity of the plasma; the result is 
2 


1 m 
Uy ~N 


0272k4E, + = Kk ekzEx 9, 
j, = —0 aa el z dk (IV. 163) 
277k? + g kt 


which shows that the effect is again small. 


5. TURBULENT PLASMA IN EXPERIMENTAL CONDITIONS 


During plasma experiments one encounters turbulence effects at every 
step, and there is now available very extensive experimental data referring to 
situations in which turbulent motions occur to various degrees. Unfortu- 
nately, the turbulent motions themselves have not yet been investigated in any 
detail, the presence of a turbulence being deduced as a rule from its macro- 
scopic manifestations. It is clear that the study of only macroscopic or 
averaged characteristics of the turbulent motion and not the spectrum and 
amplitude of the oscillations will not permit an unequivocal determination 
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of the type of turbulent motion which is responsible for any given macro- 
scopic effect. At the present stage of research into the collective processes 
taking place in plasmas, it is not therefore possible to carry out a complete 
comparison between theory and experiment, and we must limit the discussion 
here to brief survey of the main experimental data. 


(a) Anomalous Diffusion 


The very first experiments with electric arcs in magnetic fields (5) showed 
that in a strong magnetic field fluctuations of fairly large amplitude are 
regularly excited in a plasma, and these considerably increase the effective 
diffusion coefficient. Bohm, in analysing the results of these experiments, 
concluded that the diffusion coefficient decreases only as H~ +, while according 
to classical theory it ought to decrease as H~*. Simon (108, 109) showed, 
however, that these experiments admit of a very much simpler interpretation 
in terms of the so-called “‘ short circuiting ” of the electron and ion currents 
through the electrodes; the results do not then disagree with classical theory. 
However, this explanation did not eliminate the problem, since subsequent 
experiments have revealed enhanced diffusion in a great variety of conditions. 

So far the greatest progress in the understanding of the origin and charac- 
ter of the anomalous diffusion of a plasma across a magnetic field has been 
made in the case of the positive column of a glow discharge. The first 
measurements of the characteristics of the positive column in a magnetic 
field, made by Rokhlin (112), Cummings and Tonks (113), Reikhrudel and 
Spivak (114), and Bickerton and von Engel (115), showed that for small 
magnetic fields the diffusion is classical. Nedospasov (116) arrived at the 
same conclusion by measuring the length of the cathode region of an arc 
discharge, and so also did Vasileva and Granovskii (117), who measured the 
diffusion flux to the wall of the discharge tube directly. However, extending 
the experiments to high magnetic fields, Lehnert (31) observed unexpectedly 
that at a critical value of the magnetic field H, the diffusion flux began to 
increase with the magnetic field, finally reaching a saturation level. This effect 
has been studied in some detail by Lehnert and Hoh (32) and later by other 
authors (34, 35, 118-123). It has been explained theoretically by Kadomtsev 
and Nedospasov (33) on the basis of the current-convective instability. As we 
showed above (Section I.2), this instability drives the discharge into a helical 
form for H > H.. When the magnetic field is further increased the helical 
discharge in turn becomes unstable, and the motion of the plasma becomes 
turbulent, as we described in Section IV.4. The results of a calculation based 
on mixing-length theory are in good agreement with experimental data 
obtained from glow discharges at pressures of the order of 1O~'—1 mm. Hg. 

At smaller values of the neutral gas pressure, when the mean free path 
of the charged particles is of the order of the tube radius, the diffusion 
approximation is not convenient. According to present theoretical ideas, 
instabilities of the drift type must arise in this case. This conclusion is con- 
firmed experimentally. Figure 27 shows, for instance, the dependence of the 
ion current J, to a wall probe on the magnetic field in a mercury vapour 
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discharge at low pressure (p = 3.6.1074 mm. Hg), measured by Artsimovich, 
Nedospasov and Sobolev. The critical magnetic field at which the relationship 
I„(H}) deviates from the H~* law is in satisfactory agreement with the 
theoretical condition for the formation of the ion-sound instability (IV.88). 
The fact that the ion current in the turbulent discharge (H > H.) is of the 
order of magnitude of the current in the absence of a magnetic field also 
agrees with the ion-sound mechanism of the build-up of oscillations. 

The enhancement of the diffusion flux as the magnetic field increases 
above some critical value has been observed also in Penning-type discharges 
(discharges with oscillating electrons), where the longitudinal current is 
absent (126-128). In this case, together with low frequency oscillations, 
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Fic. 27. Dependence of ion current to a wall probe on magnetic field in a mercury 
discharge at 3:6 . 10-4 mm. Hg pressure 


noise at frequencies of the order of 10? cycles/sec is also observed, but this 
obviously is not directly related to the enhanced diffusion of ions. Chen and 
Cooper (128), using probes arranged along one magnetic line of force, have 
shown that the longitudinal phase velocity of the oscillations developing in 
such a turbulent plasma is of the order of 10’ cm/sec, which considerably 
exceeds the velocity of sound c, ~ 10° cm/sec. The frequency of the corre- 
sponding oscillations is of the order of 10’ cycles/sec, and the relative 
amplitude of the density fluctuations is up to 50%. To judge by all these 
data, these oscillations can be considered as drift oscillations, and the corre- 
sponding instability as of the drift-dissipative type. Moreover, it has been 
observed by Simon (129) and Hoh (130) that in a Penning discharge there is an 
additional reason for an instability of the drift type related to the small 
difference between the drift velocities of the electrons and ions in a transverse 
electric field (an analogous effect in the case of a positive column has a 
stabilising influence (90) ). 

The experimental work of Geller (131), who studied the diffusion of a 
plasma in a high frequency discharge also belongs to this group of investiga- 
tions. Geller also observed an increase of diffusion for H > H.. In this case 
oscillations were excited in the radio frequency region. According to Geller, 


126 PLASMA TURBULENCE 


the value of the critical field is inversely proportional to the tube radius a. 
Points corresponding to the observed values of the critical field are shown 
in Fig. 23, where they fall in the region of the ion-sound instability. 

As we mentioned above, one of the earliest experiments on the diffusion 
of a plasma in a magnetic field was performed by Bohm, Burhop and others 
(5) in the study of ion sources. In these experiments the plasma was set up 
in a conducting anode chamber A (Fig. 28) by a beam of primary electrons 
with energy approximately 200 eV, which ionised the neutral gas at a pressure 
in the range 1074-107? mm. Hg. The plasma so formed diffused slowly 
across the magnetic field, at the same time spreading along the lines of force 
and recombining at the chamber ends. The transverse diffusion coefficient 


Fic. 28. Diffusion of ions in an arc discharge 


may be determined from the characteristic length of the radial fall-off of the 
plasma density. Simon (103) showed that because of the “ short circuiting ”’ 
of the electron and ion currents by the conducting chamber, this length is 
determined by the ion diffusion coefficient which exceeds the electron diffusion 
coefficient by two orders of magnitude. When this effect is included, 
the experimental results can be explained on the basis of classical (laminar) 
concepts. The electron diffusion coefficient cannot be estimated from 
measurements of the radial density distribution alone. The higher electron 
current to a positively charged probe described in (5) does, however, support 
the existence of enhanced diffusion of the electrons, most probably related to 
oscillations with frequencies of 20-60 kc/s which are observed in these 
experiments. Subsequent experiments with arcs showed that the diffusion 
of the plasma in these conditions may in fact considerably exceed the classical 
value. For instance, Neidigh and Weaver (132) have found that under certain 
conditions the arc may change over into a regime referred to by these authors 
as “ mode 2” in which a rotating flare is ejected. They assumed that the 
appearance of the flare is related to variations of the neutral gas pressure 
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over the length of the discharge. However, experiments by Zharinov (133- 
135) have shown that the formation of one or several flares is also possible 
when the pressure is uniform. On the other hand, Boeschoten and Schwirzke 
(136) found that in similar conditions the diffusion was classical. 

Experiments with arcs are difficult to discuss theoretically because of the 
processes occurring at the ends. Simon and Guest (137) showed that the 
longitudinal ion current to the ends in such a discharge may give rise to a 
helical current-convective instability. The critical magnetic field correspond- 
ing to the onset of this instability agrees satisfactorily with the experimental 
data. However, this is probably only one of many possible reasons for the 
instability of the plasma and perhaps far from the most important, although 
there is scarcely any reason to doubt that effects of this general type are in fact 
responsible for the instability; these effects were discussed particularly in 
Section IV.3. 

We shall now turn to a group of experiments on the decay of a weakly 
ionised plasma. Excluding from the discussion the work of Bostick and 
Levine (29) referred to earlier, in which the toroidality of the magnetic field 
is important, we find that the papers (138-144) belong to this group. In these 
studies the plasma has a charged particle concentration in the range of 
107-10!° cm~3. If the magnetic field is not very large the measurements of 
the diffusion loss of particles agree with, or at least do not contradict, the 
classical diffusion theory. However, for high values of the magnetic field 
Golant and Zhilinskii have observed an enhanced diffusion. According to their 
results (111), at magnetic fields exceeding 300 Oersted the diffusion coefficient 
of the charged particles in helium can be approximated with an accuracy of 
30-40% by the empirical formula 
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One of the terms represents approximately the classical value D, ~ 10°p/H?. 
The first term, independent of H, may be related to volume elimination 
processes (e.g. volume recombination, electron capture by impurities with 
subsequent recombination, etc.). The second term represents some additional 
diffusion mechanism not directly related to electron—atom collisions and 
depending on the magnetic field. Since these experiments were carried out 
at magnetic fields and pressure close to the theoretically predicted region 
of drift-dissipative instability (Fig. 23), we might suppose that this anomaly is 
related to an instability of this type. 

The next group of experimental papers are concerned with the diffusion 
of a fully ionised plasma. The most comprehensive data on diffusion have 
been obtained on the B-1 and B-3 Stellarators. Ellis et al. (145) found that 
the loss of a fully ionised plasma from the Stellarator exceeds the classical 
losses by three to four orders, so that it cannot be attributed to electron-ion 
collisions. The energising of additional windings, which should stabilise the 
hydromagnetic instability of an ideal (infinitely conducting) plasma, showed 
no effect whatever on the containment time of the plasma, suggesting the 
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existence of a previously unknown instability mechanism. More detailed 
investigations by Stodiek et al. (146) showed that both during the ionisation 
stage and during the decay stage of the plasma, the charged particle loss rate 
agrees with Bohm’s diffusion coefficient, that is, the measured values of the 
diffusion coefficient could be described empirically by the expression 


T 
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where the electron temperature is expressed in eV and the magnetic field in 
kilogauss. 

To verify whether the plasma losses could be ascribed to the ion-sound 
instability (76, 147), Motley (148) studied the decay of a cold plasma in 
Stellarator B-1. In the absence of a longitudinal current the decay of the 
plasma was determined by the volume recombination if the electrons and 
ions, and the decay time was approximately 2 msec. However, when an addi- 
tional electric field (with frequency 20 kc/sec and amplitude 0:01-0:03 V/cm) 
was applied, setting up in the plasma an appreciable longitudinal current, the 
decay rate increased sharply. A critical current level was found, and the 
corresponding longitudinal drift velocity of the electrons turned out to be of 
the order of the sound velocity, 1.e. the velocity of the ions at the electron 
temperature. For small values of the ion temperature this result, it would 
seem, corresponds to the ion-sound instability (75, 76). However, the same 
values of the critical velocity were measured even when the temperature of 
the electrons was of the order of the ion temperature, so that the ion sound 
waves could not propagate because of strong damping at the ions. These 
results agree qualitatively with the theoretical concept of the drift instability 
(Section IV.4). However, it would be premature to describe this as complete 
agreement between theory and experiment, since on the one hand the turbu- 
lent diffusion theory is still in its early stages, and on the other there is as yet 
no complete understanding of the effects of small departures from equilibrium 
of the plasma column (see Section 5(c) below). 

The results obtained on the Stellarator, according to which the diffusion 
of a currentless plasma is comparable to the diffusion of a plasma at rest, 
are confirmed by other experiments with a totally 1onised plasma. Thus, 
D’Angelo and Rynn (149, 152) showed that the diffusion coefficient of potas- 
sium and caesium plasmas in a magnetic field agrees with the classical value 
up to almost 10,000 gauss. Noting that the theoretical diffusion coefficient 
(IV.139) for this case exceeds the classical value by only one order of magni- 
tude, and taking into account the approximations in the theoretical estimate, 
we might think that there is no contradiction between theory and experiment. 
It is interesting to note that in a caesium and potassium plasma device, 
spontaneously excited drift waves have been observed (151,152). The authors 
relate the excitation of the waves to the existence of ion sheaths at the ends 
of the device. 

Golant and Zhilinskii (171, 153) also found that the diffusion of a decaying 
plasma at charged particle concentrations of the order of 10'°-10!* cm7? 
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where electron-ion collisions are dominant, agrees well with the classical 
coefficient of diffusion at magnetic fields up to 1500 gauss. 

Although we have not been able to discuss individually every diffusion 
experiment performed so far, it is clear that the consensus of the results 
is that even in a homogeneous magnetic field and in the absence of external 
beams, the diffusion coefficient both in a weakly ionised and in a fully ionised 
plasma may considerably exceed the classical value obtained from binary 
collisions. In the presence of a longitudinal current the diffusion coefficient 
is considerably higher than in a currentless plasma. In the weakly ionised 
plasma of a glow discharge, it is of the order of magnitude of the diffusion 
coefficient without a magnetic field, and in a fully ionised plasma it approxi- 
mates to Bohm’s coefficient. In the absence of a longitudinal current, the 
coefficient of diffusion is considerably smaller: in a fully ionised plasma at 
moderate magnetic fields, it is of the order of the classical coefficient, but in a 
weakly ionised plasma it may still be considerably higher than the classical 
value. 

All these results either agree satisfactorily with approximate theory 
given in Section [V.4 or at least do not strongly contradict it. In order to 
achieve complete agreement between theory and experiment, however, a 
wide range of investigations is still necessary, both theoretical, to calculate 
rather than only estimate the effective diffusion coefficient, and experimental, 
to reveal the particle loss mechanism in its pure form. 


(b) Turbulent Heating of a Plasma 


By turbulent heating of a plasma—a term having a distinctly applied 
shade—we shall understand the transfer of ordcred energy, either the energy 
of charged particle beams or the energy of discrete oscillations, into energy 
of random motion and ultimately into heat, due to turbulence, i.e. the non- 
linear interaction between oscillations. The possibility of turbulent heating 
of a plasma in this sense was predicted qualitatively by Buneman (154) who 
discussed the excitation of oscillations in the plasma due to a beam instability 
under conditions in which all the electrons move relative to the ions. Such a 
state can be achieved, for instance, by means of a strong external electric 
field leading to electron “runaway ”. Buneman showed by numerical 
calculations that the development of oscillations is accompanied by the 
scattering of the electrons at the inhomogeneities of the electric potential 
and the transformation of the directed motion of the electrons into random 
motion. 

The Maxwellisation of the electron distribution function due to a collective 
interaction has been studied numerically by Dawson (155) in a one-dimen- 
sional model in which a charged plane or surface corresponds to each charged 
particle. Figure 29 shows a picture of the Maxwellisation of two oppositely 
directed electron beams obtained by that author. 

Shapiro (156, 157) studied the initial stages of the randomisation of the 
directed velocity of the electrons by means of moment equations obtained by 


integration of the kinetic equation multiplied by /, v and v’, over the velocity. 
5* 
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He showed that when the electron density in the beam is considerably 
smaller than the plasma density, an initial relatively short stage in which the 
beam spreads in velocitv is followed immediately by a longer stage of re- 
tardation, to which the quasi-linear approximation is applicable. The 
interaction between the beam and the oscillations leads to the formation of a 
plateau in the distribution function. One-third of the energy of the beam, 
which is equally distributed between the kinetic and electrostatic energies of 
the oscillations, is spent in the formation of the plateau. 

Experimental data on the interaction of beams with plasmas have accumu- 
lated over many years. Effects of the collective interaction of an electron 
beam with a plasma were observed by Langmuir (1). The related anoma- 
lously rapid Maxwellisation of the electrons has been given a special name of 
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Fic. 29. Maxwellisation of two interpenetrating beams (a one-dimensional model) 


“ Langmuir’s paradox ’’. The interaction of the beam with the plasma has 
been investigated in greater detail by Merrill and Webb (3). Using probe 
measurements they showed that during the passage of a beam localised 
regions of strong scattering are set up in the plasma. This effect was inter- 
preted by Bohm and Gross (158) as a result of the modulation of the electron 
beam by the oscillating potential of a double layer formed at the plasma 
boundary. Then, as in the operation of a klystron, the initial perturbation 
of the density increases because the faster electrons overtake the slower ones, 
and at some point the electrostatic interaction between the electrons leads to 
their scattering. The scattered electrons which pass through the double layer 
at the boundary may slightly change its potential and give rise to feed-back 
between the oscillations at the entry point and in the region of strong scat- 
tering. Moreover, the oscillations of the double layer may lead to a transfer 
of energy to the plasma electrons, contributing to their Maxwellisation. 
Large amplitude oscillations of the double layer, with a frequency of the order 
of the plasma frequency, were observed experimentally by Gabor (159) 
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which would seem to support this mechanism for the excitation of the oscil- 
lations. However, von Gierke et al. (160), who repeated Gabor’s experiment, 
observed no oscillations at all except some weak low frequency oscillations. 
This result perhaps indicates the possibility of a second mechanism for the 
interaction between the particles and the oscillations, namely a quasi-linear 
growth of the waves in space due to the instability, with the simultaneous 
formation of a plateau in the electron distribution function. The formation 
of such a plateau was demonstrated, for instance, in ref. (161). However, 
in the same paper arguments are also put forward in favour of the klystron 
mechanism of the development of oscillations. All this suggests that in fact 
both instability mechanisms may operate and the interaction between the 
electrons and the waves for small electron density in the beam is described 
satisfactorily by the quasi-linear approximation. 
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Fic. 30. Formation of plateau on distribution function 


The quasi-linear effect of the formation of the plateau at the distribution 
function of the electrons is demonstrated in Fig. 30, which shows the electron 
distribution function before and after interaction with the plasma, measured 
experimentally in (162) (see also (163)). The heating of the plasma electrons 
during the interaction with the beam has been demonstrated in (165). 

In practice, in the case of a bounded inhomogeneous plasma, the nature 
of the oscillations may be very much more complex than is usually assumed 
in theoretical discussions. In particular, when there are widely spaced 
resonance frequencies of the plasma system, strong excitation of separate 
harmonics may take place. Effects of this kind have been observed experi- 
mentally by Looney and Brown (4) in the investigation of high frequency 
plasma oscillations and by Alexeff and Neidigh (164) studying the excitation 
of ion-sound oscillations. Thus the real picture of the excitation of waves 
by powerful beams and the interaction of the waves with the particles 
and with one another may turn out to be much more complex, and to 
require a considerable further development of the theory for its complete 
description. 

The excitation of ion sound by a current in the plasma also belongs to 
the group of phenomena considered in this section. Such oscillations have 
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been observed by Nedospasov (166) in a helium discharge at a fairly low 
neutral gas pressure. Both the order of magnitude of the steady state oscilla- 
tions (~ 107? V) and the spectrum of the oscillations, in which the maximum 
occurred in the region of very low frequencies (œ < Qo), are in very good 
agreement with theoretical ideas (Section IV.2). The excitation of 10n-sound 
oscillations has also been observed by Crawford (167) who studied the noise 
ina mercury discharge. Crawford observed that as the discharge current was 
increased, separate oscillations at discrete frequencies appear first, which 
then merge into a continuous spectrum in the range 10*—10° cycles/sec. The 
phase velocity of these oscillations along the tube was 10’-10? cm/sec. These 
are interpreted by the author as ion-sound waves propagating almost across 
the discharge. The excitation of such waves agrees with the theoretical 
picture given in Section IV.2. By using additional grids, Crawford showed 
that the source of the noise is the cathode (a similar effect is described in ref. 
(166)). However, the picture of the developed oscillations is probably not 
very Sensitive to the form of the source, and therefore a comparison of these 
oscillations with a theoretical spectrum is justified. 

The excitation of ion cyclotron oscillations by a longitudinal electric 
current, a physically similar effect, has been described in ref. (150). The 
excitation of ion oscillations obviously also explains the anomalous heating 
of the ions in Zeta (see following section). As we showed in Section IV.2, the 
ion heating may be related to an increased (“anomalous ’’) resistance of 
the plasma, arising from an additional frictional force between the electrons 
and ions due to the oscillations. An ‘‘ anomalous” resistance has been 
observed on Zeta, and has also been demonstrated experimentally by Thomas- 
sen (169) in an experiment specially set up for this purpose. A similar effect 
has been observed by Adlam and Holmes (170) in the investigation of the 
dynamics of a super-fast pinch. These authors established that for a very 
fast rise of the current (tf ~ 1078 sec) the thickness of the current skin is 
independent of time and about one-and-a-half orders of magnitude greater 
than the expected value c/@,). This result can be interpreted either as the 
result of a considerable decrease in the effective density of the electrons 
participating in the current (such a situation might arise from the “ capture ” 
of most of the electrons by the potential wells @ associated with the ions), 
or because of the suggested impossibility of increasing the directed velocity 
above some small fraction of the thermal velocity. The second possibility 
might arise from the excitation of the ion-sound instability. Adlam and 
Holmes support this point of view. However, the whole effect takes place 
during an extremely small interval of time which is obviously insufficient 
for the development of the ion-sound instability, which would indicate 
rather purely electronic oscillations which may be connected with ion 
perturbations existing in the plasma before the application of the longitudinal 
electric field pulse. In this case only the small proportion of the electrons 
‘ beyond the barrier ” can carry the current. For a complete discussion of 
this extremely interesting and complex problem additional work will be 
required both in theoretical and experimental directions. 
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In an even clearer form, the “‘ anomalous ” resistance effect appears in 
the so-called “ break ” in the current during the last stage of a toroidal 
discharge (see the following section). The capture and acceleration of ions by 
powerful electron beams are associated with a similar effect, observed 
experimentally by Plutto (171). A complete theory of phenomena of this 
kind, most probably related to strong turbulence in the electron gas, has not 
yet been put forward. 

We now go on to consider experiments in which turbulent processes were 
used for heating the plasma. One of these (165) has already been mentioned. 
The anomalous heating of the ions in Zeta (188) was a pleasant surprise. 
One might think that the cyclotron heating of the electrons described 
in reference (173) would also involve the excitation of a large number 
of degrees of freedom. Turbulent heating of the ions in a transverse 
electric field also occurred in Ioffe’s experiments (174), where in addition to 
the heating due to the centrifugal instability of the rotating plasma, heating 
could have occurred due to instability of the plasma filament penetrated by 
the electron beam emitted from the cathode. The latter instability and the 
correponding ion heating has been investigated by Nezlin (175), (176) who 
observed that for a fairly high current density in the beam, when the directed 
flux amounts to approximately 1/6 of the random flux env,, periodic cut-off 
of the beam occurs. The fluctuations of the electric field associated with this 
cutting-off of the beam lead to intense heating of the ions. This effect has not 
yet been interpreted theoretically. On the one hand, it might be considered 
as the result of the development of the two-stream 1on-electron instability 
in some region where the instantaneous density was fairly small due to the 
ion oscillations (which may be excited by the Langmuir oscillations of the 
electron beam (177)). Yet on the other hand, the effect is very similar to the 
electrostatic cutting-off of a compensated electron beam (178), (179). 

Turbulent heating due to oscillations of the magnetosonic type has been 
investigated by Zavoiskti et al. (180-183). In these experiments, heating was 
achieved by superimposing on the main longitudinal magnetic field H, ~ 1000 
gauss a high frequency field with frequency f ~ 10’ cycles and amplitude 
H ~ 500 gauss. As a result of the turbulent heating, the electron tempera- 
ture increased to 500 eV and under certain conditions the ions were heated toa 
temperature of ~ 100 eV. The authors interpret the heating of the plasma 
as the result of the development of a two-stream electron-ion instability. 
The decay of the magnetosonic oscillations may also be important (see 
Section II. 1(b)). 

The experimental results summarized above indicate that the use of 
collective processes for heating a plasma holds considerable promise. 


(c) Toroidal Discharges 


A whole series of turbulent processes has been observed during experi- 
ments with high temperature plasmas in toroidal devices. We have already 
noted above that anomalously fast plasma loss is observed in figure-8 shaped 
Stellarators, the so-called “pump out” effect. The detailed investigations of 
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this effect (146, 148) have shown that the anomalous diffusion is related to 
the longitudinal current, and the effective diffusion coefficient is of the 
order of the Bohm value. These results agree quantitatively with theoretical 
predictions based on the drift instability mechanism (see Section IV, 4). 
However, subsequent experiments made on Stellarator C (184-186) showed 
that the plasma loss may also be related to small departures from equilibrium 
of the plasma column due to inaccuracy of the magnetic field. For very small 
departures from the strict equilibrium conditions, the plasma column may 
be able to preserve an equilibrium by the closure of the confining currents 
through the diaphragm limiting the plasma column. The maximum possible 
value of the current that can be closed in this way is determined by the 
magnitude of the plasma loss. Under ideal equilibrium conditions this loss 
is symmetrical around the small circumference. When the equilibrium is 
upset, a predominantly ion current flows to one half of the diaphragm and an 
electron current to the other half, and the difference between these represents 
the transverse current j; necessary to maintain the equilibrium. For a suf- 
ficiently strong perturbation of the equilibrium condition (which would be 
represented, for example, by imposing a transverse field of only a few tenths 
of a per cent of the longitudinal field), the equilibrium loss current can no 
longer retain the plasma column which starts to move about over a large 
radius, which leads to additional plasma losses. Thus, the loss of particles 
and the deviation from the equilibrium position turn out to be inter-related. 
In the earlier experiments (146, 148) the equilibrium conditions were not 
controlled and some of these results may have to be revised considerably. 

The effect of the transverse displacement of the plasma filament due to 
the perturbation of the equilibrium conditions was observed slightly earlier 
on a Tokomak device (186), where it was shown that to restore equilibrium 
additional current windings can be used setting up a transverse magnetic field. 
A similar application was also described in (184). 

Another turbulent effect is the “cut-off” of the current of runaway 
electrons in an afterglow plasma, also observed on the Stellarator. Dreicer 
(187) has shown that this effect may be related to a build-up of Langmuir 
oscillations occurring due to the appearance of a minimum in the electron 
distribution function in the region between thermal and accelerated electrons. 
This minimum appears because of the faster retardation of slow electrons. 
However, this mechanism does not explain the appearance not of one, but of 
two or more steps as the current falls. 

An unusual variety of different types of turbulent motion of a plasma are 
observed in discharges using moderate longitudinal currents. The most 
complete investigation of these processes has been made on Zeta (188-191). 
Amongst these belongs above all the anomalous heating of the ions. It was 
shown in ref. (189) that the mean energy of the impurity ions can be repre- 


sented in the form E; = Eo + — E,, where m, is the mass of the ion, mp 
m 


. D 
the mass of the deuteron, and Ey and E, are constants, Ey ~ 100 eV, FE, ~ 
2-3 eV. We can definitely conclude that the ions are in principle thermalised 
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and have a temperature of the order of 100 eV, which is several times higher 
than the electron temperature, and that they perform collective oscillations 
with a frequency considerably smaller than the cyclotron frequency, so that 
their total velocity v} = cE/H is independent of the mass, and the energy 
increases linearly with the mass. The heating of the ions may be due either 
to the current instability (for instance the excitation of magnetic sound oscilla- 
tions considered in Section IV.2) or to a secondary effect occurring during 
turbulent convection, arising from the flute or dissipative instabilities. No 
detailed theoretical discussions of this problem is as yet available. 

An effect closely related to the anomalous heating observed on Zeta is the 
excess resistance, which according to the experimental data increases linearly 
with the energy input to the plasma per particle (190). This additional 
resistance is of the same order as the ohmic resistance. Although we have 
noted effects of this type during our theoretical investigation (see Section 
IV. 2), a complete theory of the additional resistance has not yet been de- 
veloped. 

On Zeta and other devices using moderate magnetic fields, an anomalous 
particle loss is observed related to strong oscillations of the electric field (191). 
This loss may be related both to the convective (flute) instability and to the 
drift-dissipative instability. An estimate based on formula (1V.141) shows that 
in conditions in which perturbations greatly extended along the magnetic 
lines of force are admissible, the diffusion coefficient in the conditions ob- 
taining in Zeta may exceed Bohm’s value by two orders of magnitude. This 
value is quite sufficient to account for the observed losses. However, as the 
longitudinal wavelength is reduced the coefficient of diffusion according to 
(1V.143) decreases fairly rapidly. Therefore, in the outer regions of the dis- 
charge, where the effective length of a perturbation along the lines of force, 
being of the order of the pitch of the lines of force, is comparatively small, 
one would expect the primary instability to be the convective instability 
investigated in linear approximation by Suydam (192) It is possible that 
these loss mechanisms come into play and account for the “ magnetic 
number ” effect observed first on Zeta (190) and then on a smaller device 
(193). This effect consists in a prolonged conservation of a constant value 
of the pitch of the lines of force at the periphery of the discharge, the value 
being such that a line of force joins up to itself after one circuit of the torus, 
and abrupt transitions of the pitch from one such value to another as the 
discharge current is increased. One natural explanation of this effect put 
forward by Furth is to assume that on a magnetic surface with closed lines of 
force, where the longitudinal wavelength of the perturbation cannot be very 
large, the principal part in the transport of plasma is played by convective 
perturbations. During the rearrangement of the lines of force in such a 
perturbation the pitch is equalised over the radius, so that a region of constant 
pitch is generated, ultimately encompassing the entire periphery of the dis- 
charge. The enhanced diffusion of the aximuthal magnetic field into the dis- 
charge during its initial stages (194) also supports the existence of such a 
convective motion. 
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The anomalously rapid diffusion of a plasma in a discharge with moderate 
magnetic field leads to the establishment of a force-free configuration (195, 
196). This configuration, which is not quite stationary but includes a small 
amplitude random motion due to the helical instability, can be described 
satisfactorily by an approximation in which it is assumed that because of the 
transverse convective motion of the plasma the apparent transverse conduc- 
tivity vanishes (197). This approximation does not, however, account for the 
very interesting effect of the generation of reverse longitudinal magnetic 
flux at the outside of the plasma filament (195, 196, 198). This effect is hardly 
related to the toroidality, since in somewhat different conditions it is also 
observed on a straight discharge (201). The generation of this reverse field is 
probably related to the finite mixing length during the turbulent convection 
of the plasma (197). However, a quantitative theory of this effect has not yet 
been developed. 

Another turbulent effect also observed first on Zeta (190), and later on 
a device of more modest size (193) is of exceptional interest. We have in 
mind the step-wise decline of the current during the last stage of the discharge. 
As on the Stellarator, there may be several steps in the current, but in con- 
trast to the Stellarator, in this case the decline takes place in the presence of a 
longitudinal electric field which might sustain the discharge current. It 
might seem that in this case also this decline might be connected with the 
excitation of high frequency oscillations due to runaway electrons, the rela- 
tive number of which increases considerably towards the end of the discharge. 
(Such oscillations have been observed experimentally on a straight discharge 
(200).) However, the fact that during the stepped decline of the current the 
magnetic energy is completely transformed into kinetic energy of the electrons, 
approximately | keV of energy being transferred to each electron, would 
rather seem to indicate an electron-ion collective process of the type of an 
anomalous resistance. A theory of this effect is also still absent. 

We must also include in this section the work on the Levitron (199), a 
toroidal device with moderate longitudinal magnetic field and an additional 
current ring which should stabilise the hydrodynamic instabilities of the 
plasma. Oscillations and anomalous plasma loss are, however, still observed. 
These oscillations are obviously of a dissipative nature, related to the finite 
conductivity of the plasma. 


(d) Magnetic Traps 

The turbulent diffusion of a plasma in a trap with magnetic mirrors has 
been observed and studied in detail by Ioffe and others (174, 202-205). In 
these experiments a hot plasma with an ion energy of the order of 1 keV and 
particle density of the order of 10’/cm~ ° was set up by accelerating ions in a 
pulsed electric field which was applied radially between the chamber walls and 
a cold plasma filament situated along the axis of the chamber. After the appli- 
cation of a high voltage pulse, the trap fills with plasma during 10-20 usec. 
No detailed investigation of the motion of the plasma during this heating 
stage has yet been carried out, but we may suppose that the filling of the 
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chamber with plasma and the heating of the ions take place essentially due 
to the centrifugal instability of the rotating plasma. After the high voltage 
is switched off the rotation of the plasma stops and a more quiescent stage of 
turbulent convection due to the flute instability sets in. The presence of such 
a convection is evident chiefly from the anomalously rapid loss of plasma 
from the trap during a period of the order of 1074 secs. Probe measurements 
have shown that the plasma decays essentially at the side walls of the trap, 
whole tubes of plasma lying along the magnetic field being ejected simul- 
taneously towards the side walls. This result indicates that the fundamental 
mechanism of the convection is the flute instability. In a paper by this author 
(16) a semi-quantitative theory of such a convection has been developed, 
based on the mixing length concept. Further investigations of the relation- 
ship between the characteristic scale and amplitude of the turbulent pulsa- 
tions and the distance to the wall, and the dependence of plasma lifetime on its 
density (203), supply additional confirmation of the theoretical model. Final 
proof of the flute instability mechanism was obtained when by the use of 
additional current conductors arranged at the periphery of the trap it was 
possible to stabilise the plasma and achieve prolonged containment (205). 

It is interesting to note that the results on the turbulent convection were 
adequately described on the basis of hydrodynamic concepts without any 
consideration of the effect of the finite Larmor radius of the ions, i.e. of a 
collisionless viscosity, even though in the conditions of the experiments 
described in refs. (202, 203) this effect should be important. It is possible 
that it is smeared out by the strong inhomogeneity of the plasma brought 
about by the turbulent heating mechanism (see Section IV.4(j)), but the 
possibility cannot be excluded that the special features of the method used 
for producing the plasma in fact determine the entire subsequent motion. 
In fact, in the experiments of Ioffe et a/., considerable oscillations of the elec- 
tric field take place during the phase in which the plasma is set up, which 
lead to the loss through the end of ions which are mirrored very near the 
surface of maximum magnetic field H,,. As a result, the plasma boundary 
moves towards the centre of the trap to a point where the magnetic field 
H = H, < H,,. During small oscillations of the potential @ of the plasma 


tubes such that 9 < a(g — 1) all the ions continue to be contained 
e 


S 


by the mirrors, so that nothing prevents the development of the flute insta- 
bility. With other more quiescent methods of heating the plasma, we might 
reach a situation in which the plasma fills the entire trap up to the centre of 
the mirrors where the magnetic field attains its maximum value H = Hm 
In this case even small oscillations of the potential would lead to considerable 
currents to the ends due to the ejection of particles normally mirrored in the 
region of maximum magnetic field. These currents tend to stabilise the 
plasma, and this stabilisation effect is enhanced when there is a cold plasma 
surrounding the hot contained plasma. Most probably it is this effect which 
accounts for the stable containment of the plasma in adiabatic compression 
experiments in traps with magnetic fields (206, 207). 
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Amongst other collective effects observed in magnetic mirror traps, we 
may include the excitation of oscillations at very low density, where the 
transition to a quasi-neutral plasma begins (208, 209), and the excitation 
of ion cyclotron oscillations due to the anisotropy of their distribution 
function (210, 211). No detailed theory of these effects including the non- 
linear terms has been developed so far. 


CONCLUSION 


Summing up, we can state that during the past few years considerable 
progress has been made in developing a theory of plasma turbulence in the 
broad sense of this term, and in the understanding of the nature of the 
collective processes taking place in a plasma in laboratory conditions. A 
considerable number of problems concerning the development and inter- 
action of oscillations in a plasma has been considered theoretically. Together 
with the study of weak turbulence, which is described satisfactorily by the 
kinetic wave equation, methods are being developed for the description of a 
strong turbulence including particularly the weak coupling method. A num- 
ber of specific problems can be considered on a semi-empirical basis by the 
phenomenological introduction of the mixing length concept. At the present 
time, however, important fundamental difficulties are in sight which might 
delay the further development of the theory of plasma turbulence. 

The situation is less satisfactory as far as the explanation of the experi- 
mentally-observed turbulent effects is concerned. Only for a small number of 
experiments is it possible to develop the corresponding theory and to obtain 
a satisfactory description. Yet in the majority of cases, not only has no 
quantitative theory of the observed effects been produced, but often even a 
qualitative interpretation of the nature of the observed effects is missing. 
This state of the matter is due partly to the rather indeterminate nature of the 
experimental data, but also to a considerable degree to the inadequate level 
of development of theoretical concepts and quantitative methods of describ- 
ing collective processes. One must hope that it will be possible in the course 
of the next few years of hectic development of the theory, in conjunction with 
detailed and accurate experiments, to set up a complete picture of the 
turbulence of plasmas. 
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INDEX 
A 


Acoustic waves, see Waves 

Adiabatic compression, 137 

Alfvén waves, see Waves 

Anomalous diffusion, see Diffusion 

Arcs in magnetic fields, 124, 126-127 
“ mode 2,” 126 


B 


Beam—plasma interaction, 130 
instability, see Instability 
Binary collisions, 47, 58 


C 


Caesium plasma, 128 
Cherenkov effect, 47 
Collisions, binary, 47, 58 
Compression, adiabatic, 138 
Conductivity, electrical, of turbulent 
plasma, 119-123 
thermal, 79 
Convection, 7-15 
in positive column, 11-15 
Cyclotron radiation, 48 


D 


Damping, non-linear, 65-69, 76 
Debye length or radius, 47, 65, 80 
number, 45 
Decay instability, 38 
Decaying plasma, 10, 11 
Diffusion, anomalous, 1, 28-30, 48, 78, 
124-129 
turbulent, 3, 106-119, 123 
in velocity space, 48, 58 
in wave number space, 57, 66 
Diffusion coefficient, 29, 78, 96 
Bohm value, 3, 29, 108, 109, 113, 115, 
116, 128, 129, 134, 135 
classical, 115, 127 
turbulent, 107-119 
in weak turbulence, 107 
for ion acoustic instability, 107 


Diffusion coefficient, 
in rarefied currentless plasma, 108- 
110, 119 
in rarefied current carrying plasma, 
110-113 
in very low density plasma, 114 
in dense plasma, 114-116 
effect of finite length, 109-110, 116 
measurements of, 126-129 
in stellarator, 127-128 
in Caesium plasma, 128 
Dielectric permeability, 19, 21, 28, 42, 
47, 65, 69 
Discharge, glow, see Positive column 
H.F., 125-126 
P.I.G., 125 
super-fast pinch, 132 
toroidal, 133-136 
Dispersion relations, 74, 76, 80, 83, 92, 
94, 95, 97 
decay type, 34, 35, 119 
non-decay type, 34, 35, 62, 76 
Double layer 130 
Drift instability, see Instability 
Larmor, 79 
waves, 79—84 


E 


Electron Larmor radius, 98, 99 
Energy exchange time, 41 
Enhanced diffusion, see Diffusion 
Excitation, hard, 5, 6 

soft, 5, 6 


F 


Flow, non-linear, 103-106 
Flute instability, see Instability 
Four-wave processes, 35, 68 


G 


Glow discharge, see Positive column 
Green’s function, 50, 57 
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148 


H 
Hall current, 120 
Hard excitation, 5, 6 
H.F. discharge, 125-126 


Instability, absolute, 24 

beam plasma, 15, 16, 130 

centrifugal, of rotating plasma, 133, 
137 

convective, 24 

current-convective, 12, 86, 99, 102- 
103, 104, 124, 127 

drift, 4, 78-90, 102, 125, 128, 134 

drift-dissipative, 4, 79, 96-102, 115, 
125, 127, 135 

flute, 90-93, 105-106, 135, 137 

ion acoustic, 68, 94-95, 97-100, 107, 
125, 128, 131-132 

ion cyclotron, 74, 95-96 

of inhomogeneous plasma, 78-103 

of weakly ionised plasma, 96-102 

two-stream (electron-ion), 133 

Ion acoustic waves, see Waves 

cyclotron waves, see Waves 

Landau damping, 81, 85, 88, 89 

Larmor radius, 78, 81-92, 105-106, 
119, 137 

sound waves, see Waves 


J 
Jet, turbulent, 59-61, 118 


K 
Kinetic wave equation, 2, 31-55, 42-44, 
139 


Kolmogorov-Obukhov spectrum, 3, 51, 
61-62 


L 


Landau collision term, 47 
damping, inverse, 16 
ion, 81, 85, 88, 89 
linear, 41, 44, 84, 86 
non-linear, 41, 68, 76 
Langmuir paradox, 1, 130 
waves, see Waves 
Larmor drift, 79 
radius, electron, 98, 99 
ion, 78, 81-92, 105-106, 119, 137 
Levitron, 136 


INDEX 


M 


Magnetic mirror trap, 3, 136-138 

‘* Magnetic number ” effect, 135 

Magneto-sonic waves, see Waves 

Maxwell stress tensor, 28, 29 

Mixing length, 3, 59-61, 116-118, 124, 
137, 139 


N 


Navier-Stokes equation, 51 
Non-linear damping, 65-69, 76 
flow, 103-106 
Landau damping, 41, 68, 76 


P 


Penning (P.I.G.) discharge, 125 

Permeability, dielectric, 19, 21, 28, 42, 
47, 65, 69 

Perturbation theory, 51 

Plateau in distribution function, 17, 72, 
75-77, 131 

Positive column, 3, 
124-125 

Potassium plasma, 128 

Pump-out, 134 


11-15, 116-118, 


Q 
Quasi-linear approximation, 2, 5-30, 72, 
75, 131 


R 


Random phase approximation, 33, 46 
Resistance, anomalous, 132, 135, 136 
Resistivity, anomalous, 48, 77 

see also Conductivity, electrical 
Reynolds’ stress, 59 
Runaway, of electrons, 129, 134, 136 


S 


Shear, of magnetic field, 89-90 

Shock, collisionless, 2 

Soft excitation, 5, 6 

Sound waves, see Waves 

Spectrum of turbulence, for ion acoustic 

waves, 68-72 

for wind waves, 62-63 
Kolmogorov-Obukhov, 3, 51, 61-62 

Stellarator, 127-128, 133-136 

Strong turbulence, 2, 3, 49-63, 116 


INDEX 


T 


Thermal conductivity, 79 
fluctuations, 44—46 
noise level, 48 
Three-wave processes, 34 
Tokomak, 134 
Turbulence, fluid, 51 
strong, 2, 3, 49-63, 116 
weak, 2, 7, 18, 31-48, 49, 139 
Turbulent, jet, 59-61 


V 


Van de Pol method, 6 
Viscosity, 51, 61, 79 
collisionless, 105, 137 


W 


W.K.B. approximation, 23, 27 
Waves, Alfvén, 2, 39-41, 76, 83 
accelerated, 83 
decelerated, 83 
drift, 79-84, 125 
ion acoustic, 35-39, 47, 68-73, 78, 81, 
100, 131, 132 
accelerated, 81 
decelerated, 81 
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Waves, Alfvén 
ion cyclotron, 73-75, 132, 138 
ion sound, see Waves, ion acoustic 
Waves, Langmuir, 2, 3, 15-18, 26, 28, 
35-39, 42-44, 47, 51, 65-68, 134 
magneto-sonic, 2, 39-41, 75-78, 133 
sound, see Waves, 10n acoustic 
Wave equation, kinetic, see Kinetic 
wave equation 
“ Wave number,” 67-71 
Wave packets, 54-57 
Wave propagation, 2, 7, 18, 31-49 
in inhomogeneous plasma, 18-23 
Wave-particle interaction, 46-48 
adiabatic, 25-30, 48 
resonance, 17, 25-30, 48 
Wave-wave interaction, adiabatic, 54, 55 
resonance, 54, 55, 57 
Weak coupling approximation, 3, 49-58, 
116 
Weak turbulence, 2, 7, 18, 31-48, 49, 
139 
Whistlers, 76 
Wind waves, 62-63 


Zeta, 132-136 , 


